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Abstract: We propose a general correspondence between gravity and spin models, 
inspired by the well-known IR equivalence between lattice gauge theories and the 
spin models. This suggests a connection between continuous type Hawking-phase 



transitions in gravity and the continuous order-disorder transitions in ferromagnets. 
The black-hole phase corresponds to the ordered and the graviton gas corresponds 
to the disordered phases respectively A simple set-up based on Einstein-dilaton 



gravity indicates that the vicinity of the phase transition is governed by a linear- 
dilaton CFT. Employing this CFT we calculate scaling of observables near T c , and 
obtain mean-field scaling in a semi-classical approximation. In case of the XY model 
the Goldstone mode is identified with the zero mode of the NS-NS two-form. We 
show that the second speed of sound vanishes at the transition also with the mean 
field exponent. 
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1. Introduction 

There has been great progress recently in applications of holography [p], |2], [3] to 
condensed matter systems such as superconductors following the pioneering works 
of P| and [[5|. These authors managed to find a simple gravitational background 
in Einstein-Maxwell gravity coupled to a complex scalar field where a second order 
normai-to-superfluid type transition occurs at finite temperature. The basic interest 
behind application of holographic ideas to condensed matter theory (CMT) lies in the 
hope that the strongly correlated condensed matter systems may secretly possess a 
gravitational description. Indeed, computations of certain observables in the gravity 
picture, such as conductivity provides supporting evidence, see @@|§ for reviews. 
It is a considerable possibility that || the underlying dynamics behind the phase 
transition in high T c superconducting materials is a strongly coupled quantum phase 
transition at zero T. Then the hope is that, a dual gravity description of the strongly 
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coupled field theory around this critical point may also shed light over the finite 
temperature transition in the quantum critical region. 

On the other hand, there are several issues of fundamental importance in the 
proposed gravity-CMT models, such as the role of the large N limit and the notion 
of weak-strong duality, that are not entirely clarified. We have a much better un- 
derstanding in the holographic constructions of gauge theories, thanks to the basic 
example [I], Q [| of the M = 4 super Yang-Mills theory where the D3 brane picture 
provide the link between the gauge side and the gravity side. Such a "top-bottom" 
approach is missing in the gravity-CMT models. 

In this work, we entertain the possibility that such a link may be established 
under certain assumptions, at least for certain simple condensed matter systems, i.e. 
spin models, by analogy with the better understood gauge-gravity case. 

The building blocks of such a connection are already present in the well-known 
literature. First of all, we recall the famous equivalence between lattice gauge theories 
(LGT) and spin-models (SpM)[[IT], |TTJ: Integrating out the gauge invariant degrees 
of freedom in the partition function of a LGT with gauge group G, one arrives at 
an effective action for the lowest lying mode, namely the Polyakov loop P. This 
effective action is invariant under the leftover center symmetry C = C enter (G) of 
the original gauge invariance. Identifying the Polyakov loop P with a spin field s, one 
then obtains the partition function of a spin-model with the global spin invariance 
C. Using this equivalence between lattice gauge theories and spin-models Polyakov 
and Susskind were able to show the existence of confinement-deconfinement phase 
transition on the lattice, long time ago. Based on these works, than Svetitsky and 
Yaffe [0 further proposed that, if continuous critical phenomena prevails in the 
continuum limit of a certain lattice gauge theory, then it should fall in the same 
universality class as the corresponding spin-model. 

It is interesting to employ the same idea in the opposite direction in order to 
study a spin-model that is strongly coupled at criticality. In particular, one would 
like to compute the critical exponents, the transition temperature T c , certain thermo- 
dynamic functions etc., by analytic methods. If one is lucky enough to find a gauge- 
theory that corresponds to the spin- mo del under the aforementioned equivalence, 
then one may be able to study the strongly coupled phenomena by the gauge-gravity 
correspondence. 

One purpose of this paper is to emphasize that this chain of dualities may pro- 
vide a well-defined setting in understanding fundamental issues in the gravity-CMT 
correspondence. In particular, if one can figure out the relevant D-brane configu- 
ration that describes the gauge theory which arises in the continuum limit of the 
LGT under question, then one may be able to take the decoupling limit and obtain 
a gravity description of the LGT — and of the equivalent spin model — around criti- 
cality. Despite being abstract, in principle this provides a top-bottom approach to 
the problem. In particular, such an approach would hopefully provide a microscopic 
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description that is long sought for in holographic applications to CMT. 

Another purpose of this paper is to provide a concrete realization of these ideas in 
a simple setting. For this purpose we consider SU(N) gauge theory in d- dimensions 
(with possible adjoint matter) in the strict N — y oo limit. In this limit the center 
C becomes U(l) . We imagine that the adjoint matter is arranged such that the 
deconfmement transition of the gauge theory is of continuous type. This transition 
is then in the same universality class with the order- disorder transition in the cor- 
responding U(l) rotor model in d — 1 dimensions — that is sometimes called the XY 
model. The XY- models — and their 0(n) generalizations — provide canonical exam- 
ples of superfluidity that arises as spontaneous breaking of the global U(l) symmetry 
in a continuous phase transition. 

To realize this phenomenon in the dual gravity setting we consider the NS-NS 
sector of non-critical string theory in d + 1 dimensions with Euclidean time direction 
x° compactified. It was shown in []l5j| that this theory in the two-derivative sector 



exhibits a continuous Hawking-Page transition at some finite temperature T c . The 
background is of the type AdS d+1 near the boundary and linear-dilaton in the deep- 



interior. Building upon the ideas in [13|, we argue that the U(l) symmetry (in the 
strict N — y oo limit) corresponds to the shift symmetry J M B — y j M B + const where 
B is the NS-NS two-form field and M is the (r, x°) subspace of the background 
geometry. The only objects that are charged under this symmetry are string states 
winding the time circle. In the thermal gas phase these states have infinite energy and 
cannot be excited, hence the symmetry is unbroken and this phase corresponds to the 
normal phase of the spin system. In the black- hole phase on the other-hand they have 
finite energy (with an appropriate regularization) and the black-hole corresponds to 
the superfluid phase. 



It was further observed in [|15] that the geometry becomes exactly linear-dilaton 
in the transition region. Therefore, we argue that the transition region of the XY 
model is governed by the linear-dilaton CFT on the string side. Although in gen- 
eral the a' corrections can not be ignored in the type of backgrounds that we will 
consider in this paper 2 , one can still perform calculations in the critical regime, pre- 
cisely because the linear-dilaton background is known to be an a'-exact background 
in non-critical string theory [0. In particular the calculations that involve probe 
strings can be performed by employing the exact CFT description of the linear- 
dilaton background, (in the limit g s — y 0). 

The spin operator s(x) is related to a fundamental string that wraps the time- 



1 This idea in the AdS/CFT context was considered before]13|, see also jlj] and |12| for earlier 
discussions 

2 We recall that in the case of Af = 4 sYM theory the a' corrections can be ignored both for the 
bulk and the string computations at strong 't Hooft coupling A. In the theories we consider here we 
do not have a similar modulus that serves as a parameter to suppress the a' corrections. Generally, 
the string scale and the scale of the background geometry may be of the same order. 
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circle and connected to the boundary at point x. Consequently one can compute 
correlation functions of the operator s by studying the string propagation in the 
linear-dilaton CFT in the (single) winding sector. We perform such calculations in a 
semi- classical limit where we only take into account the contribution of the lowest- 
lying string states. It is shown that in this approximation, one obtains mean-field 
scaling near T c . We find that the "magnetization" behaves as 

M ~(T-T C )2, as, T^T C . 

A similar calculation with string propagation connecting the points x and y on the 
(spatial) boundary corresponds to the spin two-point function. We show that the 
expected behavior of the spin-system arises in the large \x — y\ limit near T c indeed 
arises from this calculation in a non-trivial manner. In particular, in order to show 
that the correlation length £ diverges at T c , one has to identify the transition with the 
Hagedorn temperature where the lowest lying single-winding mode becomes massless 
[ pp . With such an identification one indeed finds the expected behavior 

£~(T-T C )"5, as, T^T C , 

again in a semi-classical approximation. 

One can also study scaling of the speed of second sound that is associated with 
the Goldstone mode in the superfluid phase. This mode is identified with fluctuations 
of the zero- mode of the NS two- form field B. We find that the speed of sound indeed 
vanishes at T c precisely with the expected mean-field scaling, 

C 2~(T-T C ), asT^T c 

in a second order Hawking-Page transition. We also argue that this finding is not 
altered by possible a' corrections. 

The identification of spin operators with the F-strings suggest a similar identi- 
fication between the vortex configurations — that play an important role in the 2D 
XY model — with D-strings in the gravity dual. We study correlation functions of 
such D-string configurations and find that they exhibit the expected behavior in the 
spin-model. 

The paper is organized as follows. In the next section, we review basic ideas 
in the past literature which indicate a general duality between spin-systems and 
gravity. We first focus on the case of SU(N) in the N — > oo limit and postpone the 
general discussion to section 6. Section 3 reviews the Einstein-dilaton system that was 



studied in [15]. In section 4 we argue that the IR limit of the model is described by a 
linear-dilaton CFT and review basic features of such CFTs. Section 5 contains main 
technical results of this paper. We first review the basic statistical mechanics results 
that are relevant in what follows. Then we propose the precise identification between 
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the F-string configurations and the spin correlation functions. We calculate the one- 
point and two-point functions near criticality in the semi-classical approximation 
making use of the linear-dilaton CFT. Finally we present calculations related to 
vortex configurations. In section 6 we take a first step in formulating a gravity spin- 
model duality in general. In the last section we discuss various issues and possible 
future directions of research. 

Several appendices detail our presentation. In appendix A, we review the sim- 
plest example od the equivalence between lattice gauge theories and the spin-systems. 
In appendix B, we review the connection between non-critical string theory and the 
linear-dilaton background. Appendix C provides some basic background material in 
statistical mechanics of the XY models for the unfamiliar reader. Finally, Appendices 
D and E contain details of our calculations in section 5. 



2. Gravity - spin model duality 

Our goal in this section is to propose a particular approach to the gravity-CMT 
correspondence that relates the spin-models in CMT to gravity by a two step pro- 
cedure: The first step is to employ a well-known equivalence between spin-models 



and lattice gauge theories |T0|, | I J followed by a second step that is to utilize the 
gauge-gravity duality to relate the (continuum limit) of the lattice gauge theory to 
a dual gravitational background. 

2.1 Correspondence between gauge theories and spin systems 

Existence of the confinement-deconfinement phase transition in lattice gauge theories 
at strong coupling is rigorously proved [0, [TT]] long time ago. The proof is based on 
an equivalence between lattice gauge theories (LGT) and spin systems with nearest- 



neighbor ferromagnetic interactions, [^, O, |T^]. In the original papers of Polyakov 
and Susskind, this equivalence was established for the cases of U(l) and SU(2) 
gauge theories. Subsequently it was generalized to general Lie groups 3 . We shall 
refer to this equivalence as the LGT-spin model equivalence. We review how the spin 
systems arise from the lattice gauge theories in the Hamiltonian formalism, and in 
the simplest case of U(l) gauge group in Appendix A. 

This equivalence has profound implications in the continuum limit: As ar- 



gued and verified with various examples by Svetitsky and Yafle ||12| , the critical 
phenomena — if exists — in the continuum limit of the LGT, should be in the same 
universality class with the corresponding spin model. Therefore, a continuous order- 
disorder type transition in a d — 1 dimensional spin-model with global symmetry group 



3 See Em for a recent presentation of how the map works in a general case. 
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C is directly related to a continuous type confinement- deconfinement transition of the 
gauge theory with gauge group G where C = C 'enter '(G). 4 

Let us briefly review the argument of fL2f . The basic observation is that the 
magnetic fluctuations are always gapped both in the high and the low T limit of 
the lattice gauge theory. Therefore, they are expected to be gapped for any T on 
a trajectory crossing the phase boundary in figure |I|. This means that the mag- 
netic fluctuations should not play an essential role at criticality in the vicinity of 
a continuous confinement- deconfinement transition. Integrating these short-range 
fluctuations, one indeed obtains an effective theory that only involves the Polyakov 
loops, which in turn can be mapped on a spin model. Therefore the critical phenom- 
ena, e.g. the critical exponents etc. of the lattice gauge theory around a continuous 
transition should be governed by the corresponding spin model. 



00 








9, 



00 



Figure 1: Typical phase diagram of a lattice gauge theory with non-trivial center. Low 
T phase is confining with vanishing expectation value for the Polyakov loop P and high 
T phase is de-confined. We assume that (at least a portion) of the phase boundary that 
separates these phases is of second or higher order. Then the critical phenomena around 
the phase boundary is determined by the corresponding spin-model. 

The magnetic sector is gapped at low T by assumption. We assume that the 
(bare) coupling constant is large enough (see figure P so that the low T theory is 
confined. The argument at high T is as follows. In the Lagrangian formulation of 
the LGT one can take the action to be, 



A lgt = J2 Re \ fo J2 Trf/ ^ + & J2 TrU ^ \ 

r (. i ij ) 



UfjfAV U r,i±U r-\-[i,,vU ft+j) ,,. h 



u,fi r,u 

(2.1) 



4 0f course, not all of the spin-models exhibit continuous transitions. See Q for a list of 
examples. 
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where r labels sites on the square lattice, Uf,\w are the product of link variables on a 
plaquette with corner f. The first term in the action above corresponds to the electric 
contribution and the second to the magnetic. The electric and magnetic coupling 
constants are related to the bare coupling constant of the LGT and the temperature 
as follows: 

2 4 



9 2 




a "v^ T= x /£— (2.2) 



where a is the lattice spacing and N t is the number of the lattice sites in the Euclidean 
time direction. One observes that, at fixed coupling, f3 t ~ T. Then, for sufficiently 
high T, only configurations with vanishing electric flux contributes in the partition 
function. This reduces the system to static configurations at high T, hence the 
theory can be thought of a d — 1 dimensional LGT at zero T, with a coupling 
constant g\_ 1 = g^T. Any such LGT with a non-trivial center is confined and exhibits 
magnetic screening at strong coupling fll2"| . As mentioned before, the equivalence to 
the spin-models can be shown exactly at strong coupling, [TO, ITT . 



Svetitsky and Yaffe |T2| were able to make reliable predictions concerning the 
critical phenomena of a wide range of lattice gauge theories making use of this con- 
nection and the well-known results on the critical phenomena of the corresponding 
spin-models. 

First of all, they correctly predicted that the 2nd order transition in the pure 
SU(2) theory in 4D is in the same universality class with the 3D Ising model (see e.g. 



[ [TTJ| ] and references therein). As another check of these arguments [|I2|] presents the 
example of SU(N) theory for d — 1 = 2, N > 4 where the dual spin model is again 
Z^ symmetric and exhibit a BKT type continuous transition. In this case, it was 
argued that for large N , the theory approximates that of a U(l) LGT in 2+1 and the 
corresponding spin-model should be the XY-model in 2D. It was explicitly checked in 
[ [12| that, for the U(l) LGT the critical phenomena is in the same universality class as 
that of the 2D XY model. More generally, if the d — 2 + 1 SU(N) gauge theory — or 
a suitable deformation with additional adjoint matter — involves continuous critical 
phenomena than it should be the BKT type. 

A particularly interesting case concerns SU(N) gauge theory in d— 1 > 2 spatial 
dimensions with N > 4 (that includes the large N) where the dual spin model is Z^ 
symmetric. We then consider the large N limit that is most relevant for the gauge- 
gravity duality. It is reasonable to believe that in the strict N — > oo limit (with or 



without adjoint matter), the center Zn is promoted to U(l). See [|T3| for an argument 
in favor of this, in the case of J\f = 4 SYM at strong-coupling 5 . Another indication 
that this happens in Z^ invariant LGT at d = 2 is explained in |]12||. Therefore, 



5 See however p(| which shows that the U(l) symmetry is expected to arise only in the strict 
N — > oo limit. 



by the universality arguments above, if there exist a continuous phase transition it 
should be governed by a U(l) invariant spin model 6 . 

We review how the equivalence of LGTs and spin-systems work at strong coupling 
in Appendix |A] for the unfamiliar reader. Here we shall mention two salient features. 

• The temperature of the spin-system is inversely related to temperature in the 
original gauge theory: 

T s ~ T~\ (2.3) 

Consequently, the low temperature (confined) phase of the gauge theory corre- 
sponds to the high temperature (disordered) phase of the ferromagnet, whereas 
the high temperature (de-confined) phase of the gauge theory corresponds to the 
low temperature (ordered) phase of the ferromagnet. 7 

• Quite generally, the LGT-spin model equivalence can be generalized to incor- 
porate (adjoint) matter. This is mainly because the basic ingredient in the 
calculation i.e. the center symmetry of the LGT remains intact upon addition 
of adjoint matter. See [E3| for a related recent discussion. 



2.2 Holographic superfluidity 

Here and until section 6, we specify to the particular case of U(l) invariant spin- 
models. Continuous critical phenomena in such models include the interesting case 
of superfluidity, that requires spontaneous breaking of the global U(l) symmetry. As 
reviewed above this transition is directly connected to the confinement- deconfinement 
transition in the gauge theory. In the original derivation of |II| [|ll|] , the de-confined 



phase of the U(l) invariant LGT was understood as an ordered phase of the U(l) 
spin model. This is clear from the discussion of section |2|, as the center of U(l) is 
17(1) itself. 

Instead, here we shall adopt an alternative approach where the U(l) factor arises 
from the large N limit of an SU(N) gauge theory (pure or with adjoint matter). In 



6 One can ask whether there is any evidence, for or against criticality at N — oo. There are 
two independent arguments that argue for a second order transition |2l] [Q in the case of pure 
YM in 3+1 dimensions. On the other hand, there is the usual argument against a continuous 
transition at large N that claims, since the number of degrees of freedom in the system changes 
from 0(1) to 0(N 2 ) in a confinement-deconfinement transition, latent heat should be finite. In 
Jig ] we presented a counter-example to this reasoning, albeit in a gravitational setting: although 
the degrees of freedom change abruptly as the graviton gas deconfines in the black-hole phase, 
the entropy difference may vanish at the transition. See |14j for other examples of second order 
transitions at large N. Finally, even if the transition is first-order for pure YM, the situation may 
change when one adds adjoint matter. 

7 In d = 2, IR divergence of the spin waves prevent ordinary long range order. Instead, a 
topological long-range order in terms of the vortex-anti- vortex pairs arises p3[|24]. The gauge 
theory partition function is capable of describing the vortex configurations [n2l . 
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this case the deconfinement transition can be understood in the gravity dual as a 



Hawking-Page transition by a generalization of the arguments in |L3| . Assuming that 
the following assumptions hold, 

• There exists a suitable SU(N) lattice gauge theory with coupling to adjoint 
matter chosen such that, at large N it flows to an IR fixed point with a contin- 
uous confinement- deconfinement transition, 

• Gauge-gravity correspondence holds and maps this to a Hawking-Page type 
transition, 

then one should be able to map the normal-to-superfluid transition in the XY model 
to a continuous Hawking-Page type transition on the gravity side. 

The U(l) symmetry of the spin-model follows on the GR side from the shift 
symmetry ip — > ip + const, where ip is the flux of the B-field 



iP = B = const. (2.4) 

J M 

on the subspace M of the BH geometry that is spanned by the coordinates r and xq. 
In this paper we consider gravitational set-ups where the B-field is either constant or 
pure gauge B = dt; so that it does not back- react on the solution with H = dB = 0. 
Of course such a B-field has no visible effect on the gravitational solution and in the 
second case it can be removed by a gauge transformation. This ceases to be the case 
in presence of objects that are charged under this shift symmetry. 

In the classical approximation where one keeps only the low-lying gravity fields, 
there are no bulk fields that carry the extra B charge. However, strings that wind 
around the time-circle couple to the B-field through the term tip, thus they are 
charged under the shift symmetry with the identification ip ~ ip + 2tc. We shall 
denote this topological U(l) symmetry as U(1)b 8 Therefore a non-vanishing string 
one-point function signals a breakdown of the U(1)b symmetry. On the spin-model 
side this corresponds to an order- disorder transition upon identification of the U(1)b 
symmetry with the U(l) spin symmetry of the spin-model. Below, we would like to 
review these ideas in more detail. 

2.3 Spontaneous breaking of U(1)b, the Goldstone mode and the second 
speed of sound 

For simplicity, let us consider the (critical or non-critical) bosonic string theory on 
a background with U(l) x E(d — 1) isometry where the U(l) corresponds to the 



8 This symmetry should be broken down to Z^ for finite N by quantum effects, see EQ]. However 
we only consider the N — ¥ oo limit in this paper. 
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temporal S 1 , and E(d — 1) to translations and rotations on the spatial part. The 
general background with these symmetries is of the form 9 

ds 2 = A(r, tydxl+B(r, tydK 2 +C(r, tydr 2 +D(r, fi)dfi; $ = $(r,fi) H = dB = 0, 

(2.5) 
where Xq ~ Xq + 1/T, i^ is the d — 1 dimensional transverse part, and f2 is some 
internal compact manifold. There can be additional bulk fields but we are only 
interested in the NS-NS sector. 

Most of the following traces the arguments in |13| . The order parameter for the 
transition is the vev of the Polyakov loop, (P[C]), where C is a loop isomorphic to 
the time-circle. This maps to the expectation value of the F-string path integral, 

(P[C\) ex (W f )sg (2.6) 

where Wp denotes the F-string path integral over all of the string configurations 
with the boundary ending on C, and the final averaging is path integral over the 
super-gravity fields that couple to the string. The string path integral is 

W F = fvX^Vh ab e -/(G-HB+M (a) ) j (2.7) 

where PS 2 ' 1 is the Ricci scalar on the sub-manifold M that the F-string wraps and 
X^ denotes the matter fields. We also use the short-hand notation 

G = Vdet h ab h ab d a X' M d b X' / G ia/ , B = Vdet h ah e ab d a X»d b X v B^ v . (2.8) 

One has to make sure that W is finite by an appropriate regularization of infinite 
volume of the space time 10 and factoring out diffeo- Weyl gauge volume a la Faddeev- 
Popov. 

In the original discussion of |L3] W is dominated by the classical saddles that 



minimize the action in fl2.7p . The boundary condition for these classical strings is 
such that at r = 0, X M (a, r) ends on the temporal circle xq, some point x in K and 
at the cut-off of the radial coordinate r = e. 

The string path integral is dominated by classical saddles when £/£ s 3> 1 where £ 
is the typical curvature of the target space and £ s is the string length. In the original 
AdS/CFT correspondence this ratio is proportional to the t 'Hooft coupling of the 
dual M = 4 SYM theory, £/£ s ex A* and indeed the classical strings dominate in the 
limit of strong interactions. In the general case here one has to consider the full path 
integral. 



9 In order to distinguish the dilaton and the scalar field that appears in the Einstein- frame 
potential, which is related to the dilaton by some rescaling, we denote the former (dilaton itself) 
by $ and the latter (rescaled dilaton) by $. 

10 A cut-off in r that we call e close to the boundary would suffice for the sake of the discussion 
here. We elaborate on this regularization in appendix p.l|. 
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The vev of P[C] is given by the path integral of Wjr over the super-gravity fields 
that couple the F-string, weighted by the SG action. The non-trivial SG fields are 
the space-time metric G^, the B- field B^ v and the dilaton $. Thus one has, 

(P[C\) oc / VG^ VB^u V§ e- A °°W F , (2.9) 

where A sg is the gravity action. As we are interested in the large N limit of the 
dual field theory, we can send the string coupling g a — )■ and the SG path integral is 
dominated by the classical saddles of A sg . For given asymptotic boundary conditions 
of G, $ and B, the saddles of interest involve only two type of solutions, the thermal 
gas (TG) and the black-hole (BH). At an arbitrary temperature T (that partially 
determines the asymptotic boundary condition for G), only one of these saddles will 
dominate the SG path integral as a result of the classical limit g s — > 0. 

Let us assume that TG dominates at T < T c and BH dominates at T > T c . 
Let us also assume that the BH solution only exists above a certain temperature 
T min . Backgrounds that exhibit confinement generically satisfy T c > T min E5I, BB]. 



As explained in detail in [fL5|] and reviewed in the next section, only in the case 
Tmin = T c the transition is second or higher order. 

On the TG phase, the classical world-sheet M has infinite area. Therefore the 



string path- integral WV, hence {P[C}) in (|2.9| ) vanishes. One concludes that the TG 



solution is U(1)b symmetric and the center in the dual gauge theory is unbroken. 
This means that the dual spin-model is in the normal (disordered) phase. This 
is precisely as one expects from the behavior of the dual spin model in the high 
temperature phase, recalling that the temperature of the spin model is inversely 
proportional to the temperature on the gravity side T s oc T _1 . 

On the BH solution T > T m i n however, the classical string saddle M has finite 
area and one has to evaluate Q2.9|) carefully. One has to include all of the configu- 
rations over the classical fields G, B and $ with the same on-shell value of the SG 
action. 



The path integral over G and $ in (|2.9|) is replaced by the classical solution (|2 



that is a BH in this case. Sum over these saddles include the following important 
contribution from the B-field. In the black-hole case the sub-manifold M has finite 
area 11 and the B-field has a flux if) = J M B. if) in (|2.7|) has angular nature because 
it appears with a factor of i and it can attain any value in the range if) ~ if) + 2n. 
This identification yields the U(1)b invariance 12 . The sum over classical saddles then 
should include various different values of ip. As dB = all different values of if) yield 
the same on-shell gravity action. 



11 The divergence near boundary is regularized in the familiar way, cf. appendix D.l. 

12 In the critical IIB theory this identification arises as a result of discrete gauge transformations 



that shift the value of ip by a multiple of 2ir |13] . 
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We can thus write, 

P[C\ oc j V^e~ s ^ ] fvX fl Vh ab e^ e~ Jm(G+*h (2) ) , (2.10) 

where S sg now is evaluated on the saddle solution and is only a functional of ip. On 
the other hand the ip path integral includes the classical saddle ip = const and the 
fluctuations 5ip(K) around it. 

When K is non-compact and dim(K) > 2, then the fluctuations 5ip(K) viewed 
as a massless bosonic field on K has long-range order, hence ip should condense 13 . 
Thus, on the black-hole solution the U(1)b symmetry breaks down 14 . This happens 
exactly at the point where the black-hole forms, right above T min . As a result, the 
fluctuation Sip in (\2.10j) becomes a Goldstone mode on the transverse space K . 



Considering the wave equation for 5ip one expects to find, 

u; 2 = 4(T)q 2 + 0(q 4 ), (2.11) 

where c^ is the speed of sound of ip and there is no mass term for ip for T > T c . It 
is well-known that (see appendix for a review, and section |5.6| for a holographic 



derivation in gravity), the speed of sound c^ of the Goldstone mode vanishes con- 
tinuously as one approaches the transition temperature T c from above, only if the 
transition is of continuous type. This is exactly what happens in super-fluidity, where 
the " second speed of sound' , i.e. the speed of sound associated with the entropy waves 
vanish as one approaches T c of the XY model from below (recall that temperature in 
gravity and in the XY model are inversely related). In order to mimic this property 
of the spin model, we should require that the Hawking-Page transition in gravity is 



of continuous type, hence T c = T min [0. In section [5]6] we show by an explicit grav- 
ity calculation that indeed the second sound vanishes with the expected mean-field 
exponents. 

Our conclusion is: whenever a second order (or higher order) Hawking phase 
transition occurs in the gravitational background, it is natural to associate it with 
super-fluidity Here the thermal gas phase is dual to the normal phase of the system, 
and the black-hole phase is dual to the super-fluid. The "first speed of sound" i.e. 
the sound of the density waves is associated with the graviton fluctuations (that 
we are considered in [fOJ), and the "second speed of sound" is associated with the 



fluctuations of the B-field that we consider in section 5.6. 



13 The situation at dim(K) = 2 exactly parallels the analogous situation in the 2D dual field 
theory, where IR divergences kill long-range order. 

14 As a technical aside, in the computation above, one should check that the dilaton term in the 
action does not spoil the arguments. In the particular case of the geometries considered in this 
paper, $ diverges in the deep interior, hence this check especially becomes important. We check in 



appendix D.l that this term indeed remains finite in our case. 
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T, 



T 



Lattice gauge theory 


Gravity 




Spin model 


Deconfined, fJ(l)c 


Black-hole, $1) B 




Superfluid , f[l)s 


Confined, l{l)c 


Thermal gas, IXX)b 




Normal phase, L(l)s 



T 






One can summarize the various phases of the theories by the table above. The 
various U(l) factors in this table are as follows: The U(1)b is the dual symmetry 
that arises from compactifying the B field on the temporal circle. The U(l)c is the 
center symmetry of the corresponding lattice gauge theory that is proposed to arise 
in the large N limit of SU(N) (with or without) adjoint matter. Finally the U(l)s 
is the spin symmetry of the corresponding XY model. The arrow of increasing T is 
the same for the LGT and gravity picture and opposite in the spin model picture. 



3. A model based on Einstein-scalar gravity 

The arguments put forward in favor of a gravity-spin model correspondence above 
are general. In this section we would like to introduce a simple set-up which allows 
for computations of quantities such as the scaling of magnetization and spin-spin 
correlation function on the gravity side. The model is inspired by non-critical string 
theory and it becomes precisely non-critical string theory in the interesting regime 
near the continuous phase transition. 

3.1 The model 

The action in the Einstein frame reads, 

A = isk / d " xV= ~ 9 ( R - ^T (9 * )2 + VW ~ A 6 "*** + • ■ ') + G "■ 

(3.1) 
where the kinetic terms of the dilaton 15 and the B-field H = dB are inspired by 
non-critical string theory in d + 1 dimensions. The ellipsis denote higher derivative 
corrections. The last term in ( [3.1|) , that we shall not need to specify here, is the 
Gibbons- Hawking term on the boundary. 

We allow for a non-trivial dilaton potential V(&) that should be specified by 
matching the thermodynamics of the dual field theory. In the case of non-critical 
string theory in d + 1 dimensions the potential is given by, 

V nc ($) = J 2 ^\ (3.2) 



15 



The scalar field $ here is related to the original dilaton of the non-critical string $ by some 



rescaling that is defined in section 3.3. By $ we will always mean the "rescaled dilaton" throughout 
the paper. 



-14- 



where £ s is the string length and 5c is the central deficit, see section [3]3] for more 
detail. Gn in ( |3.1|) is the Newton's constant in D = d + 1 dimensions. It is related 
to N of the dual field theory 16 by, 

M^N 2 , (3.3) 



16ttG n " p 

where M p is a "normalized" Planck scale, that is generally of the same order as the 
typical curvature of the background L The limit of large N corresponds to classical 
gravity as usual. One should be careful in attaining this classical limit: The correct 
way of achieving this is described in section |3.3| . On the gravity side the parameter 
N arises from the RR-sector, where it is the integration constant of a space-filling 
F(d+i) form, i^d+i) oc N . Then the large N limit is defined as sending this value to 
infinity and sending the boundary value of the dilaton $ to — oo such that iVexp($ ) 
remains constant and yields M p in (|3T3"|). We refer to section EOfor details. 



In what follows we shall only consider solutions with either constant or pure- 
gauge -B-field whose legs are taken to lie along r and xq directions: 

B^ u = B r0 , (3.4) 

In this case H = in (|3.1| ) and the B-field contributes to neither the equations of 
motion nor the on-shell value of the action. However, it contributes the F-string and 
D-string solutions as we study in section 5. 

There are only two types of backgrounds at finite T (with Euclidean time com- 
pactified), with Poincare symmetries in d — 1 spatial dimensions, and an additional 
U(l) symmetry in the Euclidean time direction. These are the thermal graviton gas, 

ds 2 = e 2Ao{r) (dr 2 + dx 2 d _ x + dx 2 ) , $ = $ (r), (3.5) 

and the black-hole, 

ds 2 = e 2Mr) {f-\r)dr 2 + dx\_ x + dx 2 J(r)) , $ = $(r). (3.6) 

We define the coordinate system such that the boundary is located at r = 0. For 
the potentials V that we consider in this paper, there is a curvature singularity in 
the deep interior, at r = r s . In ( |3.5|) , r runs up to singularity r s . In (|3.6j ) there 
is a horizon that cloaks this singularity at r^ < r s where fifh) = 0. Xo is the 
Euclidean time that is identified as Xq ~ Xq + 1/T. This defines the temperature T 
of the associated thermodynamics. In the black-hole solution, the relation between 
the temperature and r^ is obtained in the standard way, by demanding absence of a 
conical singularity at the horizon: 

4nT = -f\r h ). (3.7) 



16 As explained above, N may either be the number of colors in SU(N) gauge theory or the 
number of spin states at each site in a Z^ spin-model. 
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This identifies T and the surface gravity in the BH solution. 

In the r-frame defined by (|3.5|) and (|3.6| ) one derives the following Einstein and 
scalar equations of motion from ( |3 . 1| ) : 

A" - A' 2 + ^L$'2 = o, (3.8) 

o — l 

f" + (d-l)A'f = 0, (3.9) 

(d - 1)A' 2 + A'f + A"f - -^-e 2A = 0. (3.10) 

One easily solves (|3.9| ) to obtain the "blackness function" f(r) in terms of the scale 
factor as, 



f 
1 _. JSL 



-(d-l)A 



f^ = l - J fr he _ {d _ 1)A . (3.H) 

Then the temperature of the BH is given by eq. ( |3.7| ) : 

T- 1 = 4n e «-Wk) / e ^ d ^ A ^dr. (3.12) 

Jo 

The difference between the entropy densities of the BH and the TG solutions is given 
by the BH entropy density up to 1/N 2 corrections 17 that we ignore from now on[15]: 



The difference in the free energy densities can be evaluated by integrating the first 
law of thermodynamics, [p6| : 



where r c is the value of the horizon size that corresponds to the phase transition 
temperature T(r c ) = T c , at which the difference in free energies should vanish. 

3.2 Scaling of the free energy 

In [|15|] we showed that there exists a continuous type Hawking-Page transition be- 



tween the TG and the BH solutions when the black-hole horizon marginally traps a 
curvature singularity: r^ = r c — y oo. This happens only when the IR asymptotics of 
the dilaton potential is chosen such that, 

V($) -> V^ e^r* (1 + V sub {$)) , $ -> oo (3.15) 

where V^ is a constant and V su b denote subleading corrections that vanish as $ — y oo. 
It is also shown in |15|] that the transition temperature T c that follows from ( |3.12| ) 



with r h — y oo stays finite. 

17 We choose to normalize the thermodynamic quantities by an extra factor of l/N 2 so that the 
entropy on the BH becomes 0(1) and on the TG it becomes 0(1/N 2 ). 
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Given the asymptotics in (|3.15|) one solves the equations of motion Q3.8Q and 
~M to obtain the IR behavior, as r — y oo, 

$(r) -> ^-^r + • • • (3.16) 

A(r) -> -|^r + ■ ■ • (3.17) 

where the sub leading terms vanish in the limit. 

Depending on V su b there are various different possibilities for types of transitions. 
We consider only two classes of potentials with: 

(3.18) 
(3.19) 



Case i : V sub = C e K$ , k > 0, 


$ -)■ oo 


Case ii : V sub = C $~ a , a > 0, 


$ ->■ oo 


Defining the normalized temperature, 




T-T c 
t- 





(3.20) 

the scaling of thermodynamic functions with t can be found from the following set 
of formulae: The reduced temperature directly follows from the subleading term in 
the potential, 

t = V mb ($ h ), (3.21) 

where $^ is the value of the dilaton at the horizon. Then the free-energy as a function 
of t follows from by ( [3.140 as, 

AF(t)(x f di e {d ~ l)A{i) . (3.22) 

Jo 

Here, the dependence of the scale factor on t should be found by inverting ( |3.21| ), and 
comparing the (leading term) asymptotics of the scale factor A(r) with the dilaton 
$(r) [|TJ. In the cases ( gig ) and fl37Lj) one finds that, 



Case i : A(t) = — r r hg(t/C) + ■■■ , t -y + (3.23) 

K{d — 1) 

Caseii: A(t) = 2 (t/C)~° + ■ ■ • , t -)• 0+ (3.24) 

The free energy then follows from ( |3.22| ) as : 

Case i : AF{t) oc t^ +1 , t -> + (3.25) 

Case ii : AF{t) oc e c '^*t 1+ -, t -> 0+, (3.26) 

where C = 2C* in the second equation. We see that F vanishes, as it should, for 
arbitrary but positive constants £, k and a. Other thermodynamic quantities such 
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15 



as the entropy, specific heat, speed of sound etc, all follow from the free energy above 

In the special case of 

(3.27) 



n — 1 
in ( p. 25 ) one finds an nth order phase transition. On the other hand, the special case 



of a = 2 in ( |3.19|) corresponds to the BKT type scaling 18 . 

One can also obtain the value of the transition temperature T c in terms of the 



coefficient of the dilaton potential in the IR as [15 



T c = ^~. (3.28) 

Finally, we should note the following issue. As mentioned above, the transition 
region t fa generically coincides with the singular region $^ 3> 1 in this setting. We 
do not need to worry about the a' corrections because they vanish in the interesting 
region r 3> 1 in the interesting limit r^l [15 . However, one should worry about 



the string loops. In a generic situation the higher string loops cannot be ignored 
near the transition region. We are however interested in the situation with g s — > 
(N — > <x>) that corresponds to the U(l) invariant spin-model. This can be achieved 
by sending the boundary value of the dilaton to — oo. We will now dwell on this 
point in more detail. 

3.3 The large N limit and string perturbation theory 

The effective Einstein frame action in fl3.1|) is supposed to arise from a (fermionic) 
non-critical string theory which also involves an RR-sector. The string frame action 
is, 

* = ^ J i* V^"" (*. + ^W + % - ^h) -sjsTIJf *?«>+■ ' ' 

(3.29) 
The ellipsis denote higher derivative (a/) corrections, subscript s denote string-frame 
objects and 5c is the central deficit that — depending on the fermionic or the bosonic 
string theory — reads 19 , 

5c = Cf (9 — d), fermionic; 5c = Cb (25 — d), bosonic. (3.30) 



18 Very recently holographic realizations of (quantum) BKT scalings were obtained in pTJ and 

19 The constants c/, Cb depend on the particular CFT on the world-sheet as there are various 
possibilities for the boundary conditions and GSO projections on the world-sheet fermions possible 
twisted or shifted boundary conditions for the scalar matter X^ [g9| . In the case of bosonic world- 



sheet with periodic scalars, one has Q, = 2/3 which is indeed what one obtains from solving (3.10) 



with the asymptotics (3.17) and (3.16). See the next section for details of the IR CFT. 
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F[d+i) is a space filling RR-form whose presence is motivated by holography: it should 
couple to the D^-i branes that are responsible for producing the SU(N) gauge group. 
As it is space-filling, its effect in the theory can be obtained by replacing it in the 
action by its on-shell solution [30]. This solution in general will be very complicated 
as the higher derivative corrections will also depend on F^d+i)- Let us ignore these 
higher derivative solutions for the moment in order to be definite — the following 
discussion will not qualitatively depend on the higher derivative corrections. 
The equation of motion for F^a+i) is d* F^+i) = 0. The solution is 

c F N e (d+1) 
F {d+1) = -5-7=, (3-31) 

where eu+\\ is the Levi-Civita symbol in d + 1 dimensions and cp is some 0(1) 
constant. We chose the integration constant to be proportional to N motivated by 
the fact that F should couple to N Dd-i branes before the decoupling limit. Inserting 
the solution in the action gives (we ignore the NS-NS two-form in the following 
discussion), 

A s = -^ J cf* V^"* (r. + Am 2 + |) + ||^ 2 + • ■ ■ (3-32) 

Now we define a shifted dilaton field 

$ = $ + logiV, (3.33) 

and go to the Einstein frame by 



We obtain, 



4 i 

SVi/ = e~ g^. (3.34) 



a= ^fJ dd+lx ^~ 9 { R - dh {d ® )2 + v{ ^ ) + ■ • • (3 - 35) 



where the dilaton potential becomes, 

V(*) = i (Sc e^ + le^ 1 * + ■ ■ ■) (3.36) 

We denote the corrections coming from the higher- derivative terms by the ellipsis. 
This is what one would obtain by ignoring the higher derivative terms 20 . 



20 In the phenomenological approach that we adopted in the previous section, one assumes that 
there exist a string theory that would produce a potential of the form ( 3.1 5| ) instead of ( 3.36 ). In 



particular the leading term with exponent 2(d+ l)/(d— 1) should either be absent or renormalized 
to4/(d-l). 
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On the other hand the solution of the dilaton equation of motion follows from 
( |3.29|) generically involves an integration constant that we shall denote as $o- For 



example in the kink solutions of [15| this corresponds to the boundary value of the 



dilaton on the AdS boundary. One can write 

$ = $ + «J$( r ) (3.37) 

to make explicit the integration constant. Now, we are ready to define the large-N 
limit. We send iV — > oo, $ — > — oo such that the shifted dilaton in ( [3.33J ) remains 
constant 

e ¥ «JV -> A =* e* = Ae 5 * (3.38) 

where A is some 0(1) constant. 

The shifted dilaton $ is the one that we used in the previous section to discuss 
thermodynamics and it is what we will refer in the next sections to study the ob- 
servables of the spin-system from the gravity point of view. Whether $ is large or 
small does not matter neither for the loop-counting of strings nor for the strength of 
gravitational interactions: The latter is determined by the coefficient in the action 
( p.35|) . Identification with fl3.1|) yields the Newton's constant 



2gd-l 

which shows that the gravitational interactions among the bulk fields can be safely 
ignored in the large N limit. This equation also defines the "rescaled" Planck energy 
that was introduced in (|3.3|) in terms of g s and £ s as, 



M p = £; 1 g s d - 1 . (3.40) 

The string loops on the other hand are counted by the coupling of the original dilaton 
$ to a world-sheet M with genus g as, 

e -h Im v^? (2) * = e -*ox(M) e -^ f M VhRWs* = NX ( M ) e -±f M vT^ 2 ^ ( 3 _ 41 ) 

where x{M) = 2(1 — g) is the Euler characteristic. We observe that the above 
definition of the large N limit does the job and suppresses the strings with higher 
genus. 

One might still worry about the viability of the string perturbation expansion if 
the additional term proportional to J M yhR^ 2 ^ in ( |3.41|) becomes very large in some 
limit. Indeed, as we argued above the interesting physics concerns the region $ 3> 1 
which corresponds to the vicinity of the phase transition. We check in appendix p.l 



and |D.2| that for all of the string paths that we consider in this paper the world-sheet 
Ricci scalar suppresses the linear divergence in $. For example in case of ( |3.18|) one 
finds that in the transition region yhR^ ~ exp(— k&). 
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All of the discussion we presented above can be understood in the following 
equivalent way. To be definite let us consider the simplest effective (rescaled) dilaton 



potential that corresponds to case 3.18 



V($) = Ko e^ # (1 + Ce~ K *) . (3.42) 



It was shown in [15) that this potential has a kink solution that flows from the AdS 
extremum at 

e*° = C^a (3.43) 

— where a is some number independent of C — to the linear dilaton geometry in the 
IR $ — > oo. Then the sub leading term in the potential can be written as, 

V sub {$) = a - K e K *°- K$ = a- K e K¥o e- K * (3.44) 

where we used ( |3.33| ). The statement that "the transition region corresponds to large 
dilaton" now can be quantified. What we really mean by this is that the reduced 
temperature t ( |3.2U| ) is small enough, so that the scaling behavior of observables set 
in. Now, from ( p. 21 ) we see that this is given as, 



t = a - K e K *°e-^. (3.45) 

On the other hand the large N limit ( |3.38| ) involves $o ~~ > — oo, therefore we see 
that in order for t to be small, one needs not the actual dilaton $ but the difference 
5$ = $ — $o to be large. The same reasoning can be generalized to general potentials 
that involve an AdS extremum. 

To conclude, we can safely ignore higher string loops in the computations below. 

3.4 Parameters of the model 

In the model constructed above there are various parameters. Here we shall list the 
parameters without derivation and refer to [Q for a detailed discussion. 

• Parameters of the action: In the weak gravity limit, Gn — > 0, N — > oo and 
Mp~ d = 16ttGnN 2 = fixed., there are two parameters in the action: M p and 
the overall size of the potential I. The latter fixes the units in the theory. 
One can construct a single dimensionless parameter from the two: M p £ which 
determines the overall size of thermodynamics functions in the dual field theory 
and it can be fixed e.g. by comparison with the value of the free energy at high 
temperatures, see | 26fl . In the present paper we are only interested in scaling 



of functions near T c , thus this parameter will play no role in what follows. 

Parameters of the potential: We have not specified the potential apart from its 
IR asymptotics. The IR piece will be enough to determine the scaling behaviors 
and also the transition temperature through equation (|3.28|). Therefore we 
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have only three (dimensionless) parameters: V^ 2 , C and /tora that appear 
in (|3.15|) , ( p,18|) and (|3.19|) . The first determines the (dimensionless) transition 
temperature T c i through ( |3.28| ), the second one is related to the boundary value 



of the dilaton (cf. the discussion above), and the third one determines the type 
of the transition. For example k = 2 for a second order transition, equation 

(EH- 



Integration constants: In |L5| we solve the Einstein-dilaton system and work 
out the thermodynamics in the reduced system of "scalar variables" that is a 
coupled system of two first order differential equations. One boundary condi- 
tion can be interpreted as the value of T, and the other is just regularity of the 
solution at the horizon. Therefore the only dimensionless parameter that arise 
among the integration constants is TL 21 



4. Non-critical string theory and the IR CFT 

4.1 Linear-dilaton in the deep interior 



The leading asymptotics ( 3.15 ) of the dilaton potential which follows from the re- 
quirement of a continuous Hawking-Page transition is precisely the same as the 
potential that follows from d + 1 dimensional non-critical string theory. This is 



easily seen by transforming (|3.1|) with the potential ( |3.15| ) to string frame with 



9s,iiv = exp(2$/(<i — l))g^ u . Not only that but we also have the asymptotics (|3.16|) , 
which imply that the asymptotic solution in the IR corresponds to a linear-dilaton 
background that is — very conveniently — an a' exact solution to ( |3.29| ) and corre- 



sponds to an exact world-sheet CFT. Indeed, in [15 it is shown that, with the 



subleading terms of the form ( |3.15| ), the string- frame curvature invariants both on 



the TG and the BH backgrounds vanish in the deep interior region near criticality 
i.e. for rh — > oo, (T — > T c ). Hence the higher derivative terms denoted by ellipsis in 
( p. 29 ) become unimportant in the IR theory. 



This implies that the dynamics in the transition region should be governed by the 
linear-dilaton CFT. More precisely, we expect that quantities that receive dominant 
contributions from the deep interior region near criticality should be determined by 
the linear-dilaton CFT. 
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In the fifth order system of (3.8-3. 9|) it is a little harder to work out the non-trivial integration 



constants. There it works as follows |26| : In (3J3), one requires / — > 1 as r — > 0. This fix one 



constant, and the other is gives T. In the rest, one is fixed by requirement of regularity at the 
horizon, one is just a reparametrization in r, and the last is fixed either by the asymptotic value of 
the dilaton in the case $(r) = <t>o is constant at the boundary, or the integration constant A that 
determines the running of the dilaton near the boundary $ ~ log(— log(Ar)) near r — > 0. In either 
case the thermodynamic functions can be shown to be independent of this constant pq] . 
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In the next section we shall make use of this observation to argue that the various 
observables in the corresponding spin-model scale precisely with the expected critical 
exponents near T c . 

Another implication of this is that an asymptotically linear dilaton geometry 
(with corrections governed by the subleading terms in ( |3.15| ) develops an instability at 
a finite temperature T c into formation of black-holes. It is quite reasonable to expect 
that in the limit of weak g s this point coincides with the Hagedorn temperature of 



strings on the linear-dilaton background flTT]] . We have more to say on this in section 



Finally, we note that in the case when the model is embedded in non-critical 
string theory, all of the parameters in the model are entirely fixed. To illustrate this 
let us assume that the entire potential is given by the leading term, ignoring the sub- 
leading terms etc. Then the coefficient V^, in (|3.15| ) and the transition temperature 
would be given as, 

W = Cbi2 \~ d \ T Ctnc = -L^ Cb (25-d), (4.1) 

in the case of bosonic world-sheet CFT and 

Vocnc = ^A T c , nc = J^ C/ ( 9 -d), (4.2) 

S S 

in the case of fermionic world-sheet CFT. These results follow from ( |3.3U| ) and ( |3.28| ). 
Of course, in reality these numbers should be renormalized because the theory is not 
just given by the leading piece: a potential with only the leading exponential behavior 
do not possess any phase transition. The corrections will depend on the UV physics 
where the a' corrections kick in and renormalize these coefficients. We shall argue 
for another way to fix these numbers in section [5| We will also show in that section 
that the scaling exponents are also determined completely, once the CFT is fixed. 

4.2 The CFT in the IR 

The arguments presented above point towards the conclusion that, on the string side 
the criticality of the dual spin-system should be governed by a linear-dilaton CFT. 
Here we want to spell out some of the salient features of this IR CFT. We start with 
the bosonic case and then mention generalization to fermionic CFT in the end. 

We reviewed the intimate connection between non-critical string theory and the 
linear-dilaton background in appendix [FJ. Utilizing this relation one can obtain the 
stress-tensor of the (bosonic) linear-dilaton CFT as 



T{z) = --■. dX»dX„ : +v^d 2 X^ (4.3) 

a 

for the left-movers, with an analogous expression for the right movers, v^ are the 
proportionality constants in the dilaton solution $ = v^X^. The indices are raised 
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and lowered by the flat metric. The total central charge of the theory (including the 
ghost sector) vanishes for v^ satisfying ( [B.5|) . In our case we have, 

«Ai = — ^ <W ~ m <W- ( 4 - 4 ) 

The reason for denoting this constant itlq will be clear when we analyze the spectrum 
of fluctuations in this geometry, see appendix [E|. The total central charge of the 
theory (including ghosts) vanishes only for, 



25-rf 
6fl 



bosonic 



< = 8 _ ( ■ • (4-5) 

1 ^# fermiomc 



for the bosonic and fermionic CFT's, [31 



Now we discuss the spectrum in the case we are interested in: The Euclidean d+1 
dimensional world sheet with (|4.4j ) and the Euclidean X° dimension compactified on 
a radius R = 1/2-kT. There are various ways to obtain the spectrum. Both the 
light-cone and the covariant quantization is discussed in |31|, ^9|. Here we trivially 
extend these results in our case. 

The Virasoro generators are now 

1 .°°. r~' 

L m = ^ Yl '■ a ™~n a n,» ■ +i ^f( m + 1)^0^- (4.6) 

n=— oo v 

The center-of-mass momenta are related to the zero mode oscillators as usual, p^ L = 
^ojq and ]f R = y Jr«o- Decomposing into components one has, 

w w 

p L = 2irTk + , PoR = 27rTk , (4.7) 

F0 ' 2rrTa r F ' 27cTa r V ; 

Pi,L = Pi,R = Pi, Pr,L = Pr,R = Pr- (4.8) 

In the first line the integer k denotes the Matsubara frequency and the integer w 
denotes the winding number on the time-circle. As a result of the linear piece in the 
zeroth level Virasoro generator (f4.6|) one obtains the following mass-shell conditions 



(we adopt the definition of mass in f29j) in the light-cone gauge 

/ w \ 2 2 



-m 2 d+1 =p{+p 2 r + 2tm oPr + (27rkTy+(^- f - ; ) = ~(N + N - 2), (4.9) 
= kw + N -N, (4.10) 

where p±, N and A^ denote the center-of-mass momentum, the left (right) number 
of oscillations in the space transverse to motion, respectively, 

oo oo 

N = ^2 a±- n ■ a_L, n , N = ^ ®±-n • &x,n- (4.11) 

n=l n=l 
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In (4.9) m d+l denote the d + 1 dimensional mass. One important difference 



between the linear-dilaton and the flat case is that the definition of the mass of the 



string excitations in terms of their momentum gets modified^] due to the linear 
oscillator piece in ([Op. The flat case follows by setting m = 0, hence sending 
dilaton to constant. 

Once the modified definition of mass is attained, the physical spectrum of the 
linear dilaton is exactly the same as the critical string: the lowest level N = N = 
is a tachyon with mass —4/a', the next level is massless and corresponds to the 
fluctuations of the metric, the B-field and the dilaton, etc. p9| . 

All of these results are readily extended to the fermionic case with N = 1 world- 
sheet supersymmetry [ 29| . In the light-cone gauge, one obtains the following spectra 



for the NS and Ramond sectors, 



,„:,_< =N + J2<ib-X-l> ( NS ) ( 4 - 12 ) 



g>0 



m d+1 



N + Y.^-'Pt ( R i ( 4 - 13 ) 

<3>0 



where N denotes the number of bosonic oscillations in the transverse space ( |4.11| ) 
and geZ for the R-fermions and q G Z + | for the NS-fermions. This is again the 
same spectra that one finds in the critical super-string. 

However one finds crucial differences at the one-loop level: Modular invariance 
does not allow for NS-R fermions except in particular dimensions given by multiples 
of 8. This is quite convenient for our holographic purposes, because we do not want 
any fermionic operators in the dual spin- model. Thus in a generic dimension d+1 < 8 
one has only two sectors R-R and NS-NS. Furthermore, in the generic case, there is 
no analog of the GSO projection of the superstring. Therefore the tachyon in the 
NS-NS sector survives. 

Existence of tachyon in the physical spectrum is a very generic feature of the 
linear-dilaton CFT in any dimensions. The mass of the tachyon changes depending 
on which particular CFT chosen. With the definition of mass adopted above it is 
given by m\ = —4/a' for the bosonic case, m\ = —2/a' for the NS-NS fermions, 
m 2 = — 15 /4a' for an orbifold in the r-direction, etc., but we stress that the ground 
state for k = w = in linear-dilaton CFT in arbitrary dimensions is always a 
tachyon. 22 

This fact renders the linear-dilaton theory unattractive from many perspectives. 
In our case however, it is a desired feature of the IR CFT. We recall that the back- 
ground geometry becomes asymptotically linear-dilaton only in the transition region 



22 With a more conventional definition of mass PJJ , one finds a tachyon only for d > 1 in a d + 1 
dimensional theory. In our case, the equivalent statement is that if we consider propagation of the 
tachyonic mode, we find a smooth propagation for d < 1 but oscillatory behavior for d > 1. 
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T ~ T c and only in the large r region. We do not expect that the complete sigma- 
model which corresponds to the black-hole for an arbitrary T have tachyon as a 
ground state. This would imply that the black-hole geometry is unstable at any 
temperature. Instead the linear dilaton CFT describes the physics near the transi- 
tion and we do expect instability in this region. In fact, as we show in sections ( |5.3.2| ) 



and (5.4.2), it is the presence of the tachyon which guarantees vanishing of magneti- 



zation as M ~ (T — T C )P and divergence of the correlation length as £ ~ \T — T c 
at the transition! 



5. Spin-model observables from strings 

F-strings and D-branes are important probes in the standard examples of the gauge- 
gravity correspondence. In case of the holographic models for QCD-like theories, the 
phase of the field theory at finite temperature, the quark-anti-quark potential, the 
force between the magnetic quarks etc, can all be read off from classical F-string and 
D-string solutions in the dual gravitational background. In this section we argue 
that the probe strings constitute indispensable tools also in the spin-model-gravity 
correspondence. In particular, the Landau potential, the correlation length, the 
various critical exponents, the scaling of order parameters near the transition, the 
phase of the system, spin-spin correlation function, etc. can all be computed from 
the probe strings in the dual background. In this section we discuss how to obtain 
the various observables of the spin-model from the probe string solutions. 

5.1 What can we learn from the Gravity-Spin model duality? 

In order to answer this question, one has to identify the Landau and the mean-field 
approximations on the gravity side. The Landau approach is based on integrating 
out the "fast" degrees of freedom in the spin-model in order to obtain a free-energy 
functional for the "slow" degrees of freedom, i.e. the order parameter M. We refer to 
appendix |C] for a review of the statistical mechanics background and in particular a 
description of the Landau approach. This is exactly analogous to integrating out the 
gauge invariant states to obtain an effective action for the Polyakov loop on the LGT 
side, as illustrated in appendix 0. In the context of the gauge-gravity correspondence, 
this is, in essence, very similar to keeping only the lowest lying degrees of freedom 
in string theory, i.e. the supergravity multiplet. It is tempting to think that the 
complicated step of integrating over the spin configurations in ( |C.4j ) to obtain ( |C.5| ) 
can be side-stepped by use of the gravity-spin model duality 23 . 



23 We note a very interesting paper |32| that dwells on these issues. In this paper Headrick argues 
that one can generate the Landau functional at strong coupling in terms of classical string solutions. 
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In the most general case, the correspondence between the spin model and gravity 
should relate the Landau functional flC.5|) with the string path integral 24 : 



J st- 



(5.1] 



As in the original gauge-gravity duality we expect that there is a simple corner of 
the correspondence where both sides of ( p.l[ ) become classical and one approximates 
the path integrals by the classical saddles. 

This is the large N limit: On the LHS this is given by the Landau approximation 
( P.6|) . On the RHS, it is given by the saddle-point approximation to the string theory 



where one can ignore string interactions g s — > 0. Then ( |5.1|) reduces to, 

e -PF L = e -Ast^ ( 5 _ 2 ) 

where the action on the RHS is the full target-space action including all the a' cor- 
rections, evaluated on the classical saddle. It is the effective action for all excitations 
of a single string 25 and in principle it can be obtained from the sigma model on the 
world-sheet. 

At this point, it is clear that scaling of any quantity on both side of fl5.2|) near 
T c should be characterized by the mean-field scaling. This is just a consequence of 
the saddle point approximation. Therefore, any operator in the spin-model that is 
given by a fluctuation of A st should obey the standard mean-field scaling. We shall 
refer to these operators as local operators. The only possible exceptions to this — 
within the classical approximation of ( |5.2| ) — are operators that can not be obtained 
as fluctuations of A s t- These correspond to non-local operators on the gauge theory, 
they are governed by probe F-strings or D-branes on the string side. Yet, as we will 
show in the next section, they can correspond to quite ordinary quantities such as the 
magnetization on the spin-model side. Thus, magnetization is an example of a non- 
local operator. Even for the "non-local observables" though the mean-field scaling is 
expected to hold in a semi-classical approximation, where one only keeps the lowest- 
lying string excitations in string path integrals. These excitations correspond to bulk 
gravity modes (levels N = and N = 1 of the string spectrum). We confirm this 
expectation in the sections ( |5.3.2|) and (|5.4.2j ) below. 



In practice, it is usually very hard to reckon with ( |5.2| ), and one further considers 
the weak-curvature limit where one can replace the RHS with the (super)gravity 
action: 

-F L « TV^xCgr. (5.3) 

Here, C gr is the (super)gravity action evaluated on-shell, on the classical saddle. We 
also assumed a trivial dependence on the spatial volume and made use of the fact 



24 We shall be schematic in what follows. 

25 It is important to note that this is not a string field theory action, the excitations governed by 
A s t are particles, rather than strings. 
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that the temperatures on the spin-model and the gravity sides are inversely related, 
cf. appendix |A|. 

Influenced by the standard lore of the gauge-gravity correspondence, we expect 
that the weak curvature limit corresponds to strong correlations on the spin-model 
side. On the other hand, one quantifies "strong correlations" by the Ginzburg crite- 
rion in the spin-model, as reviewed in section |C| Quite generally, the system will be 
in a regime of strong correlations around the phase transition where the mean-field 
approximation usually breaks down. Therefore, one may hope that the gravity side 
provides a better description in Q5.3| ) precisely within this interesting region. This 
can be checked explicitly by computing curvature invariants in the string frame. 
Even though one show that the Ricci scalar (and the various contractions of Ricci 
two-form and the Riemann tensor) vanish in the limit (see fll5| ) there exists invari- 



ants such as d$ 2 that asymptote to a constant that is generically the same order as 
the string length scale £~ 2 . Therefore, in a generic case one is forced to include the 
higher derivative corrections. Luckily this can be done precisely in the interesting 
critical regime, because the background asymptotes to a linear-dilaton theory. 

What observables can we actually calculate on the gravity side? Because going 
beyond the large N limit is very hard, one can (at present) only hope to obtain 
results in the Landau approximation. The main observables then include the Landau 
coefficients 26 chq(T), Qi(T), ot2(T), the basic scaling exponents (3, u, r], 7 etc., and 
the spin correlation functions. Moreover, the scaling exponents of operators that are 
dual to fluctuations of the bulk fields in C gr in ( [5.3D are necessarily given by the 
mean-field scaling. Therefore one can only hope to obtain results beyond mean-field 
in the scaling exponents of operators dual to stringy objects, such as magnetization 
or the spin-spin correlator. 

Once again, we would like to emphasize the distinction between "mean-field 
scaling" and the "mean-field approximation". The former is unavoidable for local 
operators in the Landau approximation (large N). On the other hand, gravity de- 
scription is expected to go beyond the latter. Therefore for quantities such as T c , the 
Landau coefficients at T c , etc., and correlation functions of the non-local observables 
we expect gravity to provide better answers than the mean-field approximation. 

One may still ask the question, what is the use gravity-spin-model duality if one 
can compute all of these quantities by employing Monte-Carlo simulations, or RG 
techniques? First of all, the RG techniques are limited in the case of strong corre- 
lations. Secondly, the calculations on the gravity side are much easier to perform, 
much easier than the Monte-Carlo simulations, and one can usually obtain analytic 
results. However a more fundamental reason is that, there are situations where ap- 
plicability of the Monte-Carlo simulations are limited. The well-known examples are 
the computation of real-time correlators or spin-models with fermionic degrees of 



26 We refer to appendix for a definition of these coefficients. 
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freedom. By the gravity-spin-model correspondence, one expects to overcome such 
fundamental difficulties. 

5.2 Identification of observables 

The duality between the lattice gauge theories and spin-models [[10], [[Up relate the 
magnetization directly to the Polyakov loop. On the other hand, the Polyakov loop 
is related to the classical F-string solution as discussed in section |2.3| . Therefore we 
propose the following chain of relations: 

(P{x)) <-> (m(x)) <R> e ~ SNa[Cx] . (5.4) 

Here the boundary condition C x for the string is just a point x in the spatial part 
and a loop on the temporal circle. 

The spin field is valued under U(l)s- Similarly the Polyakov loop is valued under 
the center C = Zn that becomes a U(l) in the large N limit. One should think of 
this as the exponents becoming angles in the transformation, 

P^e 2lH ~*P, k = l,2---N 

at large N. We shall denote this U{1) as U{l)c- Similarly, as discussed in section 
O at length, the F-string that winds the time-circle is charged under the U{1)b 21 , 



because it couples to the B-field. Thus one should identify 

U(l) s = U(l) c = U(l) B . (5.5) 

as in table in section 12TB. 



One should work out the identification in ( |5.4D carefully. In particular the first 
entry is a complex number and the second entry is a vector in 2D spin space. The 
precise identification of the two is provided with the standard isomorphism between 
U{1) and 0(2) representations. We imagine the vector fh on the XY plane repre- 
sented by the magnitude \m\ and the phase ip. Then the simplest option is to set 
m x = Re(P) and m y = Im(P). There is a little complication though, because in fact 
the identification should depend on the value of ip. This is because the physically pre- 
ferred reference frame is set by the direction of the magnetization vector Vi in ( |C16| ). 
All of the correlation functions should be decomposed into components parallel and 
perpendicular to Vi. Represented by the phase, the direction of magnetization reads 

v = (cos(^), sin(-?/>)). (5.6) 

Thus, the naive identification mentioned above is correct only for ip = 0. For a 
different value of ip one should obtain the correct identification by a U(l) rotation: 
P = exp(iif))(m x + irriy). Thus, in general we have, 

Re(P) = m\\ = cos(ip)m x —sm(ip)m y , Im(P) = m± = sin(^)m a; +cos('0)m y , (5.7) 



27 The charge is determined by the winding number. Here we are only interested in strings that 
wind the time circle once. 



-29- 



where 

m\\ t i = v ■ m Vi, m± ! i = (Sy — ViVj)rrij. (5.8) 

The identification of the second and the third entries in ( p.4|) is straightforward. 
One can schematically write, 

(P[C]) = (e- /G+j/B >, (5.9) 



using Q2.71) and ( |2.9|) , where we dropped the dilaton coupling 28 . Thus the magnitude 
of P is determined by the space-time metric and the phase is determined by the 
B-field as explained in section |2.3| and one should identify the if) angle defined in 
( |2.4j) and the angle defined in ( |5.6| ). 

In this picture the spin correlation function should be given by a fundamental 
string solution that ends on two separate points x and y: 

{mi{x)mj{y)) «-> e - Sira[Cm * ] . (5.10) 

The boundary condition is such that the string ends on the points x and y on the 
spatial parts and wraps the temporal circle. 

Again, one has to be careful in the identification ( |5.10| ) and has to split the 
correlator into the parts perpendicular and parallel to the magnetization vector v^ 
as in ( |Cl9| ): 



(rrii(x) rrij(0)) = (rh\\(x) ■ m\\(0))viVj + (m±(x) ■ m±(0))(5ij — v(Vj). (5.11) 

On the other hand one has the identification, 

(m{x) • m(0)) = (P*{x)P{0)). (5.12) 

Given the identification ( p.7| ) one obtains, 

(m,l(x) • m||(0)) = \M\ 2 + (s\\(x) ■ s]|(0)) = (Re P(x)Re P(0)), (5.13) 
(m ± (x) ■ m ± (0)) = (s±{x) ■ s±(0)) = (Im P{x)Im P(0)) , (5.14) 

where we also decomposed the magnetization M and the fluctuations s according to 
( |C.14j) , assuming that M is isotropic. This is of course in the black- hole phase. In 
the thermal-gas phase \M\ vanishes and any direction Vi is identical. 

5.3 One-point function 

Having identified the observables on the spin model side with the observables on the 
gravity side, we are ready to determine the magnetization M of the spin model on 
the gravity side by a one-point function calculation. As we argued in the previous 



28 We check in appendix D.l that this contribution is sub- leading and do not contribute to the 
scaling near T c . 
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section the magnetization should be given by the real part of the F-string solution 
that wraps the time circle. As a warp-up exercise, we shall first assume that the 
string path integral is dominated by the classical saddle and obtain the resulting 
scaling law for the magnetization. After this, we will loosen the assumption and 
perform the same calculation in a semi-classical regime in section ( |5.3.2| ). 

5.3.1 Warm-up: classical computation 

The definition of the Polyakov action and the boundary conditions are given in 
detail in section [2.2| . In Appendix D we prove that in all of the cases we consider 
in this paper, the dilaton coupling term in Snq, that is given by &R^ gives finite 
contributions, hence do not alter the scaling. Also, the effect of the B-field is discussed 
in detail in section 3.2. Thus we shall only restrict our attention to the area term 
(see eq. (|5.9|) ): 

\M\ = | (P[C]> | oc (e-f G ) (5.15) 

and replace the fundamental string action with the Nambu-Goto action. 

To compute the energy of the string, we fix the gauge a = x°, r = r where (r, a) 
are the coordinates on the world-sheet, x° is the Euclidean time that is identified as 
x° ~ x° + 1/T and r is the radial variable in the coordinate system given by (|3.6|) . 
Then, 

Sng = ' j dr^deth ab , h ab = <9 a x M c\x^, (5.16) 

where £ s is the string length, g s is the BH metric in the string frame: 

ds 2 = e 2A s (r) (f-\ r )dr 2 + dxl x + dt 2 f(r)) , AJr) = Air) + — ^— $(r), (5.17) 
v ' d — 1 

and e is a point near the boundary 29 . 

In passing we review the discussion in section |3.3| . There we introduced the 
rescaled field $ in relation to the "real" dilaton as in equation ( |3.33| ) as $ = $ + 
log JV = $ + 5$ + log iV = log A + 5$ the last two lines follow from ( |3T37|) and (|3T38| ) . 
The constant A is 0(1). Thus, when large, $ just corresponds to the difference of 
the real dilaton and its boundary value $o- Large $ does not mean large $ when 
$o is chosen very small. This choice indeed corresponds to the large N limit in the 
gravity language. Therefore we can safely ignore the loop corrections here, and in 
the next sections. 

On the TG solution one replaces / — y 1, A(r) — > A (r) and $(r) — > $o( r ) m 
the above formulae. As described in section appendix |D.1| the exponential of —Sng 



29 The target space metric typically diverges on the boundary and this cut-off guarantees finiteness 
of ( 5.16| ). One can remove the dependence on e by some renormalization procedure, however we do 



not need this as we are only interested in the dependence of Sng on r^ in the limit r^ —> oo that 
corresponds to T — > T c . We provide an appropriate renormalization scheme in appendix p.l|. 
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vanishes on the TG solution. Thus, 



M 



TG 



0. 



(5.1? 



This confirms our discussion in section [2.3| , that indeed t/ie TG phase of the gravity 
corresponds to the disordered phase of the spin model. 

Let's turn to the BH phase. It is shown in ( D.l ) that (|5.16|) is finite, unless 
d = 2. In the latter case the fluctuations of the zero mode of the B-field in the 2D 



transverse space makes it vanish — see the discussion in section [2.3[ Thus one finds, 

Mm^O (d>2), Mbh = (d = 2), (5.19) 

and the black-hole solution indeed corresponds to the ordered phase of the spin-model 
for d > 2. One also confirms on the gravity side, that no long-range order in the 
spin field is possible in d = 2. 

We would also like to see how M scales near the phase transition (as one ap- 
proaches from below in the spin model). Using eqs. ( |3.17| ) and ( 3.16 ), one finds that 
the scale factor vanishes A s (r) — y as r — y oo, see ( p.8|) . Thus the integrand in 
( [5.1 6| ) becomes constant in the limit r^ — y oo. Using also the fact that T — y T c in 
this limit one finds, (see App. p.l| for details), 



T -i 

SnG ~* 2^f h ' 



T ^T r . 



(5.20) 



In order to determine the scaling of Sng with the reduced temperature t, one needs 
to find the dependence of r/, on t . This is done in App. |D.1| . Define the dimensionless 
constant: 

V. = V^l = (4ttT c 4) 2 , (5.21) 

where V^ is defined in (|3.15|) and (|3.28|) is used to relate it to T c . Then the result is, 



e -S NG K tl -v- 3 
j—Sng 



i->0 



Case i : M BH 
Case ii : M BH = e~ 0jvG oc e'^^c) ^ t-> 0. 

where the constants a, n and C are defined in ( |3.18| ) and (|3.19|) . 



(5.22) 
(5.23) 



We note that ( 5.23 ) is valid strictly for d > 2. As mentioned before, for d = 2 we 
obtain M = below and above the transition. 

Let us now specify to the case of second-order transitions. Then the coefficient 
k is given by ( |3.27| ) with n = 2: 



second — order transition. 



(5.24) 



Then, comparison of ( |5.22| ) with (|C.21| ) yields the critical exponent of the magneti- 
zation as, 



= 2V- 



(5.25) 
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where V s is given by (|5.21| ). 

Finally, we note that the mean- field scaling \M\ ~ t? corresponds to a particular 
value of the parameter V s : 

V S MF = 4. (5.26) 

It seems like a contradiction that one does not automatically obtain the mean-field 
scaling for \M\ directly from the gravity action. However, it is not. As explained in 
section |5.1| , the magnetization is a "non-local operator" which maps onto a non-local 
object in the string theory side, i.e. the expectation value of the F-string that wraps 
the time-circle. Therefore the classical string computation is not bound to produce 
the mean- field result. 

On the other hand, we remind the reader that this section is just meant to be 
a warm-up exercise. The classical computation is not at all guaranteed to be self- 
consistent. In particular it assumes that the string path integral is dominated by 
classical saddles, which only holds when £/£ s is parametrically large. This is not 
guaranteed in the backgrounds that we discuss in this paper. Below, we consider the 
semi-classical computation and argue that the classical result is altered non-trivially 
due to large quantum fluctuations. We shall observe that the mean-field scaling 
arises as a result of the semi-classical computation. 

5.3.2 Semi-classical computation 

In principle the classical saddle dominates the path integral only in a regime where 
the typical curvature radius of the geometry — that is determined by the asymptotic 
AdS radius £ — is much larger than the string length £/£ s ^> 1. This is indeed the 
case for a pure AdS black-hole geometry when the dual M = 4 theory is at strong 
coupling, because the AdS/CFT prescription relates the ratio £/£ s to the 't Hooft 
coupling A t of the dual gauge theory as £ 2 j£\ oc \f\f In the theories we are interested 
in, when there is no tunable moduli like At, this assumption will generically fail — 
unless there is some physical reason for £/£ s to be large. In this section we consider 
a full path integral computation. 

What kind of a string propagator do we want to compute? In the classical 
approximation of the previous section, the recipe |35|, [J7], |TB[] to compute (P[C]) can 
be described as follows. Consider a string that stretches between the boundary at 
r = and a probe D-brane just outside the horizon at r = r^ — e. The boundary 
operator wraps the time-circle, hence the string that couples to it also should. In 
the Euclidean BH the length of the time circle measured by an observer sitting at 
r goes to zero as r — > r h thus the string world-sheet wraps a 2D ball, and yields a 
finite answer when the UV divergence regularized properly. This string world-sheet 
is the classical saddle of the Nambu-Goto action and it provides the correct answer 
for £ ^> £ s . We can generalize this picture to the case £ ~ £ s simply by considering a 
string that stretches between the boundary and the horizon, wrapping the time circle, 
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but this time computing the full path integral including all quantum fluctuations on 
the string. This is similar to an open-string annulus diagram. It is very hard to 
compute this however because the string stretches over the entire range between 
the boundary and the horizon and one needs the full CFT that governs the physics 
everywhere on the target-space. Instead one can think of this diagram as propagation 
of a closed string that is created at the boundary, travels the distance from r = 
to r = rh on the BH and absorbed at the horizon. This is easier to handle because, 
at least we know the CFT close to r^ in the limit th — > oo, which corresponds to 
the phase transition regime. This is the linear-dilaton CFT described in section f|. 
Indeed, this CFT will prove important in determining the critical exponents of the 
corresponding spin system. 

Then the idea is to divide the closed string paths into two parts: from the 
boundary to a point r m and from r m to the horizon rn- The point r m should be 
chosen such that the string propagation from r m to r^ be governed by the IR CFT, 
see figure |3|. What is meant by "semi-classical approximation" will be to consider 
the contribution of the lowest mass string states at levels N = and N = 1 in the 
IR CFT. 

Field theory analogy: It is helpful to introduce the idea first in a similar situation 
in quantum-field theory. Generalization to the string will then be clear. First consider 
the correlator of a free massive scalar field with mass m 2 in flat d + 1 dimensions 
(0(0) 4>{y)) ■ This can easily be given a point-particle interpretation 30 , 

/■oo 

{4>(x)4>(y))= / dr(0,0\r,y) op (5.27) 



where the integrand is just the propagator of a point-particle in proper time r with 
Hamiltonian p 2 + m 2 . 

Now let us consider the more general situation of computing the propagator of 
a field in curved space time. The field should be specified with some quantum 
numbers such as momenta, charge etc. determined by the isometries of the back- 
ground. We assume that there are no self-interactions, hence no Feynman loops. We 
also assume that the back-reaction on gravity can be ignored. Finally we consider a 
background geometry of the "domain-wall type" |3.6| where one coordinate r is sin- 
gled out. We denote the d + 1 coordinates as (r,x). Now, ((f>(Q,0)(f>(r,y)) can still 
be formulated in proper-time as in ( p.27[ ) but this time the one-particle Hamiltonian 
will be much more complicated. 

However, let us consider a situation when the background geometry simplifies in 
some asymptotic region, when r ^> 1, where we know how to write down the one- 
particle Hamiltonian. Then the idea is to divide one-particle paths in ( |5.27 ) from 



:»,(] 



We consider Euclidean case for simplicity. 
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to r m and from r m to rh, where ru > r m 3> 1. For this purpose we decompose the 
correlator as, 



((j)(0,0)4>(r h ,y)) = dr drd d x m (0;0,0\T m ;r,x m )(T m ;r,x m \T;r h ,y m ) (5.28) 

~ I ™< I uT m J m (I X m yU, U, U|T m , T m , X m j \T m , T m , X m \T , T } ll ym/IRi 
JO J 

where in the second line we exchanged the integral over the intermediate point r 
with an integral over r m producing a Jacobian J m 31 . 

In the second line of ( p. 28 ) the propagator in the region < r < r m is governed 



by an unknown one-particle Hamiltonian Hjjy that is the full Hamiltonian valid 
everywhere on the target-space. The second propagator is governed by Hjr for which 
we assume the knowledge of the spectrum. The two should be continuously connected 
at r m . The approximation in the second line is to replace the full Hamiltonian in the 
second propagator with this IR Hamiltonian. The entire procedure will in general 
depend on the matching point r m . But, if we are interested in how the object scales 
as a function of the end-point r^, this dependence will be irrelevant. 

A technical but crucial point is that the division of the paths as in the second 
line of ( j5.27| ) makes sense only for the paths with r(r) > 0. This will certainly be 



satisfied for "slight deformations" from the classical path, if the classical path itself 
satisfies r(r) c i > 0. It is reasonable to assume that these are the paths that dominate 
because they minimize the kinetic energy in the one-particle lagrangian L op oc r 2 + 
• • ■ . The procedure can be extended to non-monotonic paths in an interesting way. 
Let us consider one-dimensional case for simplicity (the generalization to arbitrary 
dimensions is trivial). Suppose that we want to divide the path integrals in two 
different regions in space, for r < r m and r > r m . Then, one has to classify all of 
the paths according to their "crossing number" c m that is defined as the number of 
solutions to r(r) = r m . The monotonic paths have crossing number c m — 1. This 
is obviously an odd number and the next case have c m = 3, see figure 0. All of the 
paths from r = to r = r^ > r m are classified by c m . For non-trivial paths, with 
c m > 1 one can apply the same procedure by defining r m to be the greatest solution 
to r{r m ) = r m . Then, the procedure applies smoothly. For sake of the argument here, 
we will restrict only to the paths with c m = 1. This can be achieved by choosing 
r m to be close enough to r^. This is indeed the case in the physical situation we are 
interested in, because the region where the CFT on the string is governed by the IR 
CFT corresponds to r m < r h for r h >> 1. 

31 This is achieved by making use of the freedom to choose r m anywhere in between and r 
and inserting inside the integral 1 = J dr m { - r ^ m >~ rm > a j a Faddeev-Popov. This is not to be 

confused with the usual re-parametrization invariance of the relativistic point-particle. Here we 
describe propagation of a quantum field in the Schwinger's proper-time formulation, not a relativistic 
particle. In particular the Lagrangian that generates the propagation is not re-parametrization 
invariant. 
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Figure 2: Division of the one-particle paths that contribute to the QFT correlation func- 
tion into regions with r < r m and r > r m . The paths are classified according to their 
crossing number c m . For crossing c m > 1, r m should be chosen as the greatest node. 

At this point it is helpful to switch to canonical formulation and express the 
propagators in terms of the eigenstates of H uv and H IR that we denote as £' and £ 
respectively: 



(0(0 



/»oo p 

0)(p(r h ,y)) = dr dr m J m d d x m ^ ^(0,0)^*,( 

^0 J C'al-l T T,r 



x ^ ^(r m ,£J^(r,,y)e^^««\ 



-T m H UV (£>) 



(5.29) 



where "H denote the Hilbert spaces of the respective Hamiltonians. ^ and ^ denote 
the wave-functions of the eigenstates £ and £' of ///# and H^y respectively. 

One can now carry out the integration at the matching subspace r = r m over 
x m . This integration will produce the overlap of the wave function \& of the UV 
Hamiltonian with the wave function \t> of the IR Hamiltonian 32 : 

Cw = f d d x m ^l(x m )^^(x m ). (5.30) 

One can further sum over the UV Hilbert space, by defining the overlap function 

^(0,r m )= J2 C^(0)e- H ^>™. (5.31) 

t'eHuv 

This is the amplitude for production of a state £ of the IR Hamiltonian at r m . Thus 
one has 



(0(O ; 



rfr / rfr m J m ^ ^(0,r m )^(r h ,^e-( T - Tm ) H ««). (5.32) 



32 We do not assume that the Hilbert spaces of the UV and the IR Hamiltonians have same 
dimensionality, the overlap matrix may be rectangular. 
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Now one can carry out the r integral; this will produce Hm(£) in the denominator 
and then the sum over the IR states £ will only get contribution from the on-shell 
state £* with H UV (C*) = 0. This is the analog of the on-shell state p 2 + m? = in 
the free-field case. Therefore our final expression is, 

(0(0,0)0(7-^)) = (J dr m J m A^(0,T m )e TmHl ^A %(r h ,y) (5.33) 

Dependence on the matching point r m of the UV and IR regions is hidden in the 
Jacobian J m . Our ignorance about the UV region of the target space is summarized 
by the function A^(0,T m ). More generally this may be replaced by a sum over the 
on-shell states £* as will be the case for the string propagation below. 

Closed string case: Having outlined the procedure for the simpler case of quan- 
tum field theory, let us now consider the closed string propagation. As we argued in 
the beginning of this section, the one-point function (P[0]) is given by the propagator 
(\&3, 0, 0\^ff,rh,Xf) with some initial state \Pj on the boundary that corresponds to 
the Polyakov loop at the transverse point 0, and some final state \l// at the horizon 
at a transverse point Xf. In the end of the computation everything that is not de- 
termined by the boundary condition at r = should be summed over. In particular 
we should integrate over x/. When comes to \I// the situation is as follows: In the 
CFT language, the path integral we want to compute is a sphere diagram with two 
insertions of vertex operators V^.(o"!, 7"i) and V^ / (a 2 , r 2 ). These operators should be 
defined in the full CFT. Quite generally, the PSL(2, C) invariance of the sphere 1) 
allows to fix locations of the insertion points (<7i,Ti) and (02, T2); 2) it restricts the 
conformal weights (hf, hf) of the operator V^ f in terms of the ones of the initial state 
(hi, hi). Therefore the final state will be fixed automatically. However, we will still 
have a sum because the matching procedure described above will effectively yield a 
decomposition of \I// in terms of the spectrum of the IR CFT. 

What do we know about the initial string state \I/j? It should correspond to the 
Polyakov loop on the boundary gauge theory. Clearly it should be a winding w — 1 
state in the full CFT with zero transverse momentum p± = and zero Matsubara 
frequency k = 33 . Apart from these we cannot say much. In particular we do not 
know the precise form of the vertex operator that corresponds to this state, as we 
do not know the details of the full CFT. However, this ignorance will not affect our 
final result. 

Gauge fixing: One important complication in comparison to the QFT case above 
is that now we have two re-parametrization + one Weyl invariance on the world- 
sheet that should be gauge fixed. String paths are parametrized by the world-sheet 



33 It is also reasonable to assume that it corresponds to a state with no string excitation numbers 
but we will not assume this. 
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coordinates (a, r) with a ~ o + 2tt. In the path integral we can fix the two re- 
parametrizations by fixing the world-sheet metric to be of the form h a b = h a [,e UL with 
some reference metric h. The remaining freedom a^ is the Liouville mode, which can 
be left as unfixed. It is well-known that under quantum effects <jl becomes a space- 
like dimension and the target space becomes flat with one additional dimension plus 
a linear dilaton. It was further shown in J351 — see appendix O for a review — that for 



non-critical strings on flat d — 1 dimensional target-space to make sense at higher 
genera, one is forced to introduce an extra world-sheet field <ft that couples to the 
world-sheet Ricci scalar. Then, after fixing the re-paratmetrization invariance, the 
Liouville field o~l combines with (ft to produce two extra dimensions on the target 
space plus a linear dilaton field. Thus the end-result is non-critical string theory in 
d + 1 dimensional flat target-space with a linear dilaton $, which is exactly what 
we have in the IR in the model of section |3] for r^ 3> 1. Another option is to make 
sure that the CFT has vanishing total central charge. In this case o~l decouples 
and one has fixed the entire gauge symmetry on the world-sheet. The background 
becomes linear-dilaton in the range r>l for r^ — > oo. The two options are totally 
equivalent and for definiteness let us adopt the latter option. Then we start with 
d + 1 dimensions and we fix both re-parametrization and Weyl. 

At this point there are two options that one can choose to work with. One can 
either keep the ghosts that arise from the re-parametrization fixing or one can ignore 
them and include only the transverse string fluctuations in the canonical formalism 
and treat the propagation in the light- cone gauge. In the asymptotic linear-dilaton 
regime it is known that the ghosts exactly cancel the excitations of the string along 
the r and the x° directions in the linear dilaton background, just like the flat case 
[p9 ]. We will assume that this is true also in the more general case when we have 



the correction terms in ( |3.15|) . 



The calculation becomes more transparent in the light-cone gauge which can 
easily be generalized to the linear-dilaton background [[H], |29| . Here, one ignores the 



contribution of the re-parametrization ghosts and fixes the metric h a b = h a i ) e crL by 
hand. There is a residual freedom from the combination of diffeo-Weyl that leaves h 
invariant which can be fixed by, 

X + =p+ T + x + , X ± = —{X°±r). (5.34) 

v2 

X - is also fixed through the Virasoro constraint and one is left with only the trans- 
verse oscillators along X 1 . 

Calculation: For the purpose of identifying the contribution from the IR region we 
divide the propagation into two parts. The procedure we outlined for the field-theory 
case to separate the paths in the two different regions and to sew the propagators at 
the matching region r = r m has a direct generalization to the closed string propaga- 
tion in our background: One only has to 
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Figure 3: The string world-sheet that wraps the time-circle and hangs from the boundary 
to the horizon as the dual of the one-point function of the Polyakov loop. The Euclidean 
time is identified as indicated by the dashes. We insert a complete set of string states \ 
at an intermediate point r m chosen close enough to r^, for large r^, such that the IR CFT 
description holds for r > r m . 

1. replace the Hamiltonians in ( j5.29| )with the Virasoro generators H = Lq + Lq; 

2. extend the integration over r to a complex parameter w = r + ia whose imag- 
inary part couples to Lq — Lq, hence the integral over it produces the left-right 
matching of the on-shell states; 

3. replace the wave- functions \l/^ and ^^ of the IR and UV Hamiltonian with 
vertex operators in the CFTs. 

Formally we first decompose the path integral from r = to r m and r m to oo: 
(Psi i ,0\^ f ,oo) = (^ i ,0,0\^ f ,r h ,xf) = f vX^- a[x] - a[cjl] (5.35) 

Jbci,bcf 

= ! dX^(a) J VX^e- A[x] f VX^e'^, 

J Jbci,bcm Jbc m ,bcf 

where the boundary conditions are denned by the sets, 



ba = {X»(0, a) = X?{a) : X°(cr + 2vr) = X*{a) + -; X lfi = r = 0} (5.36) 



r 



bc f = {X^(oo, a) = Xf(a) : X f)0 = x f , r = r h }. 



(5.37) 
(5.38) 



Quantities with subscript refer to the center-of-mass positions in the first and the 
last lines. 
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The r-component of the intermediate string r m (a) has a center-of-mass piece 

/■27T J 

r(r m ,a)= / —r(T m} a)-\ = r -\ (5.39) 

Jo ^ 

Just as in the field theory computation above, we can use the freedom to choose r m 
to replace the integration over ro with an integration over r m by inserting 

1 = / dT m 5(r (r m ) - r m )J m , J~ l = -^\ ro{r)=Tm (5.40) 

inside the integral over X^(a). Now, we assume that only the paths with crossing 
number c m = 1 dominate the path integral. This is certainly the case for r m chosen 
close enough 34 to r^ for r/,>l. This means that, 

r(r) > r(r), r' > r, for t > r m , (5.41) 

where 

r(r)=J %r{r,a). (5.42) 

Therefore we achieved the division of path integrals from r^ = to tq = r m and 
r o = r m to r = r h . We will approximate the second class of paths by replacing the 
action with that of the linear-dilaton CFT: 

poo /» 

(*i, 0|* /; oo) w / dr m J m / dX'MPI(Q, T m )PI IR {r m , oo), (5.43) 

where the first path integral is 

PI(0,T m )= I Va L VX^e-^ x \ (5.44) 

Jbci,bc m 

with A the full world-sheet action and the second one is 



PI IR (r m ,oo)= VX^e- A ^. (5.45) 

Jbc' m ,bc f 

The primes in Q5.43] ) and ( ^.45| ) denote omission of the center-of-mass piece in r(T m ) 
as it is fixed to r m by ( |5.40|) . 

The approximation in ( |5.43|) is two-fold: First of all we approximate the action 
in the region r(r) > r m by the IR CFT: 



-1 pZTV /*00 

Am = da drVh h ab d a X^d b X u r]^ + Aa'v^R + 2a'b ab V a c b 



(5.46) 



34 We remind that the "crossing number" was defined in the field theory analogy above. This 
assumption can be lifted as explained in the field theory analogy above by modifying the procedure 
of dividing paths. This is an unimportant detail however, which has no effect on the final result. 
Here we just restrict the analysis to c m — 1 paths for simplicity. 
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where we idsplay explicitly the reparametrization ghosts. The proportionality factor 



in the dilaton in our case is given by (4.4). We consider the bosonic linear-dilaton 



theory for definiteness; the final result (in the semi-classical approximation that is 
to be defined below) is independent of the particular linear-dilaton CFT chosen. 

The second approximation is that we restrict the analysis to the c m = 1 paths. 
Both of these approximations become better as the point r m is chosen closer to r^. In 
fact, the second assumption is automatically satisfied in the light-cone gauge ( |5.81| ) 
where the role of r(r) is played by X + (r). 

Now, we focus on the second path integral PIir. In the canonical formalism this 
can be written as, 

PI IR (r m , oo) = Yl ( V x(Xm, r m )V*(X f , oo)) A IR ( X ), (5.47) 

XGWx 

where \ runs in the transverse Fock space of the linear-dilaton CFT and V x (X,t) 
denotes the vertex-operator for creating a closed string X at world-sheet time r in 
the x eigenstate of the Hamiltonian L + L |3l 



V X (X, t) = J ^ ^V L ( X , X l (t - a))V R ( X , X R (r - a)). (5.48) 

The eigenstates x °f our linear-dilaton CFT are labelled by the center-of-mass mo- 
menta in r and transverse directions p r , p±, the Matsubara frequency k and the left 
(right) oscillator numbers N (N). 

The propagator of a state x(Pr,P±, k, w, N, N) is then given by 

A«(x) = / ^|^o(x)-i/o(x)-i ; (5.49) 

J\z\<l \ Z \ 



where the Virasoro generators are given by ( |4.6|) . In the propagator ( |5.49 ), the 
integral over z projects the states on the mass-shell ( |4.9| ), ( |4.10|) . 



The IR path integral (|5.47|) contains the r^ dependence that we seek for, inside 



the vertex operator for the final state. It is contained in the center-of-mass position 
term in the r-direction, see (|5.36|). Let us make it explicit by factoring out 



V?(X f , oo) = e-^V*JX, oo), (5.50) 



where V contains no dependence on r^. On the other hand, the sum over x i n (|5-47 ) 



contain integrals over p r and p± and sums over k, w N and N. Noting also that the 
integral over z in (|5.49|) projects onto the mass-shell states (|4.9| ) and ( |4.10| ) one can 



directly perform the integral over p r in ( [5.47| ) and find, 

PI IR (r m , oo) = ^C x e-^^-H (5.51) 
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where the constant C x does not depend on r^ and p*(x) denotes the solution of the 
mass-shell condition (14. 91): 




m 2 



(5.52) 



w ( JV + J V_ 2 ) +pl + (2 ^ rf+ (_|_) . (5.53) 

Substituting (|5.51|) in (|5.43| ) we find that the entire r h dependence of the Polyakov 
loop becomes, 

(P[0])ocVC( X )e-^^, (5.54) 



where the states x have zero transverse momenta and C is some c-number whose 
value is independent of r^. This result can be thought of a direct generalization of 
the field theory analog in ( |5.33| ). The exponential term above is the analog of the 



wave-function \I/ X « and the coefficient C x above is the analog of the expression inside 
the brackets in ( j5.33| ). The important difference is that here we have an infinite sum 
over all possible on-shell states of the string. 

The expression inside the square-root in ( |5.52|) has the following behavior for the 
various physical states. For any state other than the tachyon (N = N = k = w = 0) 
it is larger than 1. For the tachyon it equals, 

m? 1 — d , 

1 + 4 = ^7 (5 ' 55) 

Thus it is negative for any non-zero spatial dimension d — 1. For a "winding tachyon" 
N = N = 0, k > 0, its sign is determined by the value of the temperature T. This 
latter case will prove important in the evaluation of the correlation length in section 
|5.4.2| . Here, the important point is to realize that ( |5.54[ ) is always dominated by the 
tachyonic ground state in the limit r — >■ r^ because all of the higher states result in 
bigger suppression in the exponential. 

For the tachyon on the other hand the square bracket in Q5.52 ) is oscillatory, 



thus it gives an imaginary contribution to the exponent in ( |5.54| ) and the modulus 
of (P[0]) is determined by the first term in (|5.52|) , 

lim |(P[0])|oce- mor \ (5.56) 

r h ->oo 



We can now translate the variable from r^ to the reduced temperature t using (|D.13| ) 
near the transition region, 

Consequently, we obtain the scaling of magnetization as, 

|M|oc|P[C]|oc*«. (5.58) 
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In particular, for the second order transition k = 2 one obtains the mean-field scaling, 

|M|octi (5.59) 

This provides a non-trivial check on the proposed duality. Although this is a stringy 
computation in principle, we observe that the scaling exponent is determined by the 
lowest lying fluctuation of the string that always correspond to the tachyon in a linear- 
dilaton CFT in arbitrary dimension. As this is a gravity mode, it is quite reasonable 
to expect mean-field scaling by the arguments in section |5.1| . The computation above 
confirms this expectation non-trivially. 

It can also be shown that all of the arguments that we made throughout the 
derivation readily extends to the other linear-dilaton CFTs, including the fermionic 
ones. This is because there always exists a tachyon in the spectrum for which the 
square-root in Q5.52D becomes imaginary, whereas it becomes real for all of the other 
states in the physical spectrum. This result can easily be understood intuitively: 
in the transition region T m T c where the linear dilaton CFT governs the scaling 
behavior, the existence of the tachyon signals instability, hence phase transition, 
and it is this tachyonic state which dominates and determines the scaling law of 
observables. In section ( |5. 4. 2 ) we provide another example of this phenomenon. 



There are various possible modifications of the mean-field result. First of all, it 
would be interesting to see whether going beyond the semi-classical approximation 
would modify the critical exponent. For this one has to sum over all of the string 
states instead of focusing only on the dominant tachyonic contribution. It will be very 
interesting to obtain corrections to mean-field scaling as a result of this computation. 
We hope to investigate this issue in the future. 

Second type of possible modification involves the 1/N corrections. In the calcu- 
lation above we assumed that the boundary value of the dilaton can be tuned strictly 
to zero, so that we can ignore gravitational interactions. One expects that the 1/N 
corrections modify the mean-field scaling as, 

|M|octrW a ). (5.60) 

5.4 The two-point function 



As explained in section |5.2| , the spin-spin correlation function is represented by an 



F-string solution on the gravity side, (|5.10| ) that wraps the Euclidean-time circle and 



is connected to two separate points on the boundary that we take as x and 0. As 



in (|5.3.1|) , we first compute this quantity classically as a warm-up exercise. We then 



generalize to the realistic case where one has to consider the full path integral. The 



classical computation was first carried out in the case of AdS black-holes in [[33] [34 
(see also 
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5.4.1 Warm-up: classical computation 

We will perform the computation in the black-hole phase (that corresponds to the 
super-fluid phase in the XY- model). The computation in the TG phase is very 
similar and as we are essentially interested in how the correlator scales near T c , the 
two results will yield very similar results. Although all of the NS-NS fields G^ u , B^ v 
and $ couple to the F-string, we show that the G-coupling yield the dominant term 
in App. p.2| . Thus one can replace the F-string action with the Nambu-Goto action 
in this section, unless specified otherwise. 



X1 



->■ r 



(a) 




■*■ r 



(b) 




(c) 



->■ r 



Figure 4: The classical string saddles that contribute the spin-spin two point function in 
the ordered phase of the ferromagnet. (a) The disconnected diagram. Strings are falling on 
the horizon. Its value is proportional to the square of the magnetization expectation value 
\M\ 2 . (b) Connected string diagram, (c) A bulk-mode exchanged between two disconnected 
strings. 



There are three string embeddings that contribute [35], see figure 



1. Disconnected diagram: This consists of two straight strings that connect to 
the boundary at points x\ — L = \x\ and x\ = and hang from the boundary 
to the horizon. The result is twice the NG action in ( |5.16|) , thus it is finite as 
long as T > T c . As T — > T c it diverges as in (|5.20|) , that corresponds to the fact 
that the magnetization vanishes smoothly as T — > T c . Thus the disconnected 
contribution yields, 



(m(x) ■ m(0)) dis = (\P\) 2 = \M\ 2 = finite, 



T>T r . 



(5.61) 



We are not interested in the actual value of the disconnected piece, it depends 
on the normalization of the Polyakov loop. As mentioned before this piece is 
absent in the TG phase. As it is disconnected, this piece is actually 0{g~ 2 ) 
enhanced with respected to the interesting connected contributions [^3] that 
should be handled separately. 

2. Connected diagram: This is a differentiable world-sheet whose end-points 
are connected to the boundary at points x and 0. To compute its on-shell 
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action we fix the gauge r = x%, a — Xq where x\ is the coordinate on the 
boundary on which the end-points of the string L and lie. The result is (see 
App. |D3): 

1 f r f e 2A s (r) 



s ™ = ^m^ / dr , ( 5 - 62 ) 

e 4A »W/(r) 

where 77 < r^ is the turning point of the string that corresponds to dr jdx\ = 
and e is a cut-off near the boundary 35 . A s denotes the scale factor of the metric 
in the string-frame, see (|5.17| ). 

The length between the end-points of the string L is given by, 

r r f r r f i 

L = 2 dxx = 2 dr =. (5.63) 

Je Je IJ7-S / e^s(r) f(r] 



Eqs. (|5.62 ) and (|5.63|) parametrically define the function S c ^q{L). One can 



easily show that L is monotonically increasing in r/. When the distance L 
reaches a certain value that corresponds to 77 = r^, the connected string solu- 
tion falls into the horizon, thus ceases to exist. Beyond this point this diagram 
gets replaced by the "exchange diagram" that we discuss below. 

We are eventually interested in the scaling of the spin-spin correlation function 
near T c . Therefore let us focus on the limit r/, > r/ > 1. In this limit T is 
very close to T c and at the same time L is large. One can easily compute the 
function S^q in this limit (see App. |D.2| ) and finds, 

where mr(t) defines an effective mass term, that stays finite at all T and whose 
value in the limit t — > is shown above. The ellipsis denote contributions that 
are sub-leading in L. Thus, we can identify the first contribution to the spin- 
spin correlator: 

(m(x) ■ m(0)) con ~ e~ mT L+ - (5.65) 

Comparing with ( |C.17D , one indeed finds qualitative agreement where m^ 
gives a finite contribution to the spin-spin correlation length £, that stays finite 
even at T c . On the other hand, we expect £ to diverge at T c . To see how this 
divergence arises one has to perform the full path integral computation that 



we turn in section B.4.2. 



35 For the sake of the discussion here, we do not need to renormalize the action by subtracting 
counter-terms. This can be done in a standard way, if desired. 
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3. Exchange diagram: When the curvature on the string world-sheet becomes 
strong, one should also take into account the fluctuations of the string, that 



are no-longer negligible. It was first observed in |35] that this contribution 
yields a crucial correction to the Polyakov correlator, that actually resolved a 
puzzle that was encountered in |33| |34|] : As the connected contribution above 



ceases to exist beyond rj = r^, one may naively think that the connected part 
of the correlator (P*(L)P(0)) vanishes beyond a certain value of L. This is in 
contradiction with a generic QFT as the connected piece of a generic correlator 
should be a convex function which smoothly decreases with increasing L. 

The missing contribution, in fact, comes from the world-sheet fluctuations that 
become crucial in the regions where the world-sheet curvature -PJ 2 -* becomes 
large. In App. |D.2j we show that R^ indeed becomes large near the horizon, 
hence another type of connected diagram that arise from world-sheet fluctu- 
ations become more dominant at r^ = r^. This diagram can be calculated 
in the limit L becomes large. In this limit, the contribution is given by the 
diagram that consists of two disconnected world-sheets connected by the ex- 
change of gravity modes. In the large L limit, the exchange is dominated by 
the gravity mode with lowest mass. Thus, in this limit exp —Sng is given by 
the propagator of the lowest mass gravity mode in d — 1 dimensions 36 : 

e~ SNG ~ Ld _ 3 , (5.66) 

where m min is the lowest mass bulk mode. 

Now, it is crucial to figure out which gravity modes contribute to which parts 
of the correlator. As explained in [^5| , in the gauge theory, only the CT-even 
modes couple to the real part of P[C], and CT-odd modes to the imaginary 
part. Here r denotes reflections in Euclidean time. The analogous statement 
in gravity is that only the CT + bulk modes are exchanged in the part of S^q 
that corresponds to (Re PRe P) and the CT~ bulk modes are exchanged in 
the part that corresponds to (Im PIm P). Using the identification (|5.13|) and 
(|5.14j) , we find that, 



(mi, (x) ■ m,| (0)} exc ~ e" m + L /L d ~ 3 , (5.67) 

(m ± (x) ■ m ± (0)) exc ~ e" m ~ L /L d ~ 3 , (5.68) 

where L = \x\, the m± are the lowest masses of the bulk modes in the CT ± 
channels and the result is valid in the limit L ^> 1. 

The NS-NS modes in the CT + channel with their J CT designations are, Goo 
(0++), Gf^ (2++), $ (0++), B id 1— , G tt (0++). One should solve the flue- 



36 The exchange mode propagates in d— 1 dimensions because the propagator is fixed at a certain 
value of r and it's compactified on the time-circle. 
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tuation eqs. on the BH background in order to figure out the lowest mass 
one. 

This is done in Appendix [E], where we showed that all of the masses that 
correspond to gravitational fluctuations are bounded from below and non-zero, 
i.e. there is a mass-gap in the CT + channel that is given by mo = a/V^/2, 
equation QE.6Q . Hence, in the large distance limit L ^> 1 the exchange diagram 
gives the contribution, 

(m,| (x) ■ m,| (0)) exc -jir^i L > !• ( 5 - 69 ) 

On the other hand, the NS-NS modes in the CT~ channel with their J CT des- 
ignations read, Giq (l + ~), B i0 (1~ + ) and B r0 CT + . Here, the crucial observation 
is that, as explained in section |2.2| , the zero-mode of the latter is nothing other 
than the field if) in fl2.4|) that was identified with the Goldstone model Thus the 
lowest lying mass in the CT~ channel is the zero mode of the B r0 field and it 
is zero: m_ = 0. Thus we find no attenuation term in the corresponding part 
of the spin correlation function: 

— m-(t) L j 

(m±(x) ■ m ± (0)) exc Ld _ 3 = jj-g, (5.70) 

where the result is valid for large L = \x\, and at any temperature T > T c 
which corresponds to the ordered phase in the corresponding spin-system (we 
recall that the temperature on the gravity side and the spin-model side are 
inversely related to each other). 

Combining ( |5.6ip , ( |5.65| ), ( |5.69| ) and ( |5.70| ) we arrive at the total result for the 
spin correlation function (in the large L limit): 

p -m,Q L 

(m(x) ■ m(0)> ~ M 2 + c x e~ mT L+ '" + c 2 -^ zr + -^, (5.71) 

where q are some constants. 

Comparison with the mean-field spin-model result ( |C.17[ ) shows perfect qualita- 
tive agreement for temperatures T > T c (which corresponds to the low T regime of 
the super-fluid). 

However, we also observe that the classical computation fails to reproduce a very 
crucial feature of the spin-model at T c . Namely, the longitudinal component of the 
two-point function should in fact have a vanishing exponent as T — y T c : 

g— m||(t) L+— 

(rh\\(x) • m||(0)) ~ 71^3 ' with m^(t)^t u , as t — > 0. (5-72) 

This corresponds to the fact that the longitudinal correlation length also diverges at 
T c . On the other hand, in our classical string computation we found an exponent my 
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bound from below as my — > min(m , mx) where trq and itit are given by ( |E.6|) and 
( |5.64 ). We will argue below that the full path integral computation of the two-point 



function yields the desired result. 

5.4.2 Semi-classical computation 

Now, we look at the more generic situation when the assumption £/£ s ^> 1 fails — as 
expected in non-critical string theory — and the string path integral that corresponds 
to the two-point function (P*(x)P(0)) = (fh(x) -m(0)) is given by the full string path 
integral. As in section (|5.3.2| ) we shall calculate this quantity in the semi-classical 
approximation where we focus on the dominant contribution of only the lowest lying 
string states at levels N = and N = 1. 

One can again classify the string paths according to the three classes as in figure 



The contribution of disconnected paths will be just as in section 5.3.2. In the 



TG phase they vanish because the area of the string paths are infinite and in the 



BH phase they yield the square of the one-point function found in section |5.3.2j . In 
addition one also have to consider the disconnected paths corrected by bulk-exchange 
diagrams as in figure U We shall consider the contribution of these latter diagrams 
in the end of this section. First we focus on the connected string paths, see figure [5[ 
The connected path integral of the string is given by summing over all paths the 
string can travel between the space-time points (r, Xi,x±) = (0, 0, 0) and (r, Xi,x±) = 
(0, x, 0). Here x± denote the coordinates transverse to r and x. We denote these 
points as Ij„ and I out respectively. As displayed in figure [| for the connected classical 
saddle, in the limit L ~^> 1 these paths can naturally be divided into three parts, see 
figure [| 

1. The paths between the space-time points /„ = (0, 0, 0) and /;. 

2. The paths between the points Jj and //. 

3. The paths between the points If and I out = (0, x, 0). 

This division of path integrals into separate regions in space-time is a non-trivial 



operation that is described at length in section 5.3.2. Here we will not go through 



the same derivation again but only highlight the computation. 



Just as in section |5.3.2| we write formally divide the full path integral by inserting 
complete set of states at Jj and If as, 

(P*(x)P(0)) conn « fdlidlf Y, Hx,*i,Ii)±iR{x,Ii,h)F'{x,*t,l!), (5-73) 

where \l/j and \1// denote the initial and final string-wave functions and the sum is over 
the physical string states in the Fock-space of the string. The function J-"(x, *&i, It) 
denotes the overlap of the string state ^ at the point I in and the state \ a t the 
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point Ii. Similarly J-(x, *&f, If) denotes the overlap of the string state \l// at the 
point I ou t and the state \ a t the point //. The approximation in ( |5.73|) is due to 
the assumption that the propagation in the intermediate paths between Jj and If 



are governed by the IR CFT as in section 5.3.2 . The properties of this IR CFT are 



specified in section [4~2. 




IV IRCFT 



Figure 5: Quantum mechanical calculation of the spin-spin correlator. The string paths 
are divided into three separate paths by insertion of complete set of string states at It and 
If. For large L and r^ intermediate paths are governed by the IR CFT. 



The propagation of the string in the paths 1st and 3rd class depend on the full 
CFT, hence we cannot calculate. However, in the limit r^ — > oo (T — > T c ) and for 
large values of L, the dependence of the two-point function on L is determined by 
the paths in the 2nd class where the propagation takes place in the region where the 
IR CFT can safely be approximated by the linear-dilaton CFT of section |4.2| . Let us 
therefore concentrate on the evaluation of this part of the path integral. 



As discussed in section |5.4.1| , one has to consider the two parts of the correlator 
separately: 



(P*(x)P(0)) c 



(ReP(x)ReP(0)) cmn + (ImP(x)ImP(0)) c 



(5.74) 



Only the CT + (CT~) string states contribute in the real (imaginary) parts. First of 
all, we are interested in determining the scaling of the longitudinal correlation length 
£(£) near T c . For this purpose we focus on the real part in what follows and return 
to the imaginary part in the end. 

Thermal gas phase: The propagation of closed string state \ from a point x = 
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to 37 x = L in the linear-dilaton CFT is given by 

f_z 

'\z\<l \ z \ 

where p x {x) denotes the momentum of the state x m the x direction. The exponential 



A 7R (X, 0, L) oc / d Px e - iMx)L J /oW-i^oIx)-^ (5.75) 



factor in ( |5.75|) arise from the center-of mass part of the vertex operator insertions 



at x = and x = L, see section 5.3.2 



The integral over z restricts to the mass-shell states ( |4.9|) and ( 4.10 ) on which 
the momentum p x becomes, 

i 
Pi = ~i(^ (N + N - 2) +p 2 ± +p 2 r + 2im p T + (2nkT) 2 + (3^7) J * , (5.76) 

where p± now denotes the momenta transverse to the x-direction in the spatial di- 
mensions. As mentioned before the boundary condition for the closed string restricts 
only to the states with winding number w — 1. This cannot change along the prop- 
agation in the limit g s — > 0, unless the string paths fall onto the horizon, that does 
not happen for the connected string diagram. We therefore focus on the case w — 1. 
One should of course satisfy the level-matching condition ( f4.10|) , N — N + k = 0. 



In order to read off the scaling of the correlation length from the real part of 
the propagator then one has to consider the contribution of all CT + states that has 
winding w = 1. The propagation amplitude of a state x i n this region is given by 

A IR ( X , 0, L) oc J d Pr d d - 2 P± e-^ L . (5.77) 



The momenta p r and p± in (|5.76|) are to be integrated over after substituting ( [5.751) 



in (|5.73|) . This can easily be seen to produce a factor of L d ~ 3 . The correlation 
length then can be obtained by state Xmin for which the expression — ip* x in ( |5.76| ) 
for p r = p± = is minimum. 



£(*) = ¥x(.Xmin) 



(5.78) 

p ±=Pr =0 



Consequently, in the thermal gas (disordered) phase the correlation length is 
given by 

7 2 / ~ _\ / 1 



6J= |(iV + J V-2) + ( 



2nTa 




(5.79) 



The minimum CT + state is clearly given by the tachyon N = N = 0. Level matching 
then sets k = and one obtains 




&g= -77+U^-7 • (5-80) 



.",7 



We drop the subscript in x\ for notational convenience. 
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For an arbitrary temperature this is a positive number. However we observe that 
it vanishes precisely at the temperature when the winding tachyon mode 38 becomes 
massless: 

T » = m: (5 - 81) 

For critical strings in flat space-time, this special radius was shown to correspond to 
the Hagedorn temperature of the string. Furthermore |17| argued that upon turning 
on an infinitesimal g s this point becomes a first order phase transition in the string 
partition function. 

Here, similar arguments combined with our finding section |3] that the linear- 
dilaton with mild subleading corrections has a continuous transition would imply 
that, in the case of a linear-dilaton background rather than the flat space, the same 
point indicates a continuous transition into formation of black-holes. In the limit of 
vanishing g s , this argument thus determines the phase transition temperature as, 



T, 



n = ^- (5-82) 



From Q5.80|) it is also clear that one obtains mean-field scaling for the correlation 
length, 



i 



£r G (t) -► -^=|*r*, ast->0, (5.83) 

in the bosonic linear-dilaton CFT. (We recall the definition t = (T — T c )/T c .) 

In case of the fermionic CFT a similar calculation in the NS-NS sector yields a 
limiting temperature that we again propose to coincide with the continuous Hawking- 
Page transition: 

T S H = T C = — i , (5.84) 

9s =o 2\/27r4 

with similar mean-field scaling of the correlation function. On can consider variants 

of the CFT by applying non-standard boundary conditions on the r-direction, but 

the result does not change. As long as there exists a winding tachyon one obtains 

mean-field scaling provided that one identifies the transition temperature with the 

Hagedorn temperature. 

Our result for the real part of the connected contribution is then summarized as, 



e ZtgW 
(ReP(x)ReP(0)) TG ex -jp^-, (5.85) 

where the correlation length near T c is given by fl5.83| ). 

Now, we consider the imaginary part in ( p.74j ). The only difference is that 
now we have to sum over the intermediate states with CT~ quantum numbers and 



38 This is a misnomer as the "winding tachyon" is actually massive for an arbitrary temperature 
above Hagedorn. 
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with winding w — 1. The lowest lying such state is given by the NS-NS two- form 
fluctuations in the N = N = 1 level. For the w = 1 state, level matching again 
requires k = 0. It becomes clear that this contribution is subdominant with respect 
to (|5.85|) in the large L limit: 

L 

(ImP(x)ImP(0)) TG ex ^J, (5.86) 

where £o(£) asymptotes to a constant at the transition £,o(t) — > 2-KT c d, 2 s . Using the 
value of T c obtained in (|5.82 ) this constant is, 



&(<) -> Y' flsM0 ' ( 5 - 87 ) 

Recalling that the disconnected diagrams in the thermal gas phase vanish because 
they correspond to paths with infinite area, our final result then is, 



e Stgw 
(P*(x)P(0)) TG oc-^ r ^ + ---, (5.88) 

with correlation length near T c is given by ( |5.83| ). The ellipsis denote contributions 
from higher mass states in the spectrum such as (|5.85 ). 



As an aside we also observe that the classical computation of the previous section 
gave "almost" the right answer except that it missed the tachyonic contribution. One 
can check this by observing that the contribution from the 1st level mass states is 
precisely the same as the classical result ( |5.65| ), ( |5.64j ). We already calculated this in 



(|5.86 ) for the B-field, and one has the same result for the other states at the same 



level, namely the dilaton and the graviton fluctuations. 

Black-hole phase: Let us now consider the Euclidean black-hole. Here as well, 
the winding will be protected along the propagation of the closed string as long as the 
string does not fall on the horizon. This is indeed the case for the connected paths. 
The calculation of the connected paths is very similar to the thermal gas case above, 
with one important difference: here T > T c and a naive application of ( |5.79|) would 
give rise to an imaginary correlation length. There is another important difference 
though: the BH becomes linear-dilaton strictly in the limit r^ — > oo and for any 
temperature less than T c we expect additional contributions to the mass spectra. 

From the space-time point of view there is a finite horizon for finite r^ and 
the computation of the mass spectrum — as fluctuations in space-time follow from 
applying normalizable boundary condition at the horizon. This generally gives rise 
to a discrete spectrum and one expects a correction to the mass spectra which is 
supposed to vanish only in the strict r^ — > oo limit, due to the presence of the horizon. 
The only invariant quantity that would be a candidate for such a correction term 
then is the black-hole mass itlbh- Consequently we expect that the mass spectra be 
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shifted positively by a term proportional to uibh near r^ < oo. On the other hand 



it is easy to see that the ADM mass of the black-hole [2(| is proportional to the 



string- frame Ricci scalar that was computed in [|15| and it was found that, 

m 2 BH ocR s oc e^H' 2 oc U~ 2 , (5.89) 



where the last relation follows from ( |D.12| ). Consequently, the relation ( |5.80|) should 



be modified in the limit r^ — > oo in Euclidean black-hole as, 

where qy, > is some constant that we cannot determine unfortunately 39 . Thus we 
find, 

tBH^Vq^&t*, ast^0 + . (5.91) 

It should be checked that c^ > 8 for consistency of the picture we present here. 
However we do not have any means to check this at present. 

The result ( |5.91| ) summarizes the dominant contribution in the large L limit to 
the connected string paths. As the tachyon is a CT + state ( |5.91|) yields the leading 
term in the real part of the connected two-point function: 



e (■bhWs 
(ReP(x)ReP(0)) conn oc Ld3 , (5.92) 

with ( CT) . 

How about the imaginary part? Just as in the thermal gas case above the leading 
connected contribution to the imaginary part is given by the w = 1 NS-NS two-form 
with the same answer ([5.851): 



L 

(ImP(x)ImP(0)) conn oc ^|J, (5.93) 

with £o behaving as ( |5.87] ) in the limit t — > 0. 

Another crucial difference between the black- hole and the thermal-gas is that, 
for a finite horizon (any r^ other than strict r^ = oo) there are disconnected string 



paths and similarly exchange diagrams of the sort we described in section |5.4.1| . The 
two disconnected string paths that fall into the horizon as in figure (| give square of 
the one-point function that we already calculated in section |5.3.2j : 



(P(x)P(0)) dia ex |M| 2 , (5.94) 



39 One can try to obtain an effective action for the winding tachyon on the Euclidean black- hole 
and determine its spectrum by applying normalizability both in the UV and near horizon. However, 
there will still be undetermined coefficients in the effective action and the constant would not be 
determined. The only way to determine it is to construct the small BH as a marginal deformation 
of the linear-dilaton CFT and obtain the exact spectrum. 
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where \M\ is given by ( |5.58| ) in the limit t — ¥ 0. This contribution is order 0(g~ 2 ) 
[p5|] and it should be treated separately. 

There are also exchange diagram contributions of the form figure |j. These are 
of the same order in g s expansion as the connected diagrams and their contribution 
is simply given by the sum over all possible bulk modes that can couple to the 
disconnected string paths between the boundary and the horizon, see figure |j. As 
the disconnected string paths do fall on the horizon, the string with a winding w — 1 
boundary condition at r = can unwind at r^ even when g s — > 0. Therefore 
the winding number is not conserved in the exchange diagrams in the black-hole 
phase and we should include the contribution from the non-winding w = states. 
Contribution of the exchange- diagram in figure ^ is then given by the propagator of 
a bulk-mode with a d — 1 dimensional mass m\: 

(P(x)*P(0)) exCtm± oc L-V-Qe-y/V^y+mlL^ (5 Q5) 

The dominant contribution is always given by the k = modes, hence we consider 
the case k = in what follows. Again we divide this into the real and imaginary 
parts which receive contributions from the CT + and CT~ bulk states respectively. 
We first focus on the imaginary part. The lowest mass CT~ contribution is the 
fluctuations of the B-field that are massless from the d—1 dimensional point of view: 
m± = 0. We refer to appendix [E] for a derivation. Consequently we obtain 

(ImP(x)ImP(0)) exc ex jj^ + ■■■ (5.96) 

This is the transverse part of the spin-spin two point function in the superfluid 
(ordered) phase and the fluctuation of the B-field in this phase is identified with 
the Goldstone mode. There is no analogous term in the thermal gas (disordered) 
phase because the disconnected string paths that the exchange mode couple to, have 
infinite area and this contribution vanishes. 

Now we consider the real part. The lowest mass CT + bulk mode is the tachyon 
and we have to compute its mass mi in the d—1 dimensional point of view. This 
is a dangerous mode because contribution of a negative mass state to the two-point 
function implies non-unitarity in the dual spin-system. The analysis of the bulk spec- 
trum can easily be turned into a Strum-Liouville eigenvalue problem, see appendix 
[E]. In this appendix we show that, in the interesting cases of two and three spatial 
dimensions, the tachyonic mode can be avoided only when k — 2. Quite conveniently 
the case k = 2 corresponds to a second order transition — the most interesting case! 
We conclude that consistency of the entire analysis can only be established for second 
order phase transitions in 2D or 3D (spatial). 

Let us focus on this interesting case and denote the minimum value of the d — 1 
dimensional tachyon spectrum as m g . Then, the dominant contribution to the real 
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part of the ( |5.95 ) is 

(ReP(x)ReP(0)) exc oc L" (d " 3) e~ m! > L . (5.97) 



Comparison of (|5.92 ) and ( |5.97 ) shows that the former dominates in the large L limit 



and near the transition t « 0. One then finds for the final result, 

(ReP(x)ReP(0)} BH oc L'^^e'^ 1, + . . . ( 5 . 98 ) 

with mean-field scaling of the correlation length ( |5.91| ). In the imaginary part, the 
comparison of ( |5.93| ) and (|5.96| ) shows that the latter dominates for large L and the 



final result is 

(IrnP{x)ImP{0)) B H oc ^j— + • • • (5.99) 



One should add to these two, the disconnected contribution ( 5.94J) . Comparison with 



the spin-model result ( |C.19| ), using the identifications ( |5.13| ) and ( |5.14| ) shows perfect 



agreement:/n the limit of weak gravitational interactions g s — > ; we obtained exactly 
the same correlator with the XY model, with mean-field scaling exponents rj = 0, 
/3 = v = 1/2. How the mean-field scaling can be altered will be discussed in section 

5.5 D-strings and vortices 

In the original picture of [J36J1 [l37l and their subsequent generalizations, the relevant 
field theory on the boundary is SU(N) gauge theory and the Wilson loop in question 
traces the path of an "electrically charged" fundamental field in the theory, i.e. the 
"electric quarks" . The motivation to relate the Wilson loops with the open strings is 
obvious in the D-brane picture, where the fundamental strings couple to electrically 
charged fields on the D-brane. Similar considerations also suggest that the Wilson 



loop that traces a "magnetically charged" particle, i.e. the 't Hooft loop |p9|, should 



be related to the Dl branes [|13|]. Indeed, this picture is very suggestive and the 
various computations in the context of QCD-like holographic models confirm the 
field theoretic expectations |40| . 



In the LGT-spin model equivalence, the Polyakov loops are mapped onto the 
spin operators on the spin-model side. Similarly, it is very suggestive to relate the 
vortices, with the 't Hooft loops that are dual to the Dl branes on the gravity side. 
Thus we propose that the Dl brane configurations are the right tools to probe the 
vortex dynamics in spin-models. 

In this section we perform some basic checks with the D-strings and show that 
indeed one obtains the expected qualitative behavior of the correlation functions. We 
reviewed the expected spin-model result in the case of two-dimensions at the end of 
appendix |Cj. We denote the operator that creates a vortex that is localized at point 
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x by v(x) and similarly the anti- vortex by v(x). The proposal is the following chain 
of dualities: 

(v(x)) ±+ (tr P e-f Ao ) ^e~ SD1 (5.100) 

where the second object is the 't Hooft-Polyakov loop in the dual gauge theory and 
Soi denotes the on-shell value of the dual D-string configuration. 

5.5.1 One-point function 

Both in this section and in the next we consider the classical calculation of the D- 



strings. A semi-classical calculation in the sense of section |5.3.2| and |5.4.2| is saved 
for future work. 

First of all, we note that the vortex charge in the spin-model should be dual to 
the D-brane charge on the string theory side. As mentioned at the end of section 0, 
the total charge of vortices in a configuration in 2D should vanish. The equivalent 
statement on the gravity side would be that the total number of .Dls and Dls which 
wrap the sub-manifold spanned by r, t coordinates should be equal in D = 4. This 
is indeed the case in D = 4 as the Dls are charged, and the gauge-field that couple 
them in the flat transverse space has a log-divergence. This means that when one 
considers an ensemble of Dl and Dls, the configurations with non-equal numbers 
of the two species have vanishing Boltzman weights in the partition function. This 
provides the first basic check in favor of associating the vortices in the spin-model by 
the D-strings. Note that the argument applies equally-well when the target space is 
TG or BH geometry. Thus we obtain, 

(v(x))tg = (v(x)) TG = (v(x)) BH = (v(x)) BH = 0, for d- 1 = 2. (5.101) 

In higher dimensions the argument above does not apply and one can have config- 
urations with non-trivial charge. 40 Therefore, in higher-dimensions the expectation 
value is determined by evaluation of the on-shell D-string action. 

The boundary condition for a single D-string is just as in the case of F-strings, 
above: it ends on the boundary at point x and wraps around the time-circle. The 
action is given by, 

S m = -7\ ! d 2 ae-* (det[h ab + b ab })^ , (5.102) 

where T\ is the D-string tension and we defined, 

h ab = g%d a X»d b X u , b ab = B^d a X»d b X». (5.103) 



40 The dual analogous objects to higher dimensional vortices (vortex lines, planes, monopolcs etc) 
in higher dimensions are the Dp branes with p = d — 2. 
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Here g s is the string-frame metric (|5.17 ) in the case of the BH phase. In the TG 



phase the metric is given by the replacement / — > 1 and A — > Aq in Q5.17D / 

It is straightforward to find the on-shell action that corresponds to a single D- 
string hanging from the boundary to some point rf. from ( 5.102|) , (|5.103|) and ( 5.17|) : 



S m = -%± ! ' dre^Ve^+b 2 , (5.104) 

where b is the constant value of the B-field on the r-t subspace 42 : 

b = B r0 = const. (5.105) 

We compute the on-shell action separately on the TG and the BH geometries for 
d — 1 > 2 in App. |L).3| . In the BH case string hangs down to the horizon rf = r^, 
whereas in the TG case rj = oo. The result is that the on-shell action is finite both 
on the TG and one the BH geometry. 

(v(x)) TG ^ 0, (v(x)) BH ^ 0, for d - 1 > 2. (5.106) 

This is unlike the F-string case which diverges on the TG geometry, and yields 
(m) = 0. The reason for finiteness here is clear in (|5.104|) : The only potential 



divergence 43 would be in the TG case where the upper limit of the integration is 
rj = oo. The factor A s — > in that limit, however the action is still finite because 
exponential suppression provided by the exp(— $) term in the action (|5.102|) ; note 
that $ grows linearly near the singularity, ( |3.16|) . 

This provides a second non-trivial check on the proposal. It matches the dual 
statement in the XY model is that the vortices play no role in determining the phase 
of the system for higher than two-spatial dimensions. We will thus consider the case 
of d — 1 = 2 below. 

5.5.2 Two-point function 

As the one-point function vanishes, the first non-trivial object is the two-point func- 
tion (v(x)v(0)) in d — 1 = 2. This is dual to a connected Dl — D\ configuration, 



completely analogous to the F-string case described in section |5^J. The boundary 



conditions are exactly the same as in that case. Here too, we confine our interest in 



41 The D-string also couples to the gauge field on it. In fact the only gauge-invariant combination 
(under "big" transformations A^ — > A^ + a M ) is of the form b a b + fab where f a b is the pull-back of 
the gauge- field strength on the D-brane. One can make the gauge choice A^ = 0. This choice does 
not affect the discussion below. 

42 As mentioned before we take the i?-field to be either pure gauge or constant. Here we entertain 
the second possibility. 

43 The standard UV divergence e — > can easily be cured by adding counter-terms as in App. 
fl2. 
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the classical computation in order to see whether the association of D-strings with 
vortices of the spin-systems pass the basic qualitative tests. 

The computation is non-trivial and it is presented in App. p.3| . The result is as 
follows: 

1. The thermal gas: We denote the difference between the end-points of the 



Dl — Dl configuration as L = \x\. As in section 3I, one can consider three 
contributions: 

1) The disconnected D-string contribution: This is given by a discon- 
nected Dl brane and a -Dl brane hanging from the boundary and extends up 
to the singularity at r = 00. This contribution is dual to \(v)\ 2 and vanishes in 
d — 1 = 2, for the reason described above. It is finite for d — 1 > 2. 

2) The connected D-string contribution: In App. p.3| , we show that there 
is a maximum value L max that is independent ofT and above which, there exists 
no connected D-string solution. Therefore this configuration is replaced by the 
exchange diagram directly analogous to the diagram described in section |5.4| . 

3) The exchange diagram: For L > L max this is the only non-trivial con- 
tribution. It is given by a Dl and a Dl connected by exchange of bulk modes 
that couple to the Dls. For large L it is proportional to the propagator of the 
lowest mass bulk mode exp(— mu(T)L). As in section |5]4|, this mass is bounded 
from below, iri£,{T) > for all T < T c as there are no massless modes. 

We learn that the total result on the TG geometry for large L is, 

(v{x)v{0)) TG ~ e~ mDL L>1, d-l = 2, (5.107) 

where m D should be determined by a study of fluctuations around the TG 
geometry. On the spin-model side this means that there is a finite correlation 
length between vortices and anti-vortices in the high-T phase. This is in accord 
with the expectation that one obtains a plasma of vortices and anti-vortices in 
the high T phase of the 2D XY-model. 

2. The black-hole: The computation on the BH geometry is completely analogous 
to the TG case above: 1) The disconnected configuration vanishes in d — 1 = 2 
and is non-zero for higher dimensions. 2) There exists a L max (T) above which 
the connected D-string configuration does not exist. This time, however L max 
is a function of T. One can show (see App p.3|) that L max is finite for any T. 
Thus if we are interested in the qualitative result for large L, then it is again 
determined by the exchange diagram. The exchange diagram in the BH case 
differs than the TG case above, in that, there exists a massless bulk excitation 
that couples to the D-string. It is given by the zero mode of the 5-field as in 
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section 5A. Thus one finds, 

(v{x)v(0)) BH ~ e°( logL ) L > 1, d - 1 = 2, (5.108) 

which gives a power-law with a T-dependent power. In order to determine 
the power, one should calculate the 0(logL) terms in the exchange diagram. 
We postpone this computation to future work, and content ourselves with the 
qualitative result ( |5.108|) . 



All in all, comparison of ( |5.107f ) and (|5.108|) with (|C.30|) provides a non-trivial 
check on the proposal. 

5.6 Vanishing of the second sound 

As reviewed in the end of appendix |C|, the speed of sound that is associated with 
the phase fluctuations ip should vanish linearly in the mean-field approximation as 
T — y T c in the ordered phase. The effective action for this Goldstone mode can 
directly be obtained from the gravity action ( |3.1| ) because it maps onto the B-field 
on the gravity side. We recall that in the large N limit, the Landau action is given 
by the on-shell gravity action ( |5.3|) . Therefore to obtain the effective action of the 
Goldstone mode we should consider quadratic fluctuations of the B-field around the 
on-shell value. We recall that the in the BH phase the phase of the mean-field 
acquires an expectation value given by 



$= / B (5.109) 

J M 

where M is the submanifold of the blackhole spanned by the (r, xq) coordinates. 
Goldstone mode corresponds to x-dependendent fluctuations ip — y ip + 6ip(x). It is 
crucial that this fluctuation cannot be gauged away by a gauge transformation of the 
B-field of the form 

B r0 -> B r0 + <9 r £ - d £ r . (5.110) 

The shift in ip corresponds to a "big" gauge transformation, it is a topological sym- 
metry in the problem |T3| 44 . 



A technical but important point is to regulate the divergence of the on-shell 
action near boundary by subtracting a counter-term action. One can obtain this 



counter-term action by the standard methods of holographic renomalization |^TJ but 
an easier way — that is equivalent for our purposes here — is to subtract the on-shell 
thermal gas action. Consequently the Landau functional in the large N limit will be 
given by 

F L ex AA = Abh - A TG (5.111) 



44 



I am grateful to Sean Hartnoll for a discussion on the various issues in this section. 
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This is valid only in the black-hole phase i.e. T > T c which corresponds to the low 
T superfluid phase of the spin system. The energy of the Goldstone mode should 
be computed in reference to F L that corresponds to the ground state energy of the 
spin-model: 

F L (il> + <ty>) - F L ty)) = 8^F{iP) oc AA{ip + <ty>) - AA(ip). (5.112) 

Thus, we substitute 

B^ = B r0 + 8ip (x), B r0 = pure gauge. (5.113) 

in ( j3.1|) and obtain, 

fyi^(V) oc AA + C^(r h ) f d^xdMdjdiP, (5.114) 



where g denotes the BH metric ( |3.6| ) and g denotes the TG metric ( |3.5| ) and A*4. 
is the difference between the parts of the on-shell values of the action that do not 
depend on B. For a second order transition this piece vanishes at criticality as 

AAq ~ t 2 . (5.115) 

The coefficient C^ in ( |5.114[ ) on the other hand contains the desired information on 
the sound speed (and the energy of the Goldstone mode). It is given by 

r n a , r°° a , 

C4r h )= dr^e-^gVg rr g 00 - dr ^e ^*° ■<# ' g% 'g^ . (5.116) 

Jo Jo 

On the spin model side, this is proportional to the kinetic term (?, J VSip ■ V8ipd d ~ 1 x, 
hence the speed of sound for the Goldstone mode is given by the coefficient C^(rh)- 
One finds, 

P r h poo 

c$oc(Vr h )=/ dre^*+ (d - 5)A - / dre-^ 0+id - 5)A ° (5.117) 

Jo Jo 

First of all, we observe that the energy of the Goldstone mode 

E^ oc C,p J d d - 1 xd i 8^d j 8i), (5.118) 

is finite for any rn- This is a crucial requirement to be able to associate the fluctua- 
tions of the B-field with the Goldstone mode. 

Secondly, we find that c£ indeed vanishes as t — > {r h — > oo) because the BH 
background (A, $) asymptotes to the TG background (A , $ ) i n the limit r h — > oo 
where the BH mass vanishes. This is a nice check already because it also confirms that 
c^ vanishes only in a continuous transition which requires that the saddle solutions 
coalesce in the transition region. 
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A more non-trivial check however is to see whether the scaling exponent for the 
vanishing rate is indeed the one expected from the mean-field scaling, i.e. whether 
cS, ~ t or not. We recall that 45 mean-field scaling is expected whenever an operator 
in the dual field theory is related to the fluctuations of the bulk action on the gravity 
side, in the limit of weak gravitational interactions g s — > 0. Using the asymptotics 
of the BH background functions in fl3.16| ) and ( |3.17[ ) we find that the contributions 
from A and $ in the exponent conspire nicely to produce 



c\ ex 



/ dre~ VV ^ rh oce~ VV ^ rh , r h > 1. (5.119) 

Finally, use of ( p,13|) yields, as t — > 0, 

{2_ 
t«, case i, 

2(± yi (5-120) 

e vc; , case ii. 

The case of a second order transition corresponds to k = 2 in case i, see equation 
( p. 27 ) and we indeed find the expected mean-field behavior: 



qj ~ t, as t -> 0. (5.121) 

We note that the precise form of the kinetic term for the B-field in ( |3.1|) is crucial 
in reproducing the desired behavior. This form stems from the non-critical string 
action in d + 1 dimensions after a Weyl transformation to the Einstein frame. As the 
physical results should be independent of the frame, one can of course produce the 
same result directly in the string frame. 

In fact, the calculation is more transparent in the string frame. The asymptotics 
as rh — > oo are given by the linear-dilaton background, where the string-frame metric 



becomes flat RIB" | 



9s,ufi = e d -i*^ ->■ Sp,. (5.122) 
On the other hand the gravity in the string frame is given by 

A a f d d+1 x^ge- 2 *(dB) 2 . (5.123) 

As the metric becomes flat, the scaling of cL is entirely determined by the factor 
e~ 2 * above: 

cjoce" 2 ^, $ fc »l, (5.124) 

where $^ is the value of the dilaton on the horizon. Use of (|D.12j) now produces the 



same result as in ( |5.120|) . 



-15 c 



See section 5.1 for a discussion on this issue. 
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It is also easier to investigate possible a' corrections in the string frame. Higher 
derivative corrections to the B-field can be schematically represented as 

CO „ 

A« ^ c klm Je- 2 *(dB) 2k R«\d<S>) 2m , (5.125) 

k,l,m=0 

where cu m are some unknown constants — that are supposed to be determined by 
the world-sheet sigma model — and R^ l > represents higher curvature invariants, e.g. 
RW = 1 and R w = R, R {2) ~ R 2 + R a p jS R a ^ 5 etc. Sum over non-trivial cross 
contractions between dB, R and <i$ terms are also implied to be included in this 
schematic expression. 

The linear-dilaton solution in the asymptotic region is a' exact, therefore the 
form of the background functions $ and g StlijV are not subject to a' corrections in the 
far IR. All of the curvature invariants in the string frame vanish as shown in |15| . On 
the other hand only the term k — 1 above can contribute to the quadratic term in 
Sip because the leading order piece in B is constant, see equation (|5.113|) . Therefore 
one finds, 

oo „ 

5 2 A S ^J2 Cm e- 2 *W) 2 (d$) 2m (5.126) 

m=0 •* 

Finally, we note that all of the dilaton invariants also go over to a constant in the 
linear-dilaton background, hence one still obtains (|5.117 ) with a renormalized overall 



coefficient. We conclude that we do not expect the a' corrections change the linear 
scaling in ( J5.121 ). 



6. A proposal for gravity-spin model correspondence in the 
general case 

Here, we would like to return the discussion of section 2 and promote the gravity-spin 
model duality that we advocated in the case of U(l) models to the general case. 

We are interested in employing gravitational techniques to learn about the dy- 
namics of the spin model around the transition. We want to map the spin model 
to a lattice gauge theory, which than will be related to a gravitational background 
in the continuum limit. Unlike the derivation of the spin-model from the LGT as 
reviewed in Appendix [A], the opposite map from the spin model to the gauge theory 
is non-trivial. There are two sources of complication: 

1. Non-uniqueness: The map may be non-unique. Clearly, there may be many 
gauge theories that share the same center symmetry. First of all, this may be 
due to the fact that the center symmetry of different gauge groups may be the 
same. As an extreme example, the centers of SU(2), Sp(N) with arbitrary N, 
SO(N) with odd N and E(7) are all isomorphic to Z 2 . Thus for example, the 
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critical phenomena (if exists) in any of these theories should be described by the 
universality class of the Ising model in d dimensions. Secondly, the deformation 
of the pure gauge theory by addition of any adjoint matter leaves the center 
symmetry intact. Thus, generally, the equivalence maps a spin model with 
symmetry C to a set of LGTs with various gauge groups G and matter M: 

SM C — »■ {LGT[G, M}}, (6.1) 

where Center [G] = C. 

2. Non-existence: Another source of complication has to do with non-existence 
of such a map in the continuum limit. In fact, it is not easy to find continuous 
critical phenomena in gauge theories. As an example, among the pure Yang- 
Mills theories with gauge group SU(N), only in the case of A^ = 2, and possibly 
the case of N = oo (see the discussion at the end of fl2|)) exhibit continuous 
confinement- deconfinement transition. All the rest is believed to have first 
order transitions. 

Therefore, given the map ( B.l| ), the critical phenomena in the continuum limit 
of RHS may be non-existent. Define the subspace (G*, M*) of the gauge groups 
and adjoint matter (G,M) that appear in ( |6.1| ), such that in the continuum 
limit LGT[G*, M*] — > GT[G*, M*] the criticality prevails. Then, we can extend 
the map (|6.1|) to the continuum limit: 



SM C —> {LGT[G, M]} — ► {GT[G\ M*]}, (6.2) 

where again C enter [G*] = C. 



We conclude that, the map ( |6.2|) may or may not exists, and even if exists, it 
may not be unique. We observe, however that the non- uniqueness is a positive fact, 
in the sense that, it provides us with a greater space of gauge theories to scan in 
search for continuous critical phenomena. Indeed, it may be possible to find critical 
phenomena either by changing the gauge group G (while keeping the center C intact) 
or by changing the (adjoint) matter content M. 

Now, the last step of the procedure is to employ the gauge-gravity correspon- 
dence to map the RHS of (|6.2|) onto a gravitational background GR[G, M\. Suppose 
that the gauge-gravity correspondence holds for the subspace (G**, M**) of the pairs 
(G*, M*) that appear in (|6.2|) . Then, we can extend the map as, 

SM C — > {LGT[G, M]} —> {GT[G*, M*}} — » {GR[G**, AT*]}, (6.3) 

where again Center[G**] = C. The last map is the highly non-trivial gauge-gravity 
correspondence that is assumed to hold for arbitrary N, (not necessarily in the large 
N limit). By the standard lore of the gauge-gravity correspondence, the center sym- 
metry C should correspond to a bulk gauge symmetry C on the gravity side. For 
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continuous C, this can be a continuous isometry of the gravitational background. In 
the case C is discrete it may be a continuous isometry broken down to a discrete 
subgroup C by stringy effects. 
Some comments are in order: 

1. Top-bottom approach: From the above procedure it is clear that one ar- 
rives at an operational definition of the gravity-spin model correspondence. 
Take a spin-model SMq that exhibits continuous criticality. Then one should 
scan through all of the LGTs {LGT[G,M]} with various adjoint matter M 
and gauge group G, such that the critical phenomena persists in the contin- 
uum limit. This step, in principle can be done with Monte-Carlo simulation 
techniques. The outcome of this step would be the space of gauge theories 
{GT[G* , M*]} in ( K.3| ). The next step is to construct D-brane configurations 
that correspond to these gauge theories. In general this would only be possible 
for a subspace of theories (G**, M**). The next step then, is to take the decou- 
pling limit of the D-brane configurations to find the gravitational backgrounds 
{GR[G**, M**}} that appear at the end of (|fD|) . The final step is to look for the 
black-hole solutions and study the Hawking-Page transition at finite tempera- 
ture. One can then compute observables of the spin model around criticality, 
such as the scaling of the correlation functions of order parameters, critical 
exponents, the transition temperature T c etc. by holographic techniques. This 
operational definition of the duality corresponds to the so-called top-bottom 
approach in the gauge-gravity duality, that is unfortunately unpractical. 

2. Bottom-up approach: Instead, one may adopt a "phenomenological" ap- 
proach and search for continuous critical phenomena directly on the gravity 
side. This is the approach that we take in this paper. The symmetries, the 
(bulk) matter content and various dynamical phenomena (such as spontaneous 
symmetry breaking) on the gravity side should then hint at what kind of spin 
model that the gravity theory describes. Of course, one should check by compu- 
tations on the gravity side that the theory indeed fulfill the basic expectations 
of the spin model. As a last comment the bottom-up approach is not neces- 
sarily doomed by lack of predictive power. As we argued in section 4, there 
exists a notion of universality exactly around the transition region for models 
that are based on Einstein-dilaton gravity. 

3. The "large-N" limit: The meaning of the "number of colors" and the large- 
N limit becomes clear in this approach. On the gravity side, it corresponds 
to the small Gn limit where one arrives at a classical string theory in which 
the interactions between the bulk fields can be ignored. On the spin-model 
side in the case C = Z^ for N > 4 the corresponding gauge group can only 
be SU(N). Then the limit iV — > oo corresponds to the limit where the Z^ 
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invariant spin-model becomes [/(l)-invariant. In this case — "number of colors" 
correspond to the number of spin states that the spin vector s on a lattice site 
can attain. Of course the large-color limit exists only in theories with gauge 
groups SU(N), Sp(N), SO(N) (with possible additions of U(l) factors). In 
the other cases, one cannot study the spin-model in a large-N approximation. 

4. Discrete C: One may then think that it is never possible to study a spin 
model with discrete symmetry, by a gravitational theory in the small Gn limit. 
This is not necessarily the case. As an illustration, consider the Ising model 
in d-dimensions. The symmetry group is Z 2 and one of the gauge groups that 
has this, as a center symmetry is Sp(N) with arbitrary N. Thus the center 
remains Z 2 also in the large-N limit! Therefore one can make the following 
proposal for a gravity dual of the Ising model: Consider a D-brane set up 
that is dual to YM theory with Sp(N) gauge group with additional adjoint 
matter content M chosen such that the theory exhibits continuous critical 
phenomena. Then consider the background that is dual to this configuration. 
Then the black-hole solution near the Hawking-Page transition should fall into 
the same universality class as the Ising model. By this procedure, it may then 
be possible to analytically calculate the critical exponents of an Ising model in 
any dimension d. 

7. Discussion 

7.1 Summary 

This paper has two related purposes. The first one is to advocate a particular 
approach to holography in condensed matter systems. We proposed to establish the 
link between certain spin systems and gravity through the better understood case 
of gauge-gravity correspondence and the IR equivalence between gauge theories and 
spin-models. The latter is expected to hold only around criticality (in the continuum 
limit). Therefore a gravity-spin model duality is expected to hold only near the phase 
transition region. In particular one should not rely on the gravitational description 
in the UV of the spin-system. 

On the other hand, precisely around the critical region, where the spin-system is 
strongly correlated, the dual gravity description is expected to simplify as the higher- 
derivative corrections become smaller. We showed that this expectation indeed holds 
in a specific gravity model based on non-critical string theory. 

This example also hints at a kind of universality in the dual gravity theory which 
only arises in the transition region: we found that regardless of the details of the 
gravity theory, the physics around criticality is governed by a linear-dilaton CFT. 
Moreover focusing on the lowest states in the CFT at levels N = and N — 1 imply 
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mean-field scaling — in the semi-classical approximation where only the lowest lying 

states are kept in string propagation — regardless of the matter content of the CFT. 

We emphasize that the approach formulated in the previous section in principle 

allows for a top-bottom constructions in AdS/CMT. What more can be learned from 



gravity around the critical region is described in section pA . 

A second purpose of the paper was to construct a model of holographic super- 
fluidity based on continuous Hawking-Page transitions in gravity. A duality between 
gravity and spin-models of the type described above provides motivation for this 
model but it could have been constructed with no reference to such arguments. In- 
deed all one needs from the phenomenological perspective is a gravitational model 
1) with some mechanism of spontaneous breaking of global U(l) symmetry in a con- 
tinuous transition and 2) a bulk field that is charged under this U(l) which would 
serve as a dual of the order parameter. In the model that we studied here the U(l) 
is the topological shift symmetry of the NS-NS two-form field that breaks down at a 
continuous Hawking-Page transition and the fields that are charged under this U(l) 
are the winding modes of the string around the time-circle. Viewed from this per- 
spective, one wonders if a gravity model can be obtained in a more direct fashion 
by truncating the string down to the bulk dynamics of gravity, dilaton an Abelian 
gauge field and the winding modes. In this approach one expects to study an effective 
action of the sort, 

S~ J e - 2 * (R+-(d<5>) 2 + V{$) + -\DT\ 2 --m 2 T \T\ 2 \ (7.1) 

where T (T*) is the winding tachyon with w = +1 (w = —1) and mass — m\ and 
it is minimally coupled to a U(l) gauge field through D = d — iA. The gauge 
field may arise either from reduction of the 5-field on the time-circle or gauging the 
aforementioned topological shift symmetry of the B-field. In this effective theory, 
the Goldstone mode in the superfluid phase is given by the phase of the T field. The 
system would be in the ordered phase when (the particular mode that corresponds to 
the order parameter in the fluctuations of) T becomes normalizable above a certain 
T 

Such models have the same flavor as the ones in 0, || and more recently in [f42 



One immediate future direction is to understand holographic implications of a model 
such as ( \l. 1|) . 

Another immediate future work concerns going beyond the mean-field scaling 
at criticality. We showed that the lowest mass sector of the linear-dilaton CFT 
gives rise to mean-field scaling. We named such a restriction to the lowest lying 
modes in the tree-level string path integrals as the "semi-classical" approximation. 
Then, the question is what happens beyond the semi-classical approximation? Can 
one produce exponents beyond mean-field scaling in this manner? Can one obtain 
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universal exponents of the 3D XY-model by summing up contributions of all string 
states? 

Further future directions are listed below. 

7.2 Outlook 

• Embedding in string theory 

Clearly, it is of great interest to look for examples of the proposed correspon- 
dence in a consistent truncation of critical string theory. The very recent papers 
[£3|[<l [| may be relevant for this enterprize 46 . Another relevant work is [IS, W 



where it was shown that linear-dilaton type geometry universally arises from 
non-conformal branes. 

Explicitly put, one should find a consistent truncation of string theory which 
possess small slack-hole solutions that exhibit continuous Hawking-Page transi- 
tions. We observe that the asymptotic form of the scalar potential in ( p.l5| ) and 



( |3.18| ) is sum of exponentials that quite generically appears in consistent trun- 
cations of IIB and IIA critical string theory. We shall leave this investigation 
for future work. 

An explicit D-brane set-up? 

Even if one finds examples of continuous Hawking-Page transitions in string 
theory, this would not necessarily give control over the microscopic condensed 
matter system that we want to describe. On the other hand, the prescription 
proposed in section 6 in principle goes beyond this and allows for a top-bottom 
approach. 

Therefore, one should search for examples of gauge theories with gauge group 
G and adjoint matter such that the theory exhibits criticality at some finite T c . 
There are indeed examples of this. In Hl4| , the authors studied SU(N) with 
adjoint matter on S 3 , in the large N limit and showed that for certain choices 
of the matter, the theory exhibits a second order deconfinement transition at 
finite temperature, at weak coupling. This happens when the coefficient of the 
quartic term in the effective action for the Polyakov loop is negative. Whether 
this transition prevails in the limit when the radius of the sphere becomes large 
(the case relevant here), or whether it is continuously connected to a transition 
at strong coupling is unclear, but it is probable. 

Discrete center 

We note that, the proposal advocated in this paper can also be applied to 
spin-models with discrete symmetry groups in principle. In most of the paper 



Hi 



We thank Yaron Oz for mentioning possible relevance of these works. 
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we focused on the SU(N) LGT in the large N limit. Going beyond the large 
N limit seems to be a difficult enterprize at the moment, however one may 
consider other gauge groups such as SO(N) and Sp(N) with adjoint matter, 
in the large N limit. The latter is particularly interesting, because it has the 
center Z2 for arbitrary N, hence also in the large N limit. It is very tempting 
to employ the ideas developed in this paper to this particular case to arrive at 
a gravitational description of the 3D Ising model. 

Other critical exponents 

One can also study critical behavior in other quantities. One such quantity is 
the susceptibility: 

x = ii m ^Mk „ t -7 as t^O. (7.2) 

h->o dh 

Here h denotes and external magnetic field and the expectation value M = 
(s) is taken in the ensemble with an external magnetic field present, i.e. the 
Hamiltonian of the spin model replaced by H — > H + h ■ Y2i ■?«• 

How to generate such an external magnetic field in the gravity picture? One 
hint is that the external magnetic field should break the U(1)b invariance 
explicitly by analogy with the dual spin-model. This may happen through a 
Chern-Simons type coupling in the gravity action J d d+l xB2/\Hd-i where B2 is 
the NS-NS two form and H is an appropriate RR form. The role of the magnetic 
field would be played by a constant if -form on the d-1 dimensional space part. 
It would be very interesting to investigate this issue in the future and eventually 
compute the critical exponent 7 in (|7.2j ) by gravitational methods. 

The UV geometry 

We observed that most of the interesting scaling behavior in the observables 
of the spin-model depend on the IR geometry on the gravity side. We did 
not have to specify the UV geometry so far, but we tacitly assumed that it 
becomes asymptotically AdS, for consistency in holographic applications. From 
a practical point of view, the UV geometry will be important if one desires to 



obtain the full form of the n-point functions, not just the scaling with t. In [15| 
we indeed constructed analytic kink solutions that fulfill this promise. These 
solutions interpolate between an asymptotically AdS geometry (with constant 
dilaton) in the UV towards an asymptotically linear-dilaton geometry in the 
IR. It will be very interesting to study correlation functions holographically 
obtained from these backgrounds. 

The specification of the IR geometry follows from physical requirements of 
the spin-model near the transition region T c . On the other hand, the black- 
hole with temperature T is argued to correspond to the super-fluid phase with 
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temperature 1/T. Then, one can ask whether we can also produce the expected 
behavior of the super-fluids at very low temperatures, by specifying the high 
T regime of the black-hole that corresponds to the UV geometry: One basic 
feature of the two-fluid model for super-fluidity is that the (normal) speed 
of sound, that is associated with fluctuations in the magnetization vanish as 
T — > 0. This is certainly not a behavior expected from an asymptotically 
AdS geometry in the UV which would correspond to a conformal fluid with 
c 2 s = 1/3. We conclude that the kink solution that flows from AdS to linear- 
dilaton |15| would not do the job here. One possible way to proceed may be to 
consider the non-conformal brane solutions [[15|, [|(J which on one hand allow for 
a holographic computation of observables, and on the other, there is a chance 
to find backgrounds with c 2 s — > in the UV. 

Do we really expect to find a background as a solution to two-derivative 
Einstein-dilaton theory, that would produce the desired behavior in the entire 
range T e (0, T c ) of the super-fluid? The answer is most probably negative. 
Let us suppose for a moment that such a background exists as a solution to 



the full d + 1 dimensional non-critical string theory. As we showed in |L5] , 
the curvature invariants in the IR vanish in the string frame, therefore a two- 
derivative approximation is expected to work in near the transition T rj T c On 
the other hand, the invariants away from the transition region are determined 
by the intrinsic string scale £ s . This means that in a two- derivative approxi- 
mation one deals with a background with £/£ s ~ 1. To conclude: we indeed 
expect non-trivial a' corrections in the UV region and the two- derivative ap- 
proximation presented here is expected to give reliable results only near the 
transition region. 

The two fluid model of super-fluidity: 

We only performed the computation of the speed of sound for the Goldstone 
mode. It would be also very interesting to look at the dissipation coefficients 
associated with these fluctuations. In the two-fluid model of super-fluidity, 
one deals with a coupled system of pressure and entropy waves of the two- 
component superfluid, cf. |47| for a recent review. The pressure waves are dual 



to metric and dilaton fluctuations, whereas the entropy waves are associated 
with fluctuations of the B-field. It would be very interesting to work out this 
coupled system of fluctuations in order to determine the associated dissipative 
fluid dynamics 47 . 

Spin models with non-Abelian symmetry groups 



47 It seems that one needs to turn on a Chern-Simons term of the form J B A i?_R_R where Hrr is 
a d-l-form turned on the spatial directions in order to achieve such a mixing. 
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One fundamental restriction of the approach in section 6 is that the spin sym- 
metry cannot be non-Abelian as it follows from the center symmetry of the 
corresponding lattice gauge theory. On the other hand a very important model 
for superconductors involves the 0(3) model 48 . Whether one can overcome 
this restriction in our set-up is an interesting question. In phenomenological 
models such as H8J , one can achieve this simply by considering black- holes with 



non-abelian charges. 

In our perspective, one idea is to consider the enhanced symmetries of string 
theory at special radii fl6f . When the sting is compactified on the time-circle 
one obtains U(1)g x U(1)b symmetry at an arbitrary radius. The second one is 
spontaneously broken in the BH phase. At special a radius T = T s = (2n£ s )~ 1 
(in bosonic NCST) however one obtains an enhanced symmetry SU(2)l x 
SU(2) R . If this radius corresponds to the transition temperature T c then one 
may be able to obtain a model with the desired behaviour within our set-up. 
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APPENDIX 

A. Simplest example of the LGT-spin equivalence 

Let us review how the LGT-spin equivalence works in the simplest case of the U{\) 
lattice gauge theory in d dimensions. Through this example we will illustrate that 
the temperature of the spin model is inversely related to the temperature of the 
(lattice) gauge theory which also holds in the most general case. 
In the Hamiltonian formalism, the lattice theory is defined by: 

H = f a J2 E2 ( r ^ -^J2( v F} + v^]) ■ (a.i) 

(r,h) T 

Here, g is the coupling constant, a is the lattice spacing. The first sum is over the 
links (r, n) on a d-dimensional square lattice (r denotes the lattice site, ft denotes the 
direction of the link that originates from this site) and E denote the electric fields 
residing on these links. The first sum above yields the electric energy. The second 
one is over the elementary plaquets I\ The Vs denote the Wilson lines on these 
plaquets. This gives the magnetic energy. The partition function of gauge invariant 
states at temperature T} is given by, 

Z lat (T l ) = Tr'e- H / T \ (A.2) 

where the prime reminds us that we have to impose the Gauss' law on the states in 
the ensemble. Consequently, the sum above is over the gauge invariant states \ip) 
which should satisfy, 

HW) = E\*I)), r,#) = 0. (A.3) 

Here, the second equation imposes the Gauss' condition on the states; the operator 
T r is the lattice analog of V • E on each lattice site r : 

r r = ^£(r,n). (A.4) 

In the strong coupling limit, one can drop the magnetic energy term in ( |A.1| ). 

Now, the sum is only over the electric link variables and the prime can be removed 
by a suitable Lagrange multiplier a: (at strong-coupling) : []TD|, |TT 



ziat{Ti)= j n^)nE ex p(-^ £ ( f ' A ) 2+i [ a ( r )- a ( f+fi )i E ( r '4 

•'-' K r links E \ a l ' 

(A.5) 
The integral over a imposes the Gauss' law. Using the Poisson summation formula, 

Y^e cE2+iaE oc J2e-^ a+2nm) \ (A.6) 
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the sum over the E can be performed: 



Z lat (T m ) = r l[da(r) J] ]Texp (~[a{r) - a(r + h) + 27nn r A . (A.7) 

J~ n r links m r ^ ^ ' 

This is the Villain approximation to the Heisenberg model for ferromagnetism in d- 



dimensions. It is in the same universality class with the Heisenberg model JO]]. In 
particular, for d = 2 this becomes the famous XY model in 2 dimensions, where the 
BKT scaling was first observed [^ . 



The salient feature of the Heisenberg model is that it exhibits order in the low 
T phase (that corresponds to the, high T deconfined phase of the Abelian LGT), 

(e 4a(i V Q(0 ) oc 1, as, R^oo, (A.8) 

and disorder in the high T phase (that corresponds to the low T, confined phase of 
the LGT), 

(e ia ^e iai °) oc ET R/ *, as R^oo, (A.9) 

where £ defines the correlation length. Some comments are in order: 

• The computation can be generalized to the non- Abelian case [[IT], [TT| . 

• The computation is performed in the strong coupling limit where one can ignore 
the magnetic energy in ( A.l ). This constraint can easily be loosened and the 



equivalence prevails also if one considers the magnetic piece [|ITJ, [TT| . 

One can generalize to add adjoint matter, as the center symmetry remains 
intact under addition of adjoints. 



B. Relation between non-critical strings and the linear-dilaton 
theory 



It is long known that the two theories are intimately connected |31| . The connection 
is made precise in a beautiful work by Chamseddine [^jj which we would like to 
review here. 

For simplicity we consider bosonic matter. Then the non-crirtical string theory 
in d — 1 spatial dimensions with flat Euclidean target-space metric can be defined by 
the world-sheet action 

A ws = — [ d 2 o-Vh<p(R+A)+\x(M)+n f d 2 o-Vh+— [ d 2 aVhh ab d a X l d b X\ 

(B.l) 
where we ignored coupling to S-field for simplicity. Here is a scalar field introduced 
in [OUJ in order to ameliorate evaluation of higher genera diagrams. Its presence was 
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also motivated on physical grounds ||29|| . The matter index runs from i — 1 to 
i = d — 1. The manifold M can be with arbitrary genus x{M) = 2(1 — g), and the 
constant A is determined by the asymptotic value of the dilaton. The constant A is 
a free parameter and \i and is subject to renormalization. It was argued in [B!J that 
the coupling in the action apparently overcomes the difficulties, "the c=l problem" 
encountered in the study of non-critical strings. The path integral over the world- 
sheet metric can be performed in the conformal gauge h ab = e aL h ab and results in 
the Liouville action for the field <7l producing additional world-sheet terms [ic aL and 
Ae CTL . Then the effective renormalized action after gauge-fixing involves the matter 
part as in (|B.1|) , the renormalized gravity action, the ghost part that arise from the 
reparametrization fixing and the induced Liouville action: 



A ws = -*- / d 2 aVh ( j ) (R + A h a L + Ae 2f7L ) + -*- / d 2 aVhh ab d a X i d b X l 
2vr J M n Ana' J M 



+ Xx{M) 



-*- / d 2 aVhb ab V a c b + - I d 2 aVh 

27r hi 7T J M 



a{-h ab d a a L d b a L + a L R) + /ie 2<Ti 



(B.2) 



Here a is a constant that vanishes only in the critical case (e.g. a = (25 — <i)/12 for 
the bosonic case). The world-sheet terms fie (7L and Ae°" L can be shown to correspond 
to conformal primaries, hence one can treat them as marginal deformations of the 
theory with fi = A = 0. The full theory can be shown to be free of conformal anomaly 
and has a well-defined OPE among the fields X l (z), (fi(z) and <tl( z )- 

The valuable observation of [2{| is to interpret (linear combinations of) $ and 



the Liouville-field <7l as two new additional dimensions of the target-space. The 
resulting theory is described by the new action (for the case /i = A = 0), 



A ws = - [ d 2 aVlv,X fl R+-^— [ d 2 aVhh ab d a X f *d b X u r ] „ u +— [ d 2 a\ r hb ab V 1 d 
7T Jm 47ra' J M 2vr J M 

(B.3) 
where fi runs from to d and rj^ is the flat Minkowski space metric. The additional 
dimensions Xq and X d is given in terms of ctl and (ft of the non-critical string as, 



,.„ 6a' , .. /25 — d / 6a' , ,„ , 

X ° = \ ; d>, X d = 2\ <Jl + \ -A B.4 

A/ 25-d^' V 6a' L V25-CT' V ; 



and the coefficient v M in (p.3|) satisfies the condition 



25 -d , . „ 9-d 

bosonic; v^v a = 

6a' ** 4a' 



v^V/j, = bosonic; v^v^ = fermionic, (B.5) 



where we also show the condition in the fermionic case for reference. This is nothing 
else but the linear-dilaton theory that arises in the IR limit of our geometry. There- 
fore the non-critical string theory in d — 1 spatial dimensions is equivalent to a d + 1 
dimensional linear-dilaton theory. 
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The main advantage of mapping the linear-dilaton theory to the non-critical 
string theory is that the latter provides a well-defined CFT. The spectrum as well as 
the arbitrary genus path integrals can be evaluated |2j| . Another great advantage is 



that one can generalize this construction to include fermions on the world-sheet and 
N = 1 world-sheet super-symmetry This opens the Pandora's box and a rich variety 
of linear-dilaton theories can be constructed with various possible GSO projections, 
twisted or shifted boundaries for the bosons with various combinations of NS or R 
fermions including the heterotic case. 

C. Some background in statistical mechanics 

In this section we review some standard background in statistical mechanics that we 
need in the following section. 

We take the XY model as our example although the approach can be very general. 
Consider the spin-model that is described by the Hamiltonian, 

H = -J^2si-8j + --- (C.l) 

(ij) 

where J is the interaction strength which is positive for a ferromagnet and negative 
for an anti-ferromagnet, s is the spin vector that rotates on a plane and (ij) denotes 
nearest neighbor pairs on the lattice. The dimensionality of the system only shows 
up in the number of nearest neighbors of a lattice cite. 

C.l Landau action 

The Landau action can be defined from the partition function in a formal way, with no 
resort to any approximation, see for example [ I5| . One first introduces a continuous 
spin density, 

s(x) = y^5(x — Xj)si. (C.2) 

i 

Than one defines the local magnetization as in the following formal identity 

1 = / Vm5[m{x) - s{x)}, (C.3) 

inserts the RHS of this identity in the partition function 



Tr 



(C.4) 



Vm8[m(x) — s(x)]e ? 
and finally performs the sum over the spin degrees of freedom Si. This yields, 

Z L = f Vme-P FL[rh] (C.5) 



that defines the Landau action F^. This is formal but exact. Of course, the last step 
of summing over the spin degrees of freedom is practically impossible in most of the 
spin systems. 
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C.2 Landau approximation 

This corresponds to keeping the most dominant contribution ttil in ( |C.5| ): 

Z L ^e~ pFLi - mL) . (C.6) 

The most dominant contribution itll is determined by minimizing Fl- Furthermore, 
near an order-disorder phase transition the total magnetization 

M= — *— [ d^x m(x), (C.7) 

Vd-l J 

should go to zero. Then the 0(2) symmetry of the XY model dictates the following 
general form of the Landau action: 

d d - x x (a (T)\dm(x)\ 2 + ai (T)\m(x)\ 2 + ^a 2 (T)\m(x)\ 4 + ■ ■ ■ J (C.8) 

where the ellipsis stand for higher (even) powers in m. 

Basic quantities to be computed, which determine the phase diagram of the spin 
system involve the functions Qo(T), cti(T),... 49 For example the point a.\{T c ) = 
0, OJ2(T c ) 7^ corresponds to a second order transition. The point ot\{T c ) = a2(T c ) = 
corresponds to a tri-critical point, etc. A further simplification occurs in the 
case of positive J (ferromagnet), when the ground state of the system should avoid 
fluctuations in local magnetization rh(x) because they increase the energy of the 
system, hence one can ignore the kinetic term in (|C.8| ) and one can set (m(x)) = M. 
This means that the ground state has isotropic magnetization. 

C.3 Mean- field approximation 

Mean- field approximation is a standard method to compute the Landau action ( |C.8|) . 
One expands around a mean-field Si = M + 5si where Si denotes fluctuations around 
the mean value M. One substitutes this in the Hamiltonian ( fj.l| ). Ignoring the terms 
of second and higher order in 5s corresponds to the mean-field approximation. One 
can clearly compute the partition function, hence the Landau action Fl analytically 
within this approximation. One immediately obtains, 

Z = Tr e-# = e ~* = e^^ 2 J] f dO^^ ^ (C.9) 

where N is the total number of sites on the lattice and z is the number of nearest 
neighbors. Evaluating the integrals one obtains the following Landau action: 



f m* = _jvt log(27r) + NzJ\M\ 2 - iVlog 



2zJ\M\ 

•Ai 1 jT~ 



(C.10) 



49 In a more complicated system, for example with a chemical potential, these functions depend 
on additional variables. 
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The first term corresponds to the entropy and can be ignored for our purposes here, 
as it is identical on both phases. The rest corresponds to the energy of the system. 
Expanding energy near T c where M is small, one obtains, 

Fr = \M\ 2 Nj{z- Z ^+... (C.ll) 

At the second-order transition the mass term should vanish. Thus one obtains the 
mean-field value of the transition temperature: 

T c MF = zJ. (C.12) 

For a square lattice in d — 1 spatial dimensions z = 2(d — 1). 

In the context of Landau approximation 50 , the further mean- field approximation 
means that the Landau coefficients admit a Taylor series expansion near T c . For 
example 

cci,mf(T) = a c (T - T c ) + ■ ■ ■ (C.13) 

and the basic data to determine involve the quantities T c and a c in this case. This 
linear dependence is explicit in ( |C.11| ). 

Clearly, the mean-field approximation is crude and one can compute the afore- 
mentioned observables to a greater accuracy by the renormalization group methods 
or Monte-Carlo simulations. We shall see below, how the gravitational techniques 
can go beyond the mean-field approximation. 

C.4 Gaussian fluctuations 

Fluctuations around the mean value M yield crucial information on the spin-model, 
in particular the spin-spin correlation function and the associated critical exponents 
near T c . To compute them, one substitutes 51 

m(x) = M + s(x), (C.14) 

in ( |C.8| ) and expand to second order, ignoring higher order terms in the mean- field 
approximation. 

One introduces the correlation length as the natural length scale, 



« T) = fe (ai5) 

in (|C.8|) 52 . The calculation of the spin-spin correlation function in the ordered phase 
(T < T c ) within this approximation is standard, see for example [[±9| . 



50 What we mean by the Landau approximation is summarized by eqs. (C.6) and (C.8) without 
further specification of the Landau coefficients «i(T), 012 (T) etc. 

51 We drop the S in front of s from here on, for notational convenience. 

52 Strictly speaking this corresponds to the correlation length for T < T c . In the low-T regime it 
differs by a factor of \Sqrt2. However, we are mainly interested in the scaling of £ near T c and the 
scaling is the same in the mean-field approximation, from below and above. 
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The only crucial point is that, in the case of spontaneous symmetry breaking, as 
in the XY-model, one has to decompose the correlation function in parts longitudinal 
and transverse to the direction of magnetization M. Let us denote the components 
of m by mi. Introducing the unit vector along the direction of magnetization in the 
system, 

Mi 



\M\ 
Then one finds, 



(C.16) 



T e" L/?(T) T 1 

(rm(x) mj (0)) = \M\ 2 v iVj + 4?rT Ld s vw + ——fa - v iVj ), (C.17) 

where the result only depends on the radial distance L = \x\ by rotational symme- 
try. Here the first term comes from the disconnected part of the correlator, and it 
is present only in the ordered phase. The second and third terms are the pieces lon- 
gitudinal and transverse to the magnetization respectively. The longitudinal piece 
arise from massive fluctuations in the Mexican hat potential where the typical mass 
of the fluctuations is given by mi = £~ l . This attenuation term is missing in the 
transverse correlator because the fluctuations correspond to the massless Goldstone 
mode. 

The result ( p.l7| ) is valid in the mean-field approximation where we only treat 
Gaussian fluctuations. This approximation will break down if the system is strongly 
correlated. Generally in condensed matter systems, strong correlations arise around 
a phase transition. Therefore we expect that the gravity dual becomes a good de- 
scription in the transition region. The notion of strong fluctuations is quantified by 
the Ginzburg criterion: 

£ 5 " d «^, (C.18) 

a 2 T c 

where «2 and «o are the Landau coefficients in ( |U.8| ). For fluctuations of s we see that 
the mean-field approximation unavoidably breaks down near the transition where £ 
becomes very large 53 . 

Beyond the mean-field approximation, one has to take into account non-trivial 
self-energy corrections to the correlator that generically result in the anomalous ex- 
ponent rj. Therefore the generic form of the correlator is similar to ( |U.17| ) but with 



the additional anomalous dimension, in addition to the engineering dimension in the 



53 One can be more careful by considering the amplitude and phase fluctuations of s separately. 
In the former case the coefficients ao and oli stay constant at T c and from ( |C.18 ) one finds that 



strong correlations are indeed unavoidable in the transition region where £ diverges. In the latter 
case also the constant «o vanishes near the transition, see below. In the mean-field scaling ao ~ t 
and £ ~ £~2. Therefore one finds that, only for uninteresting dimensions d > 9 the mean-field 
approximation is expected to be good for the phase fluctuations. 
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correlator: 



T e - L/ ^ T) T 



(m t (x) m 3 {0)) = \M\ v iVj + — L ^ 3+v v iVj + — ^^-((Jy - Vl v 3 ), (C.19) 

The mean-field approximation corresponds to n = 0. 

Equation ( |(J.19| ) gives the correlator in the ordered phase that is dual to the 
black-hole solution in gravity. Above the transition (M) vanishes and there is no 
Goldstone mode. Hence, the correlator is given by the second piece of (|C.19|) . This 
is dual to the thermal graviton gas phase of gravity. 

In a second order phase transition, the correlation length diverges as T — > T c , as 

£(*) ~ *-", (C.20) 

where v defines a critical exponent. In the mean-field approximation v = 1/2. 

Another important point concerns the scaling of M near the transition. It van- 
ishes in a continuous transition, as it is the order parameter. As explained above one 
can ignore the kinetic term in (|C8|) in the ground state of a ferromagnetic system. 
Vanishing of M near T c is characterized by the critical exponent /3 54 : 

M~t p . (C.21) 

In the Landau theory the expectation value M is determined from (|C8|) as, 



(m(x)) = M=^/i^l. (C.22) 



As a 2 and a stays constant at T c , we see from ( |C.15|) that, in the mean-field ap 



proximation the scaling of M and £ are inversely related, 

M oc r 1 oc t u . (C.23) 

Therefore, it suffices to determine how M scales in order to obtain the scaling of £ 
in the mean-field approximation. 

Futhermore, in the mean-field theory the coefficient a(T) in ( p.8|) is assumed to 
be linear in t near T c , ( |C.13|) . Comparison with ( G.22Q and ( |C.23|) then shows that 
in this approximation: 

Pmf = vmf = -■ (C.24) 



54 Not to be confused with the perimeter of the time-circle j3 = 1/T. We use the same notation 
for these quantities and which one is meant should be clear from the context. 



-78- 



C.5 Vanishing of the second sound 

We note that, the vanishing of M as T — > T c from below implies vanishing of the 
sound velocity associated with the phase fluctuations. To see this we represent the 
fluctuations s(x) in ( p.l4| ) as, 



m{x) -> (\M\ + p(xj) e ^-H^(*). (C.25) 

Substituting this in (|C.8|) , one obtains the kinetic term for the phase fluctuations, 



5F L ~ / d d ~ l x (-?(T)\M\ 2 (5^(x)) 2 + • • -J (C.26) 

Therefore, the speed of sound associated with the phase fluctuations vanish near T c , 

c^ ~ t 2/3 , t ->■ 0, (C.27) 

and the rate it vanishes is determined by the critical exponent associated with the 
magnetization ( |C.21| ). 

In the derivation above, we used the Landau approximation and only kept the 
leading terms in fluctuations. Therefore, within this picture the magnetization crit- 
ical exponent should be the mean- field one, ( |U.24j ). This means that in this picture 
one obtains, 

Cprvt, t->0. (C.28) 

An important check for the proposed gravity-spin model correspondence here is to 
derive the same scaling law on the gravity side. This is done in section |5.6| . 

C.6 BKT theory 

Finally, we consider the XY-model in two spatial dimensions. As well-known, in 
less than three dimensions, long-range order is destroyed by the IR divergences in 
fluctuations of the order parameter [|50|] , i.e. there are no Goldstone bosons. However 
Berezinskii, Kosterlitz and Thouless |SJ observed that the 2D XY-model still serves 
as a good model for superfluidity. The main observation is that, although there is 
no long-range order in the standard sense, there exists a topological order below a 
certain T c , where the vortex- anti vortex pairs condense. Above T c the system is the 
"deconfined" phase where the vortex anti-vortex pairs are liberated and one has a 
plasma of vortices. All of this is of course very similar to what happens in QCD, 
with the replacement of quarks with "magnetic" quarks. 

Vortices are charged objects. One assigns vortex charge Q v = ±1 for the vor- 
tices and anti- vortices respectively. The total vortex charge in a configuration should 
vanish in two-dimensions because the gauge field has an IR divergence and the en- 
ergy of an unbalanced configuration would diverge, hence its Boltzman factor in the 
ensemble vanishes. 
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What do we expect from the behavior of vortex correlation functions above and 
below T c ? Let us denote the operator that creates a vortex, localized at point x by 
v(x) and the operator that creates an anti- vortex by v(x). For the reason described 
above, one cannot have any non-trivial expectation value neither below nor above 
the transition, as it would break the vortex charge 55 : 

(v(x)) TG = {v(x)) TG = (v(x)) BH = (v(x)) BH = 0. (C.29) 

The phase of the system can be probed by the two-point function of the vortex-anti- 
vortex pair, however. One finds that the two-point function is exponential in the 
high T phase, hence there exists a correlation length, whereas it has power-law in 
the low T phase: 

(v{x)v{0)) TG ~e- LmT) ; (v(x)v(0))bh~ L p{T) . (C.30) 

where L — \x\ ^> 1 and p is some power. 

In systems with more than two spatial dimensions, one can still consider vortex 
configurations, however they would not have as significant effects on the phase of the 
system as in 2D. The relative objects in higher dimensions would be the vortex-lines, 
planes etc, that are analogous to monopole configurations in gauge theories. 

D. Fundamental string action 

Here we fill in the details of the computations in section [|. The F-string action 
involves two terms 56 : 

Sng = Sg + S<p (D.l) 

where 

Sg = i f f dadT ^ habd ^d h X"G, U) (D.2) 

S* = — [ i dadr y/det h ab # (2) $(X(a)), (D.3) 

271 J « 

where l s is the string length, r/ is some turning point of the string embedding that 
will be specified in the following, R^ is the Ricci scalar that corresponds to the 
world-sheet metric h a b, and G^ u is the BH metric in the string frame: 

ds 2 = e 2A s{ r) (/-l( r ) dr 2 + dx 2_ i + ^f^ ? A ^ r ) = ^( r ) + ^T^ 7 ")' ( D ' 4 ) 

The on-shell value of the action depends on the boundary conditions of the string. 
In this paper we consider three separate cases: 



55 The TG (BH) phase of gravity is dual to the high (low) T phase of the XY-model, whence we 
denote the vortex correlators accordingly. 



o(i 



The role of a non-trivial B-field is already discussed in the test and we shall ignore it here. 



-80- 



1. The Polyakov loop, 

2. The Polyakov correlator (The Wilson loop), 

3. The 't Hooft loop. 

D.l The Polyakov loop 

First of all, we fix the world-sheet diffeo-moprhism invariance as 

a = xq, r = r, (D.5) 

where Xq is the Euclidean time. 

We consider two separate geometries that the string is embedded: a) the thermal 
gas solution with topology S 1 x Bj where Bj is a d-dimensional ball and b) the black- 
hole with topology D 2 x Bd-i where D 2 is a 2 dimensional disk. 

The boundary of Bd is S^ 1 but we are interested in the flat limit where S d ~ 1 — > 
R d ~ l . Therefore, in both cases the boundary of space-time becomes S l x i? d_1 . The 
string ends on a curve C on the boundary where C = S 1 x P, P being a point x on 
R d ~ l that we can take as the origin with no loss of generality. In case a, the only 
string solution with C as the boundary is the semi-infinite cylinder S 1 x R where R 
is isomorphic to the radial coordinate r and P is the point that corresponds to the 
endpoint of the line R at r = 0. In case b, the only string solution that ends on C 
is isomorphic to D 2 , hence it wraps the entire D 2 part of the bulk geometry. 

Clearly, the action ( |L>.2| ) diverges in case a because 07 = r/ = 00. Therefore 
( p.l| ) will diverge unless there is some cancellation between ( |D.2j ) and ( |L) . 3| ) . In 
the following we show that (|D.3[) is finite in all of the cases under consideration. 



Therefore the result in case a is that the Polyakov-loop vanishes. 

Now, consider the case b, i.e. the black-hole geometry. As explained above, 
the string wraps a D 2 . Then, the radius of D 2 is given by tv Clearly, both (p.2| ) 



and ( p.3|) are finite hence contribute to the energy for an arbitrary but finite r^. 
However, we are interested in the limit T — >■ T c i.e. r^ — > 00 and we ask how do 
( |TX2|) and (|LT3|) scale with r h . 

Let us first consider the scaling of S$ with Vh in the limit r^ — > 00. The world- 
sheet metric h in the gauge ( |L) . 5| ) is given by, 

ds 2 ws = e 2A ° (dx 2 f + dr 2 /f) (D.6) 

One finds, 

VhR^ = -2A' S f - 2 / A" s - /", (D.7) 

where prime denotes d/dr. Using eqs. ( p,16|) and (|D.8| ) one finds that, 

^(r)oc( e r a V ' CaSei : (D.8) 

r , case 11, 
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where k! = kVoo/2 and a are positive constants. Thus we find that, even though $ 
diverges linearly as r — > oo, the integrand in S$ vanishes in this limit. Therefore the 
contribution of S$ is finite in the case b, also in the limit r^ — > oo. 

Now we consider the scaling of Sq with r^. The metric in the string frame is 
given by ( p.4| ): From ( p.2| ) one immediately finds that, 

where e is the UV cut-off that should be removed by a proper counter-term action, 
but it will be irrelevant for the scaling of the Sq with r h . 

In passing we note the trivial result on the TG solution. This corresponds to 
( p.9| ) where rh replaced with oo and A s replaced with A StQ . As A St0 goes to zero 
in the limit r — > oo { \U.9\ j is divergent and the magnetization that corresponds to 
exponential of — S vanishes. 

In the limit rh — > oo in case of the BH, one has Sq oc limr^oo e 2As<rh V/i, whence 
it diverges linearly, whereas S$ remains finite as we showed above. Thus, using the 
fact that A(r h ) -)• (|D.8| ), one finds, 

-s NG „ — ( T ~\ r h) 



P[C] oc e- bNG oc exp ( — ^ r h ) , r h -> oo (D.10) 



The next task is to express r^ in terms of the normalized temperature t ( |3.2U| ). We 



know how $^ can be expressed in t from ( |3.21| ). Thus if suffices to find r^ in $^ in 
large $^ limit. This is given by |3.i6| : 



Now, ( P^TD yields, 



Therefore ( [D.ll| ) gives, 



lim $ h = ^ r h . (D.li; 



, — - \og( t/C), case i, 

*a=h r A-l .: ( d - 12 ) 

(^J , case n, 



/T^« 1 °s( t / C '). case i> 



" ' ' ■ '■ case n. 



Substitution of ( |D.13| ) in ( p,10[ ), use of the fact that T — > T c in the limit r^ — >■ oo 
finally yields 

Case i : e~ SNG oc t^ t -> (D.14) 

Caseii: e" &G oce^^) ", t -)- 0, (D.15) 

In the arguments above, we ignored the issue of renormalizing the UV divergence in 
the action (p.9|). The renormalization can be done by subtracting the self-energy 
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of the single quark that corresponds to a single disconnected string solution that 
hangs from e to r/. This renormalization is considered in detail at the end of App 
D.2\ below. The counter-term action is the same (up to a factor of 2) as there. The 



same conclusion reached there — that ignoring the renormalization does not affect the 
leading term in L in the large L limit — is also valid here. 

D.2 The Polyakov loop correlator 

We compute the on-shell string action that corresponds to the Polyakov-loop cor- 
relator here. We will consider the BH geometry, and the same problem on the TG 
geometry can be obtained from our result below, see below eq. ( p.29|) . 
The fundamental string action is 

Sfi = 2^J dTda ^ [( hah 9% + ^ ab B,v) d a X»d b X» + £ 2 MX)RW] , (D.16) 

where h is the induced metric, g s is the target-space metric of the BH geometry in 
the string-frame 

ds 2 = e 2A ^ {r\r)dr 2 + dx 2 _, + dx 2 J(r)) , A s (r) = A{r) + ^y*(r), 

(D.17) 
and R^ is the world-sheet Ricci scalar. The string that corresponds to the to 
Polyakov-loop correlator (P*(x)P(0)) is a connected string solution with end-points 
x and on the boundary. With no loss of generality, we take these points to lie on the 
same axis that we call X\. The string that connects these points extends towards the 
deep-interior of the d + 1 dimensional target-space in the r-direction. Thus, a good 
choice of the gauge-fixing is given by a = xq, r — X\ where t is the Euclidean time 
coordinate with perimeter 1/T. The string should also wind-around the time-circle. 
As we look for a solution that only depends on r, the r- integral factors out and yields 
a multiplicative factor of 1/T. 

First of all, the 5-coupling cannot arise here because it yields an imaginary 
contribution, whereas the Polyakov-loop correlator is manifestly real. Thus we can 
drop the second term in ( |L) . 1 6| ) . Secondly, one can show that the contribution of the 
dilaton-coupling is sub-dominant with respect to the first term in ( p.!6| ). One can 
see this as follows. Let us assume that indeed the dilaton-coupling is sub-dominant. 
Then, the induced metric is solely determined by the first term in ( |D.16| ) 



h ab = d a x»d b x v g%- (D.18) 

Given this one can compute the Ricci scalar: 



VhR^ = {f + Z 2 ) 



Zf' 2 + 4f 2 (A' s Z' + ZA" S ) + 2Z\2A'J' + /") 



+2f(Z'f + 2Z 3 A': + Z(A'J> + /") 



(D.19) 
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where we defined, 



In this paper, we are interested in how the Polyakov-loop correlator scales near T c 
and for large L. The latter is given by, 

7 r r f 

d Xl = 2 Z-\r)dr, (D.21) 



where e is some cut-off near the boundary 57 . The limit T ^- T c and L large corre- 
sponds to 

r h >r / >f. (D.22) 

In the limit r^ ^> I one can show that the blackness function can be replaced by 
/ « 1 everywhere except r = r^ where it vanishes (see [Q for a derivation, that 



immediately carries over here). At the same time for 77 3> 1 the scale- factor A s (r/) — > 
0. Then, using the sub-leading terms in A(r) (|D.8|) , one finds, 

Z~e- K ^ /4r (Casei); Z ~ r~ a/2 (Case ii), (D.23) 



The constants « and a are defined in ( |3.18| ) and ( p. 19 ). 



In the limit 77 — > 00 the dilaton $(77) diverges as in (|3.16|) , therefore one may 
worry that the last term in ( p,16|) contributes significantly. On the other hand, we 
see from (|D.8|) that the A' s factors (and similarly v4")factors in (p.!9|) are suppressed 
exponentially (case i) and with power-law (case ii) in the region r 3> 1: 

4(r) ~ { C _ a ' X : (D.24) 

r , case 11, 



and the Z factors in ( p.!9| ) are suppressed as ( p.23|) . Using the latter in (|D . 2 1|) we 



find, in the region L ^> 1: 



L~e KVV " /4r f (Casei); L ~ r7 o/2+1 (Case ii), (D.25) 



Thus, we conclude that, in the regime r^ — > 00, and 77 ^> 1, the dilaton contribution 
to the string action (the last term in (p,16|) ) scales as 



. L , case i, 
5$ ~ { 2-I/2 (D.26) 

L 1+Q / 2 , case ii, 



57 We will comment on how to remove the cut-off by appropriate renormalization in the end of 
this Appendix. However, as we are interested in the limit rf — > 00 in this paper, the e dependence 
can be kept, it will not contribute to the results. 
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in the region L ^> 1. 

Below, we show that — with the assumption that one can drop the dilaton con- 
tribution in ( p,16| ) — the string action scales linearly in L: 

S G ocL, L > 1, (D.27) 

where by Sq we denote the on-shell contribution of the first term in (p,16|) with the 
assumption that the last term can be dropped. Thus, we can safely conclude that, 
our assumption in the beginning of this discussion, namely that the dilaton- coupling 
in ( \D.1Q ) does not contribute to the string-solution in the limit ( \D.23j ) is valid in 



case i, and also in case ii, unless a < l. 58 In the case of a = 1 our assumption 
above is violated as both Sq and 5$ and one should solve the full action in ( |U ■ 1 6| ) . 



In the case a < 1 the metric term Sq is sub-dominant to S$ and one can turn the 
aforementioned argument in favor of S$, i.e. one can assume that S$ is the leading 
contribution in the limit L ^> 1. These cases provide interesting examples that 
the ^-coupling becomes crucial in determining the behavior of the Polyakov loop 
correlator (similarly the quark-antiquark potential in holography), however they are 
not of direct interest to us in this work. 

Thus, we can drop the last two terms in ( p,16| ) and the solution for h a b is given 



by ( p,18|) . The on-shell action then becomes a Nambu-Goto action. Given the 



target-space geometry ( p,17| ) finding the on-shell NG action is a standard exercise 
(see for example [pT|): 



1 pr f p 2A s (r) 

S F1 = — — / dr _. (D.28) 

e 4A s (r)j( r ) 



Then the action as a function of L is given by the parametric solution of ( p.21|) and 

(EH- 



We are interested in the limit ( p.22| ) where we can replace / ~ 1 throughout 



the entire range of r, up to r/ and up to the value of r/ slightly smaller than r^ (in 
the limit r^ — > oo, T — > T c . Comparing ( |L).21| ) and ( p.28| ) and using the fact that 



A s — > in the numerator of ( p.28[ ) we conclude, (assuming a > 1 in case ii) 

Sfi -> ^jr L, as L^ ex). (D.29) 

The on-shell action of the same type of connected string solution in the TG geometry 
is given by replacing the metric functions above by / = 1, A s — y A Q S and $ — > $o- 
The same arguments above then directly carry over to this case. 

Finally, let us discuss the renormalization of the action ( |D.28| ): We regulated 
the action by inserting a cut-off at r = e close to the boundary. As one removes 



O.S 



The case of BKT scaling corresponds to a — 2 and it is in the safe region. 
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the regulator, e — > 0, the on-shell action diverges, due to infinite area of the space- 
time metric near the boundary 59 . One should remove the regulator by adding a 
counter-term action designed to cancel the divergences. In the case of Polyakov-loop 
correlator in AdS, this is done by subtracting the self-energie that corresponds to 



two disconnected strings hanging from e to r^ [|3"3"] , |34 



*- = sab/' «"*'*■■ (a30) 



Clearly, subtraction of S ren from Spi in ( |k>.28| ) removes the divergence at e = 0. 



However, as criticized in [j35 |, this makes the counter-term action temperature de- 
pendent. Moreover, ( p. 30 ) is divergent in the limit r^ — > oo, that is the limit that 
we are interested in. Instead, here we propose to regulate the action by subtracting 
the self-energy of two quarks of the Polyakov-loop correlator as in |51| : 



^»~/ f e 2A ^dr. (D.31) 



Clearly this removes the UV divergence in ( |D.28| ) in the limit e — > and yields a 



finite result for finite L. It is also apparent from ( p.31| ) that it diverges in the limit 
L — > oo but this divergence is physical. The only point that we have to worry about 
is that, it does not diverge faster than ( |L).29| ). This is required for the consistency 
of our discussion above, where we ignored explicitly regulating the action. 

Let us now determine precisely how it scales with L in this limit. The region 



L ^> 1 corresponds to r/ ^> 1. As A s — > in the integrand in this region, (|D.31| ) 



scales linearly in r/. One can convert rj to L by using (p.25|) . We consider case i 
first. In this case one finds, 

s - ^WJk- ioel - l>>i - < d - 32 > 



From ( |3.28[ ) one finally obtains, 



Slen -► tt log L, L » 1. (D.33) 

V S K 

As it scales like log(L) our arguments above by neglecting the counter-term action in 
the large L region is thus justified. However, we note that this term does contribute 
the final result in the spin-spin correlator as it affects the sub-dominant terms that 
are denoted by ellipsis in ( |5.71|) . For example, in the mean-field approximation ( |5.26| ) 



for a second-order phase transition n = 2 one finds that the coefficient in (|D.33| ) is 
precisely 1. In order to determine the exact power of the sub- leading terms in ( |5.71|) , 



59 No where in the paper, we explicitly specified the form of the metric near the boundary, however 
we assume that it is asymptotically AdS. 
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one should take this contribution into account in addition to possible other log L 
terms that may arise from the expansion of the leading piece ( p.28|) and possible 
quantum fluctuations of the string. 

In the case ii, a similar calculation shows that 



i 



S? en ~LT^. (D.34) 

As this is always sub-dominant to the linear behavior in ( p.29| ) for a > 0, our 
discussion above by neglecting the counter-term action above is again justified. 

The same calculation should be done for dilaton-coupling in ( p,16|) . Of course 
the counter-term action has the same physical origin. It comes from the dilaton- 
coupling in the string solution that corresponds to two quarks in the Polyakov-loop 
correlator. This is given by the dilaton-contribution for two disconnected strings 
hanging from e to r/: 

S ren ,* = -^ r (2A' s (r)f(r) + 2f{r)A!' s {r) + f"(r)) <$>{r)dr. (D.35) 



This is obtained by calculating the world-sheet Ricci scalar R^ in ( p,16[) for the 



disconnected string solution. We ask whether this cancels out the UV divergence in 
the dilaton-term in ( p,18|) in the limit e — > 0. The integrand of the latter is given 



by $(r) times ( p,19| ). Near r rj the function Z in (p.20|) diverges, thus only the 
terms proportional to Z 3 in the square brackets in ( |D.19| ) survive. This is precisely 
the form that one has in ( p.35| ), thus one indeed see that the counter-term given in 
( p.35| ) does the job. As it scales exactly the same way in L as ( p.26| ), we see that, 



also the renormalized dilaton-coupling is subdominant to the area term Sq, hence 
our arguments by neglecting the renormalization above are still justified. 

D.3 The 't Hooft loop 

Here we detail the computation of the 't Hooft-Polyakov loop correlator which we 
propose to correspond to the vortex-anti vortex pair (v(x)v(0)). This is represented 
by a Dl — Dl brane pair ending on the boundary at points x and 0, and wrapping 
the time circle. Thus we want to compute the RHS of 

(v(x)v(0)) (xe~ SD \ (D.36) 

with 

^-T./^e-^ + M)', (D.37) 

where T\ is the D-string tension and we defined, 

h ab = g%d a X»d h X v , b ab = B^daX'dbX". (D.38) 
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Here g s is the string-frame metric 



ds 2 = e 2A s (r) (/-i( r ) dr 2 + dx 2 di + d xlf(r)) , A s (r) = A(r) + j^jHr), 

(D.39) 
in the BH phase. In the TG phase the metric is given by the replacement / — > 1 and 
A ->■ A in (ETTTp . 

We choose the gauge, a = x$ and r = a?i. Here X\ is one of the spatial directions 
of the spin-model on which the points x and lie. From (|D.38| ) and ( p.39|) one finds, 



/ dr \ 1 \ , , , / dr 



h 00 = e 2A °f; h u = e 2A - j 1 + l—j - I ; 6 10 = -fo i = b 1 — 1 (D.40) 
/iio = ^oi = &oo = 611 = 0, (D.41) 

where 6 is a constant given by the (r, xo) component of the B-field b = B t q. From 
(|D.37|) we find the action for the Dl — Dl pair, 



Sm = ~ £ dre-*e 2A Jf + (^j (1 + Ve~^). (D.42) 

The action only contains derivatives of r explicitly, thus the corresponding Hamil- 
tonian should be a constant of motion. Let us define the following functions for 
notational simplicity: 

p = e 4A s -2* /; ~2 = e 4A s -2* (1 + & 2 e -4A s) _ (D ^ 3) 

Then the Lagrangian in ( |D.42| ) is 



C=^Jf 2 + f 2 g 2 . (D.44) 

The canonical momentum that corresponds to r and the Hamiltonian is given by, 

dC a 2 r f 2 

p r = - r = - r y H=p r f-C = -- F . (D.45) 

dr V P + r 2 g 2 L 

We are interested in the connected Dl - Dl pair. This is given by a curve on the 
(xx,r) plane that ends on the points 60 (x, e) and (0, e) and has a turning-point at r/, 
at which f = 0. We assume that the curve is symmetric around the turning point 
which corresponds to the point (x/2, 77) on the (x\,r) plane. As the Hamiltonian is 
conserved and independent of x\ it can be computed at 77. Changing variable to r 



(>o 



e is the boundary cut-off. We shalf keep it explicit here, and it can be removed by renormalizing 



the D-string action by adding counterterms as in section D.2. 



one has, T-L = — /(r/). Then the Lagrangian is given by C = P{r)/f(rf) and one 
obtains from this, the first-order e.o.m that determines the shape of the curve: 



-i 



^P-) =± / / (r ) * Jf 2 (r)-P(r f ). (D.46) 

dr J g(r) /( r/ ) V ^ ^ " ^ ' 

One chooses the plus sign for x\ G (0, x/2) and the minus sign for X\ G (x/2,x). 
Substituting this in the Lagrangian (p.44j) one obtains the on-shell action, 

*„'/"**>=. (D .47) 



TiT 



1- 



/ 2 (r/) 



/ 2 (r) 

The distance L = \x\ between the end-points on the boundary is given by, 

dx x = 2 f f dr~ 9{r)Krf) , l , (D.48) 



P(r) I /2 (r/) 

V / 2 M 

where we used (p.45|). The on-shell D-string action, hence the vortex-anti- vortex 



correlator in ( p.36|) is given by the parametric solution of ( p.47|) and (|D.48| ) in favor 
of S D1 (L), 

In fact, this solution will only be valid for particular values of L less than some 
Lmax both for the BH and the TG background. The reason is that, the connected 
D-string solution ceases to exist beyond this value. To see this we note, first of all, 
that the integrand in ( p.48| ) is positive definite, hence L increases with increasing 



Tf. Then, in case of the BH, when rj = r^ at which the Dl and the Dl falls into the 
horizon. This corresponds to the maximum value of L for the connected D-string 
solution. At the technical level, one can see this by observing that the integrand in 
( p.48|) vanishes for 77 = r^, hence L is bounded by the value L max = L{rh). Beyond 



this point, the RHS of (p.36|) is dominated by the exchange diagram as explained in 



section 5.5. 



Let's now consider the TG solution. The distance L is again given by ( |D.48j ) 
but this time the metric functions are given by A = A$ and / = 1. One also has, 
$ = $ : 



V W) 



(D.49) 

In order to see that (p.49|) is bounded from above in the entire range rj = e to 



r 



f 



00, one divides the range into two parts (e, ri) and (ri,r/) 61 , where r\ is large 



61 This argument is first given in E3|. We warn the reader that there are typo errors in that 
reference. Here we prefer to present the argument independently. 
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enough so that we can assume that in the second range the background functions 
are approximately given by their asymptotic forms: 

$(r) « $oo r, A{r) « -i M r. (D.50) 



Let us denote the two contributions in ( |D.49| ) from the ranges (e, r\) and (?T,r/) as 



Li and L 2 respectively. One first shows that Li is bounded from above: 

u = 2 f »*£> , ■ = 2 iM f *& (r) 



fo(r) I /g(r) 1 /o(ri) 7 £ / (r) / /g(r) /gfr, j" 

< 2&M f rfrl^l , X = 2^L(r 1 ). (D.51) 

/o(ri) A / (r) / f§(r) _ 1 /o(ri) 

V /o('-i) 

In the second line we used the fact that fo(r) is a monotonically decreasing function 62 . 
Therefore the only possible divergence in L can come from L 2 : 

l 2 = 2 r *m - i « ^±? r~ (r/ " n) ^ (^ 

A /o(r) / jg(r) , $00 7o e»-l' 

V /o(-/) 

where we used that, by assumption the background functions are given by the asymp- 
totic forms ( p.50[ ) in this range of r. We see that this is bounded from above. In the 



limit rf — y 00 one finds, 

lim L 2 (r f ) « ?rV J TF . (D.53) 

In fact, with little more effort, one can show that the RHS of (p.53|) is the upper 
bound on L 2 . Thus L is bounded from above and there is a maximum value L max 
that is reached at some point rj = r max . It is clear from the calculation above that 
this point is independent of temperature in the TG phase. Beyond this point, the 
connected D-string solution ceases to exist and the vortex correlator is determined 



by the exchange diagram, cf. section |5J) 



As a last comment, we observe that the calculation above could be applied 
directly to the BH case, just by replacing the functions / and g by / and g. The 
crucial point about the monotonicity of the function / is guaranteed just like in the 
footnote below, given that / is also monotonically decreasing. Thus, also in the 
BH geometry, one has a point L' max (r h ) above which the connected diagram does 



(i2r 



This is clear from the definition ( D.49J ). The exponent is 4Ao(r) + 2$o(r)/3 where Aq is 



the Einstein frame scale factor. By assumption, Aq is a monotonically decreasing and $0 is a 
monotonically increasing function. We also know that the combination AA s $(r) = 4Ao(r)+8$o( r )/3 
is monotonically decreasing and having an asymptotic minimum at r = 00. Thus 4A (r) + 2$o( r )/3 
should be monotonically decreasing, with no minimum. 
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not exist. It is an interesting question, whether this point is before or beyond the 
horizon. Namely, one can ask the question whether L' max (rh) < L(rh) or not. If so, 
then the connected D-string solution would cease to exist even before it falls into the 
horizon! The answer will be determined by the precise background functions of the 



holographic model, however, this point does not change our arguments in section [575 
that only depends of existence of some L max . 

E. Spectrum of bulk fluctuations 

E.l Graviton and dilaton 

We first consider the spectrum of the dilaton and the graviton fluctuations. We shall 
carry out the calculation with Minkowski signature for convenience. The results 
can easily be translated in Euclidean time by analytic continuation. There are two 
independent modes with spin-0 and spin-2. The first one is given by a mixture 
of the isotropic fluctuations of the metric components h\\ = h 22 = ■ ■ ■ = hd-i,d-i 
and the dilaton fluctuations 5$. We denote this mode as ho- The second one is 
is just the transverse traceless shear fluctuations that will be denoted by h%. The 
equations of motion for these modes on a generic BH background (|3.6|) , are obtained 
by decomposing h Q2 (r,x) — e~ ipiXi+tuJX ho }2 (r). They can be found for example in 
[L5|] , J52] . In the Einstein frame one obtains 63 



K + K ((<* - IK + f j + 2§) + »o {^f- - ! ~) = 0, (E.1) 
hi + h' 2 (td - l)A' + £) + h 2 (^jr^-) = o. (E.2) 

where the function X is given by X(r) = 2/((d— l)v«) <$>'/A' and it asymptotes to a 
constant X — >■ —1/y/d in the IR region r — > oo. The analogous fluctuation equations 
on the thermal gas background us given by setting / = 1 in these equations. 

From ([E.1D and ( E.2|) we see that the spin-0 and spin-2 modes become degenerate 
in this far IR region, hence it suffices to consider only the latter. As argued before and 
shown in |26| in the limit r^ — > oo the function / approaches to 1 and A approaches 
to the scale factor of the TG solution ( |3.5| ) A — > A . Therefore, in this regime of 
interest we want to solve, 

h'i + ti 2 (d - 1)A' + h 2 (u 2 - \p\ 2 ) = 0. (E.3) 

One can easily transform this equation to a Schrodinger form by h 2 = h 2 exp(— (d — 
l)A/2): 

-h>> + V s (r)h 2 = (u 2 - \p\ 2 )h 2 , V S = d -^K + ^-^<- (E-4) 

63 Strictly speaking these equations are correct only when either of \p\ 2 = p l p % or w 2 vanish. 
Otherwise there may be some more complicated mixing terms. We will keep this combination and 
in the end of the computation we will be interested either of these two cases. 
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The asymptotics of the function A ( 3.17 ) imply that the Schrodinger potential asymp- 
totes to a constant in the far IR: 

V s ^ml + 0(e- K ^ r ), m\ = \ (E.5) 



and for the form of the subleading corrections we refer to 

On the black-hole background the fluctuation equation ( |E.2|) is solved by impos- 
ing normalizability near the boundary and incoming boundary condition (or normal- 
izibility in the Euclidean signature) at the horizon. Consider first the case |p| = 0. 
Then one gets a discrete spectrum of a; 2 on the BH. As r^ is taken to oo, i.e. near the 
phase transition region T — >■ T c the spectrum becomes nearly continuous with the 
lowest lying state determined by the asymptotic constant in ( E.5|) . Consequently, at 
finite but large r^ the spectrum of states are of the form, 



Uo = m = ^-^, Wl = m0 + O(e- K — rh ), etc. (E.6) 

The constant mo is the same as the one that appears in the derivative of the dilaton 



( |4.4ft . This is not a coincidence but required for the consistency of the theory. What 
we learned is that, the gravitational fluctuations on the linear-dilaton background is 
always gapped with a gap factor m , i.e. 

u 2 > m 2 , for N = N = 1. (E.7) 

In the actual calculation of the two-point function in section |5.4.1| and [5.4.2 



we need the Euclidean spectrum with compact time. One finds that the exchange 
diagram is of the form exp(— mL) where m is always bounded as in ( E.7|) from below. 

E.2 B-fleld 

We consider the fluctuations around B r0 = const that we denote as Sip. The spectrum 
of fluctuations are obtained from the equation of motion for the B-field d * dB = 0. 
We decompose 5ijj(r,xo) = e ~ l Pi x i+ luJX S^f(r). On simply obtains u 2 = 0, hence 
the spectrum of fluctuations from the point of view of d — 1 dimensions are also 
massless. This means that in the exchange diagrams of section ( |5.4.1| ) and (|5.4.2|) , 



the contribution from the lowest CT modes is massless, in accord with existence of 
the Goldstone mode in the superfluid phase. 

E.3 Tachyon 



The tachyon action is [16| (in the string frame) 



At ~ J Vale- 2 * Urd^duT - ^T 2 ^j , (E.8) 
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where the metric and the dilaton reads, 

ds 2 = e 2A ° {r) {f~\r)dr 2 + dx 2 d _ x + dx 2 f(r)) 



$ = ruor, 



(E.9) 



and mo is defined in ( |4.4|) . We shall carry out the calculation on the black- hole, and 
the thermal gas result will be obtained simply by setting / = 1. 

We fluctuate T = (T) + T(r)e~ tpZx * +UJX0 in the action. We note that these fluc- 
tuations do not mix with the dilaton fluctuations for (T) = which is indeed what 
we assume throughout the paper: the only non-trivial profiles in the background are 
the metric and the dilaton. It is straightforward to obtain the fluctuation equation 
from ( |E.8|) and (|E.9| ). As we are interested in the spectrum of u we set p l = 0: 



T" + 



(d - 1)4 -2m Q + f - 



r 



-,2A S 



4 e 2 ^ 8 ur , 

5- + ?" " 



(E.10) 



This can be transformed into a Schodinger form by the change of variable T 
exp(-(d - I) A/2 + m r - \ log f)f with the result, 



LU 2 ^ 



T" + V T T = —T 



P 



Yr 



d-1 



A"- 



v — - A s -m + ±-) 2 + -( J —- J —) 



4e 



2A S 



2f J 2 , f p, £ 2 f 

(E.ll) 
Let us first consider the simpler case of the thermal gas background that is obtained 
from ([E.ll]) by setting / = 1: 



Vn 



(l 1 A'> + (^lA' a -m ) 2 -^e 2A ^ 



TG 



2 s v 2 
In the TG phase in the IR A s vanishes as ( p.8| ) 

A s (r) -> a ie - Kmor 
with some positive coefficient at. The Schrodinger potential then becomes, 



m r 



TG 



- ■■■ "i ( '" Kmor ( ^(d - 1)(« + 2)«mg - ^ ) + 0(e- 2 '" 1 "^ 



(E.12) 



(E.13) 



(E.14) 



For a moment let us consider the pure linear-dilaton geometry. In this case the 
exponential correction term in ([E.14 ) is absent g^ = 0, and one obtains the exact 
answer as, 



Vn 



m n - — 



LD 



l-d 



X-d 



bosonic 
fermionic 



(E.15) 



where we used the no-anomaly condition ( |4.5|) . The tachyon in the fermionic case 
comes from the ground state of the NS fermions and has the mass m 2 ^ = —2/£ 2 . 
We re-derived the well-known result that the "tachyon" in dimensions 2 or less is 



-93- 



actually a stable mode (recalling that our total number of dimensions is d+ 1). This 
is of course a consistency check. 

Coming back the issue of the spectrum, the result (|E.14j) indicates that the fluc- 



tuations in the deep interior of the thermal gas geometry, in the vicinity of the phase 
transition T — y T c 64 are tachyonic. Luckily we do not need this lowest mode in the 
calculation of the two-point function in section |5.4.2| because the entire propagation 



is governed by modes with winding mode w = 1 which are non-tachyonic. 

However, we needed this mode in the calculation of the two-point function in 



the black-hole phase, cf. section |5.4.2| . Now let us inspect the spectrum of tachyon 
fluctuations on the black-hole. This is determined by the equation (|E.10|) . The 
blackness function near the horizon behaves as 

f->4irT(r h -r), r -> r h . (E.16) 

Then the fluctuation equation becomes the standard form, 

T"-T\r h -r)- l + - T « 0, r -> r h , (E.17) 

ir h - r) 1 

where Co = lo/4ttT. This can be solved by changing the variable as r^ — r = exp(— u) 
and the solution near the horizon becomes, 

T -> T + e^ u + T_e~^ u . (E.18) 

The incoming one corresponds to T + = 0. To inspect the issue of the tachyon, 
we can change to the Euclidean metric by u> — y —iuo and the incoming solution of 
the Minkowskian BH corresponds to the normalizable solution of the Euclidean one. 
This means that the Euclidean spectrum is always discrete and bounded from below. 
However we still have to see whether there is a negative mode in the limit r^ — y oo 
(T — y T c ). We recall that in this limit the BH geometry asymptotes to the TG 
geometry. In particular 

/(r) -y 1, for all r < r h . (E.19) 



Then, for any r < rh the corresponding Schodinger potential is given by ( E.14|) which 



becomes negative in the far r region for r < r^. Near r^ it becomes positive again and 
finally it diverges as (rh — r)~ 2 as r — y rh- Then existence of a negative discrete mode 
crucially depends on whether the approach to the negative minimum that is given 
by ( |E.15| ) is from above or below. This is determined by the sign of the coefficient 
of the exponential term in ( [E.14| ). Recalling that g^ > 0, we observe that the sign is 



always positive for the interesting case of k = 2 which corresponds to a second order 
transition both for d — 1 = 2 and d — 1 = 3. For a third or higher order transition it 



(>4 



This is the only regime the world-sheet CFT becomes linear dilaton 
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is negative both for d — 1 = 2, 3. We conclude that in the cases of interest, although 
there is a negative minimum in the tachyon potential, the approach to this minimum 
is from above and the potential can always be arranged (by choosing the form of the 
next-to-subleading terms in the dilaton potential ( |3.15|) ) so that there is no negative 
discrete mode in the spectrum. The same cannot be said for fluctuations on the 
thermal gas, as explained above. 

It is a reasonable question to ask whether the tachyon of the linear-dilaton CFT 
(on the thermal gas) can be extrapolated to the UV theory. To answer this question 
one has to study the tachyon potential in the UV. This can only be done in an 
heuristic way. The reason is that, in the UV we do not have an exact CFT description 
unlike in the IR and the a' corrections would renormalize the following discussion. 
Nevertheless, let us pretend that there are no a' corrections in order to see what 
possible behavior can arise. In this paper we have not specified the UV geometry, 
but in fact we always tacitly assume that the UV geometry is AdS. In ]I5|] we found 



analytic kink solutions that flow from the UV in the AdS and linear-dilaton in the 
IR. In the case of AdS the metric scale factor is A — > — log r/£+- ■ ■ . Then we obtain, 

,d-l (d-lf U 2 ,l _ , 

^«(— + L - r J "- 1 r)pj, r->0. (E.20) 

This will be bounded only when the term in the bracket is positive. In the case of 
d — 1 = 3 this gives the condition £ 2 /£ 2 S < 15/32. As mentioned above, this result is 
supposed to be corrected by a' corrections. However it is reasonable to expect that 
there will always be an upper bound on £/£ s by demanding that there is no d — 1 
dimensional tachyon in the spectrum in the UV. On the other hand, the simplest way 
to achieve this is to demand that there is no tachyon to start with i.e. the spectrum 
in the d+ 1 dimensional theory is non-tachyonic and the tachyon of the linear-dilaton 
theory only arises in the IR effective theory. 
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